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The behavior of photons in the presence of Lorentz and CPT violation is studied. Allowing for 
operators of arbitrary mass dimension, we classify all gauge-invariant Lorentz- and CPT-violating 
terms in the quadratic Lagrange density associated with the effective photon propagator. The co- 
variant dispersion relation is obtained, and conditions for birefringence are discussed. We provide 
a complete characterization of the coefficients for Lorentz violation for all mass dimensions via a 
decomposition using spin-weighted spherical harmonics. The resulting nine independent sets of 
spherical coefficients control birefringence, dispersion, and anisotropy in the photon propagator. 
We discuss the restriction of the general theory to various special models, including among others 
the minimal Standard-Model Extension, the isotropic limit, the case of vacuum propagation, the 
nonbirefringent limit, and the vacuum-orthogonal model. The transformation of the spherical coef- 
ficients for Lorentz violation between the laboratory frame and the standard Sun-centered frame is 
provided. We apply the results to various astrophysical observations and laboratory experiments. 
Astrophysical searches of relevance include studies of birefringence and of dispersion. We use po- 
larimetric and dispersive data from gamma-ray bursts to set constraints on coefficients for Lorentz 
violation involving operators of dimensions four through nine, and we describe the mixing of polar- 
izations induced by Lorentz and CPT violation in the cosmic-microwave background. Laboratory 
searches of interest include cavity experiments. We present the general theory for searches with 
cavities, derive the experiment-dependent factors for coefficients in the vacuum-orthogonal model, 
and predict the corresponding frequency shift for a circular-cylindrical cavity. 



I. INTRODUCTION 

The properties of electromagnetic radiation have 
proved a fertile testing ground for relativity since its in- 
ception over a century ago. Tests such as the classic 
Michelson-Morley, Kennedy-Thorndike, and Ives-Stilwell 
experiments [1-4] were key in establishing Lorentz in- 
variance, the foundational symmetry of relativity. The 
proposal that minuscule deviations from Lorentz sym- 
metry could emerge from an underlying unified theory 
[5] has rekindled interest in sensitive relativity tests, and 
the past decade has seen a broad variety of searches for 
Lorentz violation at impressive sensitivities [6]. 

Violations of Lorentz symmetry at attainable energies 
are described by the Standard-Model Extension (SME) 
[7, 8]. The SME is a comprehensive effective field the- 
ory that characterizes general Lorentz and CPT viola- 
tions. It contains both General Relativity and the Stan- 
dard Model, and so it is a realistic theory that can be 
applied to analyze observational and experimental data. 
A Lorentz-violating term in the Lagrange density of the 
SME is an observer scalar density formed by contracting 
a Lorentz-violating operator with a coefficient that acts 
to govern the term. The operator can be characterized in 
part by its mass dimension d, which determines the di- 
mensionality of the coefficient and can be used as a naive 
guide to the size of the associated effects [9] . 

The focus of the present work is Lorentz and CPT vio- 
lation involving photons. Numerous searches for Lorentz 
violation in electrodynamics have been performed, yield- 
ing some of the best existing constraints on SME coef- 
ficients. One major class of tests consists of laboratory 
searches involving electromagnetic resonators or interfer- 



ometers [10-18], which can be viewed as contemporary 
versions of the classic tests of relativity. Another ma- 
jor class of tests consists of astrophysical observations 
searching for tiny defects in the propagation of light that 
has traveled over cosmological distances [7, 11, 19-23]. A 
variety of other analyses involving photons lead to con- 
straints on SME coefficients [24] . There is also a substan- 
tial literature on various topics in the photon sector of the 
SME, including renormalization [25], photon interactions 
[26], vacuum Cerenkov radiation [27], the Chern-Simons 
term [28], electromagnetostatics [29], and related phe- 
nomena involving photons in other contexts [30]. Out- 
side the photon sector, the SME serves as the theoretical 
underpinning for studies of Lorentz symmetry involving 
electrons [31, 32], protons and neutrons [33-35], mesons 
[36], muons [37], neutrinos [38], the Higgs [39], and grav- 
ity [8, 40]. 

In this paper, we extend the existing treatment of 
Lorentz violation in electrodynamics to include opera- 
tors of arbitrary mass dimension d. This further develops 
and consolidates previous systematic studies for opera- 
tors of renormalizable dimension [7, 11], and it incorpo- 
rates many phenomenological models for Lorentz viola- 
tion. For definiteness, we focus on an action having the 
usual U(l) gauge symmetry and invariance under space- 
time translations, so that charge, energy, and momentum 
are conserved. If the Lorentz violation is spontaneous, 
then the SME coefficients originate as expectation values 
of operators in an underlying theory, and the requirement 
of invariance under spacctimc translations corresponds to 
disregarding soliton solutions and any massive or Nambu- 
Goldstone (NG) modes [41]. In a more complete treat- 
ment including gravity, the NG modes may play the role 
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of the graviton [42]. Alternatively, the NG modes may 
be interpreted as the photon for Einstein-Maxwell the- 
ory embedded in a Lorentz- violating vector model called 
bumblebee electrodynamics [8, 43]. The approach dis- 
cussed here can readily be adapted to these and other 
scenarios, including applications to the photon sector of 
the SME in the context of the various topics mentioned 
above. Our methods are also relevant for other sectors 
of the SME [44]. 

The motivation for this work comes in part from cur- 
rent doctrine, which regards the combination of General 
Relativity and the Standard Model as the low-energy 
limit of a unified quantum gravity theory that holds 
sway at the Planck scale, Mpianck ~ 10^^ GeV. Expe- 
rience teaches us to expect a smooth transition from the 
known low-energy physics to the underlying theory, so 
it is plausible to interpret the low-energy action as the 
zeroth-order term in a series approximating the underly- 
ing theory. Dimensional analysis suggests that operators 
with larger d correspond to higher-order corrections. For 
physics involving violations of Lorentz and CPT symme- 
try, the complete series is given by the SME action, while 
the leading corrections form the action of the minimal 
SME. Lorentz- violating operators of larger d are there- 
fore likely to be especially relevant in searches involving 
very high energies and in theoretical studies of founda- 
tional properties such as causality and stability [45] . Un- 
der suitable circumstances, nonrenormalizable operators 
may even dominate the physics. For example, the action 
of noncommutative quantum electrodynamics [46] incor- 
porates Lorentz-violating effects associated with a non- 
trivial commutator for the spacetime coordinates. When 
the action is expressed in terms of conventional photon 
fields, a subset of the SME emerges in which the lowest- 
order Lorentz-violating operators have mass dimension 
six [47]. Similarly, operators of larger d dominate in 
Lorentz-violating theories with supersymmetry [48] . 

A comprehensive investigation of all Lorentz-violating 
operators with arbitrary mass dimensions is a challenging 
task. Here, we concentrate on terms in the action that 
are quadratic in the photon field and therefore con- 
tribute to the propagator, which in practice is the quan- 
tity of immediate interest in many searches for Lorentz 
violation. Our basic approach consists of constructing 
the quadratic action and developing a scheme to classify 
the operators. Rotations are a prominent subgroup of the 
Lorentz group, and a spherical decomposition can be per- 
formed on any Lorentz-violating operator. We use this 
fact to classify all Lorentz-violating terms in the action 
using nine sets of coefficients for Lorentz violation. The 
classification scheme is well matched to the description of 
physical Lorentz-violating effects in the photon propaga- 
tor, including birefringence, dispersion, and anisotropy. 

With this classification scheme taming the infinite 
number of operators, specific analyses of observational 
and experimental data become feasible. We study a va- 
riety of methods for seeking Lorentz violation using pre- 
dictions from the quadratic action. The sharpest tests 



involve astrophysical birefringence, which involves propa- 
gation in the vacuum. Some Lorentz-violating operators 
produce no vacuum birefringence at leading order but 
nonetheless cause dispersion in the vacuum, and these 
can also be studied using astrophysical observations. In 
addition, there are many other Lorentz-violating opera- 
tors that are undetectable at leading order via astrophys- 
ical observations and hence are best sought instead in 
laboratory experiments. The analyses in this work yield 
several first measurements of coefficients for Lorentz vio- 
lation, and numerous interesting arenas emerge for future 
exploration. 

The structure of this paper is as follows. The ba- 
sic theory is discussed in Sec. II, which contains five 
subsections. The construction and counting of Lorentz- 
violating operators of all mass dimensions is presented 
in Sec. II A, while the Lagrange density and constitutive 
relations are obtained in Sec. II B. We derive; the covari- 
ant dispersion relation in Sec. II C, discuss the physics of 
birefringence in Sec. II D, and offer general comments on 
Lorentz-violating effects in Sec. II E. 

The spherical decomposition of the coefficients for 
Lorentz violation in terms of spin- weighted spherical har- 
monics is performed in Sec. III. We consider various 
special limits in Sec. IV, including the minimal SME, 
isotropic models, the vacuum limit, nondispersive non- 
birefringent 'camouflage' models, the vacuum-orthogonal 
case, and some limits providing connections to other for- 
malisms. The rotation properties of the spherical coeffi- 
cients for Lorentz violation are discussed in Sec. V. Other 
key properties of the spin-weighted spherical harmonics 
are summarized in Appendix A. 

The remainder of the paper applies the results to ob- 
servations and experiments. Astrophysical observations 
are studied in Sec. VI. Vacuum dispersive effects are 
discussed in Sec. VIA, where new measurements and a 
summary of existing constraints are obtained. Vacuum 
birefringence is considered in Sec. VI B, which contains 
three subsections. Some basic theory for vacuum birefrin- 
gence is presented in Sec. VI Bl. We apply it to point 
sources in Sec. VI B 2, using polarimetry from gamma- 
ray bursts to obtain new measurements, and also to the 
cosmic microwave background in Sec. VI B 3. In both 
cases, we tabulate some existing sensitivities. 

Laboratory experiments are discussed in Sec. VII. We 
construct a general theory for resonant-cavity tests in 
Sec. VII A and apply it in Sec. VII B to derive cavity 
factors for nonbirefringent operators and the form of the 
fractional frcriuency shift predicted by Lorentz violation. 
Explicit values of some cavity factors for a circular cylin- 
drical cavity are calculated in Sec. VII C. Section VIII 
provides a summary and discussion of the results in the 
paper, including tables compiling essential properties of 
the spherical coefficients and various limiting cases. Un- 
less otherwise stated, this paper follows the notation and 
conventions of Ref. [11]. 
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II. BASIC THEORY 

This section discusses the theory and basic features 
of the quadratic action for electrodynamics allowing for 
arbitrary Lorentz and CPT violation. We begin with a 
discussion of the SME procedure for constructing the La- 
grange density associated with the effective photon prop- 
agator. Attention is focused primarily on the case with 
conventional U(l) gauge invariancc and translational in- 
variance. The resulting effective field theory conserves 
charge, energy, and momentum. It represents an explicit 
presentation of all Lorentz-violating operators for pho- 
ton propagation that are consistent with observer covari- 
ance. Following the construction of the theory, we ex- 
tract a complete set of coefficients for Lorentz and CPT 
violation and discuss some of their basic properties. A 
technique for deriving covariant dispersion relations is 
presented and used to obtain a general covariant disper- 
sion relation for the photon in the presence of arbitrary 
Lorentz and CPT violation. The issue of conditions for 
birefringence is considered, and we offer some remarks 
about the connections between birefringence, metricity, 
and electromagnetic duality. 



A. Construction 

A low-energy Lorentz-violating theory that is both co- 
ordinate independent and consistent with current obser- 
vational data can be written as a Lagrange density con- 
taining sums of standard polynomial tensor operators 
contracted with coefficients for Lorentz violation [7] . The 
coefficients may be viewed as background fields inducing 
Lorentz violation, and they can correspond to vacuum ex- 
pectation values of fundamental tensor fields. Applying 
this general idea to source-free linear electrodynamics, 
we arrive at an action S that is a quadratic functional 
of the photon field and its derivatives. The action S 
can then be expanded in a sum of terms S'(d) of the form 

S{d) — J df^X K,"^-^^'^^ "'*j4q,j9„3 . . . 9q;^j4q,2, (1) 

where d is the dimension of the tensor operator. Each 
term S(^ii^ violates CPT if d is odd or preserves CPT if d 
is even. The coefficients /(;^^i^"2a3 - "d j^g^yg mass dimen- 
sion 4 — d. In general, they can be dynamical and can 
depend on spacetime position. We can ensure invariance 
of S under spacetime translations, and hence obtain the 
usual conservation of energy and momentum, by restrict- 
ing attention to the case of constant coefficients. Con- 
stant coefficients may arise naturally, but they may also 
represent the dominant components of dynamical back- 
ground fields or an averaged effect. 

Size estimates for the coefficients for Lorentz violation 
depend on the details of their origins. Since the effects 
are expected to be small, it is natural to suppose the 
Lorentz-violating operators are suppressed by some large 



mass. The intimate connection between Lorentz sym- 
metry and gravity suggests the relevant scale is set by 
Planck-scale physics and therefore by the Planck mass 
.^Planck- Various scenarios can be imagined, although 
the absence of a satisfactory underlying theory combin- 
ing gravity and quantum physics, and hence the lack of 
specifics concerning possible Lorentz violations, makes 
such scenarios a matter of surmise and likely naive. For 
example, one simple estimate has coefficients varying as 
/C((;) ~ C-^punck' where ( is of order 1. This means oper- 
ators of renormalizable dimension < 4 are unsuppressed 
relative to conventional physics. Observations then im- 
ply only the operators d > 5 are experimentally viable, 
so the dominant new physics is controlled by nonrenor- 
malizable terms. Another class of scenarios has Lorentz- 
violating effects suppressed by a factor involving some 
power of the ratio m/Mpianckj where m is an appropriate 
low-energy scale. In these cases, the Lorentz violation 
may be related to one or more of the known hierarchies 
in nature [9]. For example, taking m ^ 10^ GcV as the 
electroweak scale gives a dimensionless suppression factor 
of some power of 10~^^ for Lorentz-violating physics 
compared to conventional effects. 

While a general study of all possible operators of the 
form (1) would be of interest, it would be rather cum- 
bersome and introduce various features in addition to 
Lorentz violations, thereby complicating both theoreti- 
cal and experimental considerations. The possibilities 
are simplified somewhat by restricting attention to oper- 
ators that maintain the conservation of energy, momen- 
tum, and electric charge. This implies focusing on the 
case of constant coefficients and requiring U(l) gauge in- 
variance. The latter imposes certain symmetries on the 
Lorentz-violating operators, thereby reducing the total 
number of independent coefficients. 

The first step in imposing these symmetries is to iden- 
tify properties of the coefficients K^'^^^^as-.-aa ^^i^^^ follow 

from the intrinsic structure of S. One property of the co- 
efficients is total symmetry in the d—2 indices {as . . . aa}. 

Another can be displayed by integrating Eq. (1) by parts 
d—2 times. This reveals that the only coefficients con- 
tributing to CPT-odd terms are antisymmetric in the 
first two indices of jC'^^^^"^---"^ ^ while the only ones con- 
tributing to CPT-even terms are symmetric. The same 
conclusion follows by considering the contributions to the 
variation of the action: 

SS^a) = J d^x /C«J«^«3...a. (5^^ . . . 

-I- surface terms, (2) 

where the minus and plus signs apply for CPT-odd and 
CPT-even terms, respectively, and where the brackets 
[ ]± indicate symmctrization and antisymmetrization. 

With these intrinsic symmetries understood, term-by- 
term gauge invariance can be imposed and the resulting 
additional symmetries of S(^^) identified. The usual U(l) 
gauge invariance, representing symmetry of the action 
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FIG. 1: Representation decomposition for the term S(^cl) using 
Young tableaux. 
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g'J{d) — I ^ ^ /^(^) 

X Adas ■ ■ ■ ^ctd (^[ai^a2l± + (3) 

vanish for an arbitrary scalar function A. Direct inves- 
tigation of this equation is possible but awkward. To 
identify the additional symmetries and hence the coef- 
ficients of interest, it is more practical first to perform 
representation decompositions of the associated Lorentz- 
violating operators. The intrinsic symmetries drastically 
limit the number of representations that can appear. The 
total symmetry in the last d-2 indices of /C^^i^"2«3-"<i 

implies that all representations that are antisymmetric 
in any pair of these indices are absent. We can therefore 
construct the relevant irreducible tensors from a prod- 
uct of symmetric representations. It turns out that this 
limits the possibilities to only five representations. In 
terms of Young tableaux, these five representations are 
displayed in Fig. 1. 

For CPT-odd coefficients, the antisymmetry condition 
on the first two indices and the gauge variation (3) imply 
that tableau (a) is irrelevant, since it is symmetric in the 
indices {aia2}. Also, representation (d) is symmetric 
under the simultaneous permutation of 1 <-> 2 and 3 4 
and so fails to contribute. The antisymmetry in {1, 2, 3} 
of representation (e) directly implies gauge invariance, so 
it satisfies our restrictions. The remaining two represen- 
tations (b) and (c) lead to nonvanishing SgS{d) and are 
therefore gauge violating. We conclude that all gauge- 
invariant CPT-odd operators are associated with coeffi- 
cients /C"j^j"^"^'""'' belonging to representation (e). These 
are antisymmetric in the first three indices and symmet- 
ric in the last d — 3. Since dimension d < 1 operators 
are absent in a linear theory, we take d > 3 for CPT-odd 



operators in what follows. 

For CPT-even coefficients, the antisymmetry of 
tableau (e) in {1,2} implies it cannot contribute to Eq. 
(2) and therefore is absent in the theory. The antisymme- 
try of tableau (d) in {1,3} and {2,4} leads to SgS{d) = 0, 
so this representation is gauge invariant and satisfies our 
restrictions. The tableau (a) is totally symmetric and 
leads to SgS{d) ^ 0, implying gauge violation. Similarly, 
tableaux (b) and (c) are also gauge violating. We thus 
find that all gauge-invariant CPT-even operators have 
coefficients /^«^i«2«3. .Q:<i belonging to representation (d). 
Note that operators with d = 2 are gauge-violating, so in 
what follows we take d > 4 for the CPT-even sector. 

To simplify handling and to provide a convenient 
match to the usual coefficients for Lorentz violation in 
the minimal SME, it is convenient to introduce further 
definitions. For the CPT-odd case, the dual coefficients 



(fc 



AF/ 



/ii/pQia2...a<j_3 



(4) 



provide a generalization of the usual coefficients {kj^p)^, 
in the minimal SME. The synuiietries of tableau 
(e) translate into total symmetry of the coefficients 

(^af)« ^' 1*^^^ d — 3 indices, along with 

the trace condition {k^^p)ai '"^ ^' = 0. Counting the 
number of independent components in representation (e) 
using standard group-theory techniques [49] yields for di- 
mension d the result 



N^^l = ^id+l){d-l){d-2). 



(5) 



This number can also be obtained by noting that sym- 
metry in the last d — 3 indices yields 4{d— 2){d — l)d/3\ 
components while the trace condition given above pro- 
vides {d — 3){d — 2)(d — l)/3! constraints, and taking the 
difference yields N^^p. Note that the number of these co- 
efficients for CPT-odd Lorentz violation grows rapidly as 
the cube of d: the usual 4 for = 3, then 36 coefficients 
for d = 5, 120 for d = 7, etc. 

For the CPT-even case, it suffices to define 



^j^{d)-jK,\lj,i/ai...a^d-4) 



id) 



(6) 



to obtain coefficients that mimic the definition of 
(kp)'^^'^" in the minimal SME. The first four indices of 
(/c^')'''^^'^"!'--"'''-'') have the symmetries of the Riemann 
tensor, and there is total symmetry in the remaining d— 4 
indices. Also, one can show that antisymmetrization of 
(fci?''')'^'^''^"^ on any three indices produces zero. 

In this case, counting the number of independent com- 
ponents for dimension d gives 



id) 



{d + l)d{d-3). 



(7) 



This counting includes the total trace term, which is 

Lorentz invariant and represents a scaling factor. Again, 
note that the number of these coefficients for CPT-even 
Lorentz violation grows as the cube of d: the usual 20 
(19 plus a Lorentz-invariant trace) for d = 4, then 126 
coefficients for d = 6, 360 for d = 8, etc. 
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B. Lagrange density and constitutive relations 

The construction outlined in the previous subsection 
leads to a general gauge-invariant Lagrange density that 
can be written in a form similar to the photon sector of 
the minimal SME: 

£ = -^F^.F^^-' + ^e^^^-Ax{kAF).F^. 

- \F,x{kFr^'"'F^,, (8) 

where the differential operators kAF and kp involve CPT- 
odd and CPT-cvcn violations, respectively. These oper- 
ators arc given by the expansions 

{kAF)^ = J2 5ai •••5a(^_3), (9) 

d=odd 

d=even 

(10) 

where the sums range over values d > 3. The co- 
efficients (fc^^)^ ''^ are defined in Eq. (4), and 
they have symmetry properties yielding N^^p indepen- 
dent components as given in Eq. (5). The coefficients 
(^(d)>)KAM<'ai...a(d_4) are defined in Eq. (6) and have N^J'^ 
independent components according to Eq. (7). The usual 

minimal SME terms are kAF = k[^p and kp = fc^"* , where 
in the latter the overall trace is removed to leave 19 in- 
dependent coefficients. 

In principle, obtaining and interpreting equations of 
motion for a Lagrange density of the form (8) is problem- 
atic due to the infinite sum, whose action cannot be var- 
ied in the usual way, and also due to Ostrogradski insta- 
bilities [50]. However, these issues can be circumvented 
by noting that Eq. (8) represents the low-energy limit of a 
more fundamental theory, with each successive term rep- 
resenting a perturbation on preceding terms. Truncating 
the Slims at any definite value of d and restricting at- 
tention to perturbative effects therefore can be expected 
to provide a good approximation to the low-energy be- 
havior. Only at extreme energies can qualitatively new 
effects and late terms in the sum play an important role. 
At these energies, the theory must converge to the un- 
derlying fundamental physics, which presumably is free 
of these issues. Adopting this truncation, we obtain equa- 
tions of motion given by 

««r?^^a, + (fcAF).e'*''"^ + {kFr-'''^d^)F^i} = 0. (11) 

Note the explicit gauge invariance of these equations. 

In the minimal SME, the coefficients fc^^ and kp^ are 
known to produce photon behavior analogous to that 
of conventional electrodynamics in anisotropic and gy- 
rotropic materials [7, 11]. This analogy can be extended 
to the present situation involving Lorentz- violating oper- 
ators of arbitrary dimension. The first step is to define a 
field tensor 

G^'" = - 2e''''^^{kAF)aA0 + CkFr^'^^F^is, (12) 



in terms of which the equations of motion become 

d^G'"' = 0. (13) 

Note that the latter equation is gauge invariant, even 

though the definition of G^'^ depends on the choice of 
gauge and is unique only up to a term of the form 
€'^''°'l^{kAF)adfjK for an arbitrary scalar function A. 

For conventional electrodynamics in macroscopic me- 
dia, a constitutive 4-tensor % is typically introduced 
that maps the 2-form field strength F to the macro- 
scopic 2-tensor field strength G, via G^"' = x^"''" Fp<j- 
In the present context, we can reformulate the situa- 
tion in terms of a set of unconventional constitutive re- 
lations, which may explicitly depend on the choice of 
gauge. However, we must generalize the usual notion of 
a constitutive 4-tensor to encompass more general oper- 
ator constitutive tensors. We now require an operator 
4-tensor x'^^p'^ and an operator 3- tensor X^'^p, defined 
by 

X^'^P = e^''f"'{kAF)a- (14) 

Note that the 3-tensor controls CPT violation. The ef- 
fective macroscopic field strength G defined in Eq. (12) 
is then given by 

G"" = x'"'P''Fpa + 2X'"'PAp. (15) 

As in conventional electrodynamics, the new constitutive 
relations remain linear. However, unlike electrodynam- 
ics in linear media, the relations (15) inherit a nonlocal 
aspect due to their differential nature. 

Decomposing G^"^ into an effective vector displacement 
field D and an effective pscudovcctor magnetic field H in 
the usual way, the equations of motion (13) take the same 
form as the familiar source-free inhomogeneous Maxwell 
equations, 

V • D = 0, 
V xH- doD = 0, (16) 

where 

D = E + 2kAF X A + koE ■ E + kdb ■ B, 
H = B- 2{kAF)QA + 2feAF^o + khb ■ B + khe ■ E 

(17) 

with 

{KDEy'' = -2(fcF)0^■0^ 

{KHSy'' = l(fcF)'™''' (18) 
{KDBy'' = -{fiHE)''' = (fcF)°^''"e*='™. 

The latter equations are operator generalizations of the 
SO (3) decomposition of the coefficients kp^ into 3x3 
matrices introduced in Ref. [11]. 
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Using the fields D and H, the Lagrange density may 
be written as 

C = -lF^M'"' = UE-D-B-H), (19) 

which also parallels conventional electrodynamics in 
macroscopic media. We remark that the analogy with 
electrodynamics breaks down when attempting to con- 
struct a conserved energy-momentum tensor. Standard 
techniques can be used to build a conserved tensor 
that reduces to the conventional symmetrized energy- 
momentum tensor in the limit of vanishing Lorentz vi- 
olations. However, the resulting tensor takes an uncon- 
ventional form in terms of F^^, and G'*''. One possibility 
is the tensor 

= -G"^J>7 - + l{di3A^ - A^d0)G"\ (20) 

The first two terms parallel those in conventional elec- 
trodynamics, but the addition of the last term is neces- 
sary for energy-momentum conservation to hold. This 
last term is separately conserved in conventional electro- 
dynamics and so can be removed there, but a term of 
this type must be present for T"^ to be conserved in the 
presence of general Lorentz violation. 

C. Covariant dispersion relation 

Much of the propagation behavior of the photon is en- 
coded in its dispersion relation, which provides spectral 
information for the modes. While standard methods can 
be used to find the dispersion relation from the equations 
of motion, at least at leading order in the coeSicients for 
Lorentz violation, handling the gauge freedom typically 
entails the loss of observer Lorentz invariance. In this 
subsection, we present a technique for deriving the exact 
covariant dispersion relation, based on the rank-nullity 
properties of the equations of motion. As a concomitant, 
the technique provides some insight into the nature of 
birefringence, which arises whenever the dispersion rela- 
tion has non-degenerate physical solutions. 

We begin by adopting the ansatz 

A,{x)=A,{p)e-^-P. (21) 

This implies the equations of motion (11) can be ex- 
pressed in the matrix form 

Mi^^A^ = (22) 

with 

= 2r'""'PaP0 + 2iX''->„, (23) 

where it is understood that each occurrence of da in the 
operators kAF and kp is replaced with —ipa- The matrix 



M''" is hermitian, Af'"' = (M'^^)*. It satisfies the con- 
ditions for charge conservation, M'^'^p^ = 0, and gauge 
symmetry, M'^'^p^, = 0, which imply that the determi- 
nant of M^^" vanishes identically. 

The standard method to handle the gauge freedom 
and the vanishing determinant involves making a defi- 
nite gauge choice, thereby reducing the four-dimensional 
problem to a three-dimensional one. The determinant 
of the reduced linear equation yields the dispersion rela- 
tion. Within the particular gauge choice and for a given 
solution to the dispersion relation, one can then solve for 
the polarization mode Af^ of the photon. The general 
solution for this mode is the sum of this solution and an 
arbitrary pure gauge term cx p^. Typically, the gauge 
fixing explicitly breaks observer Lorentz invariance. The 
alternative method presented below focuses on the rank- 
nullity of the linear equation (22), which allows us to 
preserve Lorentz covariance throughout the calculation. 
We take advantage of the exterior product and its ability 
to determine linear independence of a set of vectors. The 
reader is reminded that a set of vectors {A°-}, a= 1,2,... 
is linearly independent if and only if their exterior prod- 
uct is nonzero [51]. 

Starting with a set of arbitrary basis vectors 
{A\A'^,A^,A'^} with A^ A A"^ A A^ A A^ ^ 0, the im- 
age space of M is spanned by the vectors = MA"". 
The dimensionality of this space is equal to the rank of 
M, and it determines the dimension of the solution space 
of Eq. (22). In particular, gauge freedom ensures that M 
is rank three or lower. This implies AB'^ AB^ AB'^ = 0, 
which is equivalent to the condition det M^'^ = 0. So to 
find nontrivial solutions, we must impose rank two or less, 
which implies A S*" A = for any a,b,c= 1, 2, 3, 4. 
This is the minimum requirement, and it leads to the 
covariant dispersion relation. 

We can translate the above discussion into conditions 
on M by noting that M generates for each n a linc^ar 
transformation between n-vectors in its domain and im- 
age spaces. The transformation for a given n is denoted 
A^M and is specified by 

A • • • A B"" = (AfA"i) A • • • A (MA"") 

= (A"M) (A"! A • • • A A"" ) . (24) 

For each n, the linear A"M transformation takes an 
arbitrary n-vector u of the domain space into an n- 
vector (A"M)(a;) of the image space. In a coordinate 
basis, we get an explicit expression for the cofactor ten- 
sors (A"M)^j^2...^^'^^''^'"'^'' controlhng the transforma- 
tion A"M: 

(A^M)^,^,...^/---" = 1m[^/^M,/= . . . M^^]''". 

(25) 

The key point for our purposes is that the rank of M 
is completely determined by the set of A"M transforma- 
tions. For any given n, the matrix M is rank n if and 
only if A"M ^ and A"+^M = 0. Note that A"M = 
implies A"+iM = 0. 
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We can now determine the size of the null space of 
M and hence the nature of the solution space for the 
equations of motion by using the rank-nullity relation. 
Since M is a 4 x 4 matrix, the rank can in principle 
range from four down to zero. The five possibilities can 
be summarized as follows. The case of rank four has 
A^M ^ and is excluded by gauge invariance, so no 
solutions exist. For rank three, which has A^M = and 
A-Hf ^ 0, the nullity is one and hence there is a one- 
dimensional solution space. In the present context, it 
corresponds to pure-gauge solutions. For rank two we 
have A'^M = and A^M 7^ with a two-dimensional 
null space. This case yields a one-dimensional non-gauge 
solution space. For rank one A^M — 0, A^M ^ 0, and 
the nullity is three, so there is a two-dimensional non- 
gauge solution space. Finally, for the case of rank we 
have A-^M = M = 0, which is trivial. 

The above discussion reveals that the covariant disper- 
sion relation ensuring at least one physical solution is 

A^M = 0, (26) 

in which case M is rank 2 or less. In a coordinate basis, 
this dispersion relation takes the tensor form 

iM[^/W^/W^3f 3 = 0. (27) 

Also, if it so happens that requiring A^M = also leads 

to A^M = 0, then a two-dimensional physical solution 
space exists. We interpret this situation as follows. Sup- 
pose wc fix the 3- momenta pj, and find a frequency po 
that solves the dispersion relation (26). If at this fre- 
quency A^M ^ and hence M has rank two, then there 
is exactly one non-gauge polarization mode associated 
with this solution. Other frequencies po that solve the 
dispersion relation (26) lead to different polarizations and 
different phase velocities. The solution is therefore 



birefringent. However, if for a given frequency solution 
we find A^M = and hence M of rank one, then there 
are two independent polarizations that propagate with 
the same phase velocity. This situation represents non- 
birefringence. We therefore obtain the correspondence 

A^M = -(-^ no birefringence, (28) 

which provides the explicit condition for the existence of 
nonbirefringent modes. The next subsection explores the 
issue of birefringence in more detail. 

The covariant tensor dispersion relation (26) can be 

rewritten as a covariant scalar dispersion relation by the 
judicious use of gauge symmetry. To see this, adopt 
the special domain basis {A^ , , ,p}. Since gauge 
invariance implies Mp = 0, the only 3- vectors oj that 
yield nonzero values of (A^M)(a;) must be proportional 
to A^ /\ A^ f\ A'^ . Consequently, the transformation 
{/V'M){uj) is one dimensional, and therefore its dual 

/\^M must also be one dimensional. The hermiticity 
of M implies that there is a vector in terms of 
wjhich the dual transformation /\^M is determined as 
(f^Mdy = Direct calculation with this result 

shows that (A^M^p^ = (A^MYpP^", and from these re- 
lations we find pPpp(f^lY" = [AH^iy pp^'p" . Finally, 
this expression implies that the tensor dispersion rela- 
tion (26) is satisfied for any nonzero p if and only if the 
trace of its dual vanishes. So we arrive at the covariant 
dispersion relation 

(A?MYp = 0, (29) 

which is a scalar density. 

In terms of the constitutive tensors (14), the covariant 
scalar dispersion relation (29) can be written as 



= - lr'"'^{{A^M,)^^,p - 3(a2M„)^„,^) 



1. 



— 3^(Ul/:i2(Ii3jli4"^l'll'2!^3'^4 



e.,„Mp.PpsPp.r'''''''''r'''''''''x'''''''''''' + 8PaP0ikAFUkAF).r''^-. (30) 



In this expression, the dual of A^M is defined by 

(A2M)^a./3 = \e^apje^0a5M>"'M^\ (31) 

while Me and Mq are, respectively, the CPT-even and 
CPT-odd parts of M. 

The covariant scalar dispersion relation (30) is a neces- 
sary condition for the existence of nontrivial plane-wave 
solutions. The first term on the right-hand side of Eq. 
(30) is CPT even and matches the result found in Ref. 
[52] in the appropriate limit. The last term contains all 
CPT-violating contributions. Note that this covariant 



scalar dispersion relation is independent of the spacetime 
metric ry^j/. Note also that the momentum dependence of 
the constitutive tensor j^'^P'^ and of (fcyiF)/^ implies that 
the dispersion relation (30) is typically a polynomial of 
degree greater than four in the frequency po. As a result, 
more eigenfrequencies typically exist in the presence of 
Lorentz violation than for the corresponding situation 
in conventional electrodynamics, and so more modes can 
propagate. However, following a reasoning similar to that 
leading to the equations of motion (11), we expect the 
solutions of interest to be small perturbations of the lim- 
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iting physical solutions in conventional electrodynamics. 
It follows that only the corresponding subset of the solu- 
tions to the dispersion relation (30) are relevant to low- 
energy physics, while the others represent high-frequency 
modes that may play a role as Planck-scale energies are 
approached. 



D. Birefringence 

Astrophysical searches for vacuum birefringence in 
photon propagation provide sensitivities to Lorentz vi- 
olation that are many orders of magnitude beyond those 
attainable via other techniques. It is therefore valuable 
to classify coefficients for Lorentz violation according to 
their birefringence effects. The analysis in the previous 
subsection has already provided some insight via the con- 
dition (28) for the absence of birefringence. In this sub- 
section, we provide a decomposition of coefficients for 
Lorentz violation that distinguishes birefringent and non- 
birefringent cases. We also offer some remarks about 
generic conditions for birefringence and their connection 
to the metric and to electromagnetic duality. 



1. Coefficients for birefringence 

In the minimal SME, leading-order birefringence is 
known to be controlled by the (i = 3 coefficients {kAF)ti 
and by a subset of the d = 4 coefficients (Aif)'*^'^" [7, 11]. 
In the present context with Lorentz-violating operators 
of any dimension, a similar pattern holds: leading-order 
birefringence is associated with all the odd-d coefficients 
itAF and with some combinations of the even-d coeffi- 
cients kp. One way to verify this is via the covariant 
scalar dispersion relation (30). 

To identify the relevant even-d coefficients, it is useful 
to introduce the definitions 





= ^{kde + khb) - ^Tr{kDE + khb) 


l^e— 


= ^{kde - khb) - IT^{kde - khb) 


ko+ 


= ^{kdb + ^<he)- 


Kq — 


= |(kdb - khe), 


Ktr+ 


= ^Tr{kDE + khb), 


Kt, — 


= ^Ti{kDE - Khb), 



which give an experimentally judicious decomposition of 
the coefficients for Lorentz violation appearing in Eq. 
(18). The ffist four of these are traceless 3x3 matri- 
ces, with ke+, kg-, ko- symmetric and ko+ antisymmet- 
ric. The last two are S0(3) rotation scalar combinations. 
Note that ktr+ can be disregarded in the minimal SME 
because it represents a simple Lorentz-invariant scaling 
factor in that context, but nonrenormalizable terms of 
this type can violate Lorentz symmetry and so must be 
included in the present context. 



Among the combinations (32), only the matrices Kg-i- 
and Ko_ cause vacuum birefringence at leading order. 
This result can be understood in terms of a Weyl decom- 
position of the constitutive tensor x^'^^'^ analogous to the 
Weyl decomposition of the Riemann tensor, 

+ v'^'icFT" - v'^^icFV") + C^"""". (33) 

In this equation, the tensor C'^'^p'^ corresponds to the 
Weyl component and is traceless, C^^'^P'^rj^a = 0. The 
term corresponding to the Ricci component involves 
(cf)"^, which is defined as the symmetric combination 

(c^)«/5 ^ {kFr''^ - \{kFr\x^- (34) 

The relations between (cf)"^, C'^"'"^ and the k matrices 
are 

(cf)°" = \{rKtr-+ktr+), 

(cfY^ = -{ke-y^ + ^{{ktr-} - («tr+))<5^', 

{cfT' = -l{ko+re^''', (35) 

com ^ 

/yjklm ^ { C- \np jkn I mp 

C-OJ" = i(K„_)^™e^'". 

We see that the ten independent components of cf are 
equivalent to the ten independent components of kg-, 
ko+, ktr+, and kt,--, while the ten independent compo- 
nents of C match those of ke+ and ko- ■ 

The decomposition (33) reveals that the coefficients 
cf play the role of a small distortion of the spacetime 
metric at leading order. In this respect, they are anal- 
ogous to the c"^ coefficients in the matter sector of the 
SME, which motivates the notation. As further discussed 
below, a small metric distortion leaves unaffected the 
usual degeneracy between polarizations and so cannot 
cause birefringence at leading order. In contrast, the 
non-metric Weyl piece of Eq. (33) breaks the degener- 
acy and causes birefringence. Note that the effects of 
c-type coefficients in the minimal SME are unobservable 
in experiments involving only one sector, since they can 
be removed by a judicious coordinate choice. However, 
in the present context where dp depends on energy and 
momentum, dispersion effects may arise that are observ- 
able. 

The above decomposition and results hold for the vac- 
uum, where the Weyl decomposition is performed using 
the Minkowski metric. However, the effective metric g^^, 
for electrodynamics in a macroscopic medium M is no 
longer Minkowski, so the above decomposition must be 
modified. The Weyl part C'^^ is required to be trace- 
less with respect to ^f^^ instead, so the assignment of k 
matrices to and differs. The attribution of bire- 
fringence effects to a coefficient can therefore be medium 
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dependent. It also follows that an experiment in a suit- 
able medium can achieve sensitivities to different coef- 
ficients compared to the same experiment performed in 
vacuo. 

As an example, consider a uniform isotropic medium 
M with refractive index n and permeability fi. Perform- 
ing a decomposition of the constitutive tensor reveals 
that the relations between , and the k matrices 
become 



+ \^{{ktr-) - {ktr+))S^'", 

[cffi = -l^[k,^f'e^^\ (36) 

('t?^ _|_ n 

//^M\jklm \"' ' ' f \np jkn Imp 

/' \np jkn Imp 

^^M^ajkl _ ij^^^ ym^klm 

For this type of medium, we see that birefringence is 
associated with the matrix well as the matrices 

Ke+ and kg-. Moreover, the matrix ke+ can now affect 
nonbirefringent phenomena at leading order, unlike the 
vacuum case. Additional coefficient mixings can be ex- 
pected in anisotropic and gyroscopic media. 



2. Birefringence, metric, and duality 

An interesting challenge is the identification of the min- 
imal set of requirements leading to birefringence. For 
simple local constitutive relations, nonbirefringence is 
known to be associated with a pure-metric constitutive 
tensor [53, 54], while electromagnetic duality plays a role 
[55]. The idea that the essential properties of electrody- 
namics rely on constitutive relations rather than the un- 
derlying metric structure of the spacctimc is a key aspect 
of the premetric approach to electrodynamics [52, 56]. 
In this subsection, we ofi^er some remarks on the role of 
the metric and of electromagnetic duality in determining 
birefringence conditions within the context of the general 
constitutive relations (14). 

Our primary conjecture is that the only nonbirefrin- 
gent terms arise from the non-Weyl component of ■xj^'^P" . 
Wc therefore sock a procedure for extracting this com- 
ponent. The structure of the Weyl decomposition (33) 
suggests the possibility of introducing an effective met- 
ric cjfj^i, in terms of which the Ricci component vanishes. 
This would imply a natural decomposition of into 
metric and non-metric components or, equivalently, into 
nonbirefringent and birefringent components. 

We therefore postulate the existence of an effective 



metric satisfying the eigenproblem 



(37) 



where the scale factor for gf^^, is chosen so that the propor- 
tionality constant in this equation matches the value for 
standard electrodynamics. Note that in the present con- 
text the effective metric g^p typically depends on the 4- 
momentum and therefore cannot be interpreted as a con- 
ventional spacetime metric. The existence and unique- 
ness of solutions to Eq. (37) is an interesting open issue, 
but here we suppose all physically reasonable constitutive 
operators x^^''P'^ lead to a solution that is unique up to a 
sign. For practical purposes, this issue is moot because 
it suffices to use Eqs. (14) and (37) to find a perturbative 
expansion for g^^. 

The effective metric provides a unique Wcyl de- 
composition of the constitutive tensor ^'^'-'P'^ , given by 



X 



Xn 



(38) 



where Xw is the trace- free Weyl component, 
iXw)'^'''"^ gi^a — 0. Our conjecture now identifies 
Xw as the birefringent component. We can verify this in 
special limits. For vanishing Xw and no OPT violation, 
the dispersion relation (30) reduces to 



(39) 



as expected. Also, for this case we find that that the co- 



factor tensor A M is proportional to {g 



which 



demonstrates the absence of birefringence. Another limit 
of interest involves small Lorentz violation. Solving at 
leading order in kp, we find 



ixv.r^^'' 



(40) 



This is also consistent, since it demonstrates that C 
corresponds to the leading-order birefringent component 
and c_F to the nonbirefringent part. 

In conventional electrodynamics, electromagnetic du- 
ality ensures no birefringence occurs in the vacuum. 
Given the 2-form field strength F and the dual 2-form 
H = *F, the Maxwell equations in vacuo take the com- 
pact form dH = 0, dF = 0. Here the *-operator is 
defined as = \^iiv'"^ and obeys ** = — /. In 

this language, electromagnetic duality can be under- 
stood as the statement that if F is a solution then so is 
F' = exp(a*)F = Fcosa + H sma. This can be viewed 
as a rotation between F and H. The chiral components 
F± = i(l±i*)Fof the field strength are irreducible 
representations of these duality rotations and also of the 
Lorentz group. The relevance to our discussion is that 
duality symmetry excludes birefringence because it im- 
plies the space of plane- wave solutions is two dimensional. 
There are two independent polarizations that mix under 
duality transformations, and the duality symmetry en- 
sures that both polarizations propagate with the same 
phase velocity, so no evolution of polarization can occur. 



10 



This conventional duality can be generalized to the ef- 
fective metric Qfj^i, relevant in our context. For an arbi- 
trary effective metric, we define 



* = sJ-deig*{A^g-^), (41) 

which is normalized to satisfy H = —I as usual. This 
operator induces a chiral structure with respect to 
The projections of the field strength F onto the chiral 
subspaces are now 

F± = i(l±ii)F. (42) 
Neglecting CPT-odd terms, the dual field strength 

H = *G = *xF (43) 
can also be decomposed into chiral components, 

H± = *xF± = U'^±ii)H. (44) 

According to our primary conjecture, the constitutive 
tensor x provides a natural decomposition of the solu- 
tion space. It is therefore plausible that the absence of 
birefringence is associated with a duality symmetry gen- 
erated by *, with any nonzero Weyl piece Xw or nonzero 
CPT violation breaking this symmetry and so causing 
birefringence. We can explore this idea for the case of 
constitutive relations of the metric type, x = ^^5^^- 
This form of the constitutive tensor yields the closure 
relation *x* x = 1/ det^. Under these circumstances, 
if F represents a solution to the Maxwell equations 
dF = dH = 0, then it can be shown that F *xF 
is also a solution. This generalizes the result for con- 
ventional duality rotations to the *-chiral subspaces. It 
follows that every polarization is associated with a sec- 
ond polarization that is also a solution and hence that 
birefringence is absent, as expected. 

A rigorous derivation of the above results is an open 
problem of definite interest, although it lies beyond our 
present scope. Note that our arguments hold inside 
any simply connected source-free region. We anticipate 
that they can also bo applied to more general scenar- 
ios with inhomogeneous constituent tensors and curved 
spacetimes, in which duality would be defined locally. 
It is also plausible that duality breaking could arise via 
boundary conditions. This could lead, for example, to 
a foundational understanding of degeneracy splitting in 
resonant cavities. Generalizations of the framework may 
also merit investigation. For instance, more complicated 
closure relations of the form *x* x = a + b * x for scalar 
a and b also result in duality symmetries, and they may 
be related to improved decompositions of the constitu- 
tive relations. One could also consider a duality at the 
potential level, involving mixing of two 1-forms A and B 
obeying F = dA and H = dB, which can be found for 
any vacuum solution. Such approaches may make it pos- 
sible to incorporate CPT-violating effects in this picture. 



E. Effects of Lorentz violation 

It is of interest to categorize the types of effects pro- 
duced by Lorentz violation in various physical situations. 
A number of schemes are possible. In this subsection, 
we offer some remarks about the categorization used 
in later sections, which is based on identifying Lorentz- 
violating effects in terms of birefringence, dispersion, and 
anisotropy. 

In the previous subsections, we have defined birefrin- 
gence as the existence of only one low-energy eigenmode 
for a particular solution to the covariant dispersion re- 
lation (30), and we have conjectured its connection to 
a breakdown of duality. This definition contains more 
than the notion of rotation of the polarization of light in 
the vacuum. The breaking of eigenmode degeneracy can 
cause effects in circumstances other than vacuum prop- 
agation, such as the splitting of resonant frequencies in 
cavities. In the remainder of this work, the term birefrin- 
gence is used in the sense of eigenmode nondegeneracy. 
In vacuum propagation, this reduces to the usual notion 
of rotation of polarization and can be termed vacuum 
birefringence. 

A similar dichotomy appears in the definition of dis- 
persion. In what follows, we adopt the term dispersive to 
refer to Lorentz- violating operators in the Lagrange den- 
sity that appear with other than two derivatives. Only 
operators with mass dimension = 4 are nondispersive 
in this sense. In the momentum-space covariant disper- 
sion relation, these operators contribute terms that are 
non-quadratic in the momentum. However, this defini- 
tion implies more than merely a nonlinear relationship 
between the frequency po and the momentum p for vac- 
uum propagation. As shown explicitly in Sec. IV E, some 
operators involving non-quadratic derivative terms and 
hence labeled as dispersive according to our usage in fact 
produce no leading-order dispersion in vacuum propa- 
gation, although they can produce analogous effects in 
other situations such as cavity resonators. Throughout 
this work, we use the term dispersive to mean d 7^ 4, and 
we reserve the term vacuum dispersion for modifications 
to the usual relation po = |p| for vacuum propagation. 

In any physical situation, the properties of electro- 
magnetic waves are determined by the coefficients for 
Lorentz violation, the medium, and the boundary con- 
ditions. Some useful intuition can be gained by consid- 
ering the interplay between these and the resulting bire- 
fringent, dispersive, and anisotropic effects. 

Consider first the presence of a macroscopic medium. 
We have seen in the previous subsection that a coeffi- 
cient controls birefringence if it is associated with CPT 
violation or if it contributes to the Weyl piece of the con- 
stitutive tensor x^'^'"^ ■ The presence of a medium can 
affect the Weyl decomposition, as shown in Eq. (36), so 
the manifestation of Lorentz-violating birefringence de- 
pends on the medium. In contrast, a coefficient controls 
Lorentz-violating dispersion according to the derivative 
structure of the corresponding operator in the Lagrange 
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density, which is unaffected by the medium. Similarly, 
Lorentz-violating anisotropic effects are determined by 
the rotation properties of the relevant operators, and so 
they too are independent of the medium. 

The role of the boundary conditions is different. A 

given choice of boundary conditions determines which 
coefficients for Lorentz violation are measurable. This 
feature is similar to a property of conventional electrody- 
namics. A given solution, such as a vacuum plane wave or 
a cavity eigenmode, is determined by both the Maxwell 
equations and by the boundary conditions. Changing 
the boundary conditions reveals distinct sets of cigcnsolu- 
tions with different physical properties. One cannot, for 
instance;, expand vacuum plane waves in terms of cavity 
eigenmodes. Similarly, in the presence of Lorentz viola- 
tion, a specific choice of boundary conditions fixes certain 
eigenmodes as solutions. However, only a suitable subset 
of the coefficients for Lorentz violation affects a given set 
of eigenmodes. For example, far-field boundary condi- 
tions suitable for vacuum radiation yield solutions that 
depend on a much reduced subset of coefficients, as is 
explicitly identified in Sec. IV C. 

Despite their role in determining the observability of 
coefficients for Lorentz violation, the boundary condi- 
tions have no impact on the associated physical effects. 
For example, the birefringence properties of a given oper- 
ator are unaffected, essentially because the duality sym- 
metry is a local property and therefore is independent 
of the boundary conditions. Dispersive properties are 
also unaffected because they are associated with position- 
space derivatives on the Lorentz-violating operators and 
therefore are basically a local feature in position space. 
The associated modification of the momentum-space re- 
lation between po and p reflects an impact on the eigen- 
solution space as described above, rather than a change 
in the underlying dispersive properties of the operators. 
Similarly, the anisotropy properties of Lorentz-violating 
operators are fixed by the Lagrange density and are in- 
dependent of the boundary conditions. 

This feature of the boundary conditions has the in- 
teresting implication that a comprehensive search for 
Lorentz violation requires multiple observational and ex- 
perimental methods, since any one method typically ap- 
plies one type of boundary condition and so cannot be 
expected to access the whole coefficient space. It also 
means there are multiple approaches for categorizing the 
coefficients. For example, one can split the coefficient 
space into a subset selected by boundary conditions ap- 
propriate for vacuum propagation and its complement. 
This turns out to be an apposite splitting for several rea- 
sons, and we develop it in some detail beginning in Sec. 
IV. However, one could in principle consider alternative 
splittings of the coefficient space using other boundary 
conditions, such as ones for resonant cavities. 



III. SPHERICAL DECOMPOSITION 

In the previous section, we have demonstrated that 
the effects on photon propagation of Lorentz-violating 
operators of arbitrary dimension are specified by the 
Lagrange density (8) and are determined by the co- 
efficients for Lorentz violation (fc^^)^ and 
^^(d))KA^x^ai...a(i_4), A classification of these coefficients 
that characterizes the key effects of the corresponding 
Lorentz-violating operators is both useful and convenient 
for more detailed investigations. The ideal scenario is to 
establish a minimal collection of independent coefficients 
associated with operators having physical properties of 
direct relevance to observation and experiment. 

One natural classification scheme takes advantage of 
the role of spatial rotations to perform an SO (3) decom- 
position of the coefficients for Lorentz violation. This 
technique has been applied to obtain first measurements 
of certain coefficients relevant to vacuum photon prop- 
agation from astrophysical observations of active galax- 
ies, gamma-ray bursts, and the cosmic microwave back- 
groimd [22, 23]. The SO (3) decomposition uses spin- 
weighted spherical harmonics [57, 58] , which are angular- 
momentum eigenstates and so obey relatively simple 
transformation rules under rotations. The method has 
the advantage that spin-weighted spherical harmonics are 
commonly used in some areas of astrophysics [59] and are 
well understood. A summary of some properties of spin- 
weighted spherical harmonics is provided in Appendix A, 
which also derives several mathematical relations used in 
what follows. 

In this section, we discuss the decomposition of the 
coefficients Icaf and kp into spin- weighted components. 
For the analysis, the coefficient kAp is separated into the 
pseudoscalar (fcAF)o and the pseudovector Uaf, while 
kp is separated into the tensors kde and khb and the 
pseudotensor kdb defined in Eq. (18). Each of these 
five components is expanded in a helicity basis and de- 
composed into spin-weighted spherical harmonics. The 
symmetries of kAp and kp then permit extraction of a 
minimal collection of spherical coefficients for Lorentz vi- 
olation. 

The results in this section demonstrate that the mini- 
mal collection of spherical coefficients includes nine sets 
of coefficients. For convenience, we summarize the no- 
tation here. Three sets are extracted from kAP and are 

denoted {k^Al)ffrl^ i^^f^j^^ and (^AF)nf2- Six sets 

emerge from A:f, denoted {k^p )nfl, i^F )nj2^ 

(^F )njrn' (^F )njm > and ( fcp. Y^jJ. The Symbol c Spec- 
ifies coefficients associated with nonbirefringent opera- 
tors, while k specifics birefringcnt ones. A negation di- 
acritic -> denotes coefficients that have no leading-order 
effects on the vacuum propagation of light, a property 
derived in Sec. IV E. The subscripts n, j, and m deter- 
mine the frequency or wavelength dependence, the total 
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angular momentum, and the z-componcnt of the angu- 
lar momentum, respectively. The allowed ranges of these 
three indices and the counting of the coefficients are given 
in Tables I, III, IV, and V. The label d gives the mass 
dimension of the corresponding operator for Lorentz vi- 
olation, while the numerals 0,1, or 2 preceding either E 
or B refer to the spin weight of the operator. The su- 
perscripts E and B refer to the parity of the operator, 
with parity 1)-' labeled as E and parity (— 1)-'+^ as B. 
The phases are chosen so that each spherical coefficient 
)Cjm for Lorentz violation obeys the complex- conjugation 
relation 



ilCjmT = (-1)™/C 



j(-m)- 



(45) 



Note that this implies that the pair of complex coeffi- 
cients JCjm and lCj(^m) are codependent but represent 
two real degrees of freedom. 



A. General CPT-odd coeflScients 

We begin the decomposition of (fc^i?)^ by performing 
an expansion of the pseudoscalar component {kAF)o in 
spherical harmonics. The first step is to separate (fcAF)o 
into pieces with definite frequency and 3-momentum de- 
pendence. Denoting the frequency by w = p° and the 
components of the 3-momentum p hy p'^, we obtain the 
expression 

d n=0 



0...0/ci...fe„ 



-3— n„fei 



p'^'.-.p"". (46) 



Writing the magnitude of the 3-momentum as p = \p\, 
we see that each term in the sum involves a factor of 
^d-s-n^n 'pj^g frequency and wavelength dependence 
of each term is therefore controlled by the new index n. 
The direction dependence introduced by the components 
of p is characterized through the expansion in spherical 
harmonics below. 

To determine the relevant angular-momentum eigen- 
values j for the spherical-harmonic expansion, we can 
break the coefficients for each n into a series of three- 
dimensional traceless symmetric tensors of rank n, n — 2, 

n — 4, This implies the range of eigenvalues I = n, 

n — 2, n — 4, . . . for the orbital angular momentum. Since 
the spin of (fc^F)o is zero, the eigenvalue j must also span 
this range. An alternative and more elegant approach 
makes use of parity. For a given n, the orbital angular 
momentum is limited by n. Since the spin is zero, we 
must have j < n. Given that (^af)" is a pseudoscalar 
and p is a vector, each term in the expansion must have 
parity (—1)"+^ with only S-type parity occurring. This 
imposes the relation (—1)"+^ = (—1)-'+^, which implies 
that j — n is even. The conclusion is that j = n, n — 2, 
n — 4, . . . > 0, as before. 



by 



The resulting expansion in spherical harmonics is given 



(W)o = E '^'"'">" oYjmip) {kfXf2- (47) 

dnjm 



The spherical coefficients for Lorentz violation {k^AF)n^^ 



/njm 

are nonzero for the n and j values listed in Table I. As an 

example, consider d = 5. In this case, there are two j = 

singlets, {k'"^p)lx)o ^ ^^'-^ (^af)20^^- There is also one j = 
1 triplet (fc^^f')iim , with m = —1, 0, 1. Finally, there is 
one j = 2 quintuplet {k^^p)^2m ^ with m = —2, —1, 0, 1, 
2. The total number of coefficients is 1 + 1 -|- 3 -|- 5 = 10. 

Next, we perform a spherical-harmonic expansion of 
the radial component {kAF)r — P ■ ^af — Sr ■ ^af- This 
component is also a pseudoscalar, and the decomposition 
follows the same basic steps as above. However, the addi- 
tional p factor implies that the total angular momentum 
is now limited by j < n-l- 1. Also, the parity of each term 
is given by (—1)" = (—1)-'"'"^, yielding the range of j as 
j = n + 1, n — 1, n — 3, . . . > 0. We therefore obtain the 
expansion 



{kAF)r 



dnjm 



-'-VoYjm{p){kfJ°fJ, (48) 



involving another set of spherical coefficients ik'"^p)^^\ 
However, it turns out that the symmetries of the tensors 

(^af)k ^' imply that these new coefficients can 

all be expressed as combinations of the other spherical 
coefficients occurring in the expansion of {kAF)K- Before 
demonstrating this, we first complete the expansion of 

kAF- 

The remaining components of kAF have spin weight 
±1, {kAF)± = e± • kAF- For these cases, we again find 
j < n + 1, but now coefficients with both E- and B- 
type parities occur. The parity is (—1)", which implies 
j = n + 1, n — 1, ...> 1 for the _B-type components and 
j = n, n — 2, . . . > 1 for the E'-type components. Since 
the spin weight is ±1, the index j for the total angular 
momentum is limited from below by 1 . The expansions of 
the components {kAF)± in terms of spin- weighted spher- 
ical harmonics take the forms 



ikAF)± = J2 ^^'""'V ±lYjmip) 



dnjm 



f=^(±(4%"J+^(^?i)S), (49) 



where we introduce a factor of 1/ \j2j(j -|- 1) for later 
convenience. The spherical coefficients for Lorentz vio- 
lation ik^AF)nj^ s-'^d (k:AF)njwl represent the remaining 
two independent sets. Their ranges are summarized in 
Table L 

The symmetries of {k^AF)^ ' ' described in Sec. 
II B imply certain constraints on the four sets of coeffi- 



Cienib yl^AFJnjm' y^AFJnjm ' \'^AF)njm^ '^nu [i^AFJnjm- 



'(d)^(OS') /,(d)^(li3) 



(IE) 



13 





/ , (rf) UOB) 






n 


3 


3 


3 



1 
2 
3 
4 

d-3 




1 

2 

1 3 
2 4 

2 ••• d-3 


1 

2 

1 3 

2 4 
13 5 

1 3 ••• d-2 


1 

2 

1 3 
2 4 

2 ••• d-3 


total 


id(d-l)(d-2) 


i(d - l)(d2 +d-3) 


i(d+l)(d-l)(d-3) 



TABLE I: Summary of the allowed ranges of indices n and j for the independent spherical coefficients associated with CPT-odd 
operators. The dimension d is odd with d > 3, while n < d — 3. The index m satisfies the usual restrictions —j<m< j, so 
there are 2j + 1 coefficients for each j. For a given dimension d, the number of coefficients of each type is given in the last row. 
Adding these gives the expected total of |(d + l)(d — l)(d — 2). 



To determine these constraints, the symmetries must be 
formulated in the momentum-space helicity basis. Recall 



that the tensors {k^p)^ 



(d) X "l---«(d-3) 



are totally symmetric in 



the last d—3 indices and that any trace involving the first 
index vanishes. The total-symmetry condition is implicit 
in the above spherical-harmonic decompositions, but the 
trace condition provides a nontrivial constraint. 

In momentum space, the trace condition can be written 
as the differential equation 

d - d - 

= ^ {kAF)K = -^(fcAF)o + V • kAF 

OPk OUJ 

d . 

= -g^{kAF)o + '^r{kAF)r + V +{kAF)- + V_(A;af) + 
= ^(fcAF)0 + (^ + -)(W). 



1 



+ -{j+{kAF)- - J-{kAF) + ), 
p 



(50) 



where we have used the identities (A37). The combi- 
nation J^{kAF)- — J-{kAF)+ is a pseudoscalar that is 
generated by the B component of kAF and that con- 
tains no E component. Consequently, Eq. (50) provides 
a symmetry constraint involving only the i?-type coef- 
ficients {k^AlinSi' i'^Al)nfrJ' ^nd (A:^^].)^^^ Insert- 
ing the spherical-harmonic expansions into Eq. (50) and 
making use of the identity (A35) yields an explicit rela- 
tion involving these three coefficient sets. Careful consid- 
eration of the index ranges then reveals that (A;^^)^°^^ 
may be written as 



(fc 



, I o\\'^AF)r. 



AF /njm 



/njm 



+ (d-2-n){k^fj;i%.J. (51) 



We conclude that the auxiliary coefficients {k^AF)nfm 
completely determined as linear combinations of the in- 
dependent coefficients {k^fH^f^ and {k^Al)nfrl- 



The net yield of the helicity decomposition of (kAF)K is 
therefore three independent sets of spherical coefficients 

for Lorentz violation, which are the i?-typc coefficients 
(^AF)n7m ^nd ik'^^l)^^^], and the -E-type coefficients 



Ad) 



/njr. 
(IE) 



( ^AF Jnjm ■ Thesc Coefficients completely characterize the 
CPT-odd Lorentz-violating operators of arbitrary mass 
dimension associated with the quadratic action for elec- 
trodynamics. All the operators are birefringent. On each 
set of coefficients, the label d specifies the (odd) mass di- 
mension of the operator. The three indices n, j, and m 
determine the frequency or wavelength dependence, the 
total angular momentum, and the z-component of the 
angular momentum, respectively. The allowed ranges of 
these three indices and the counting of the coefficients 
are given in Table I. The superscripts E and B refer to 
the operator parity, while the superscript numerals or 
1 preceding either E or B specify the spin weight of the 
operator. The total number of coefficients for given odd 
d is i(d-t- l)(d— l)((i— 2), matching the group-theoretic 
result from Sec. II. 



B. General CPT-even coefficients 

The spherical decomposition of the CPT-even Lorentz- 
violating operators (kp)'^^^" follows a procedure similar 
to that for the CPT-odd case discussed in the previous 
subsection. However, instead of working directly with 
(kpY^^'^ , it is more convenient to decompose the three 
matrices kde, khb, and kdb given in Eq. (18), which 
have equivalent content. 

We begin by separating the three matrices kde, i^hb, 
Kdb into their S0(3)-irreducible trace, symmetric trace- 
less, and antisymmetric parts. These irreducible parts 
can be expressed in the helicity basis via Kat = • k • e^, 
where a,b = +,r,—. Each component can then be ex- 
panded in the appropriate spin-weighted spherical har- 
monics, recalling that koe, khb are tensors while kdb 
is a pseudotensor. At this stage, we find the results can 
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. (d) aOE) /(d) ^(0i5) 
K'^DElnjm i K'^HBJnjm 5 

/^(d) ^{OS) 


(Jd) UlB) |.^(d) .(IE) 
y^DEInjrm y^HBInjrm 

(Ad) n(1B) 


/(d) .(2B) ,(d) .(2E) 
y^DEJnjrm X'^HBInjrm 

(Ad) ^(2B) 

K'^DBJnjm 


/^(d) ^(OE') , (d) .(OE') 
K'^DEInjm > K'^HBJnjm 


n 




j 


3 


3 



1 
2 
3 
4 

d-4 


2 

1 3 

2 4 

13 5 
2 4 6 

2 ••• d-2 


2 

1 3 
2 4 

13 5 
2 4 6 

2 4 ••• d-2 


2 

3 

2 4 

3 5 
2 4 6 

2 4 ■•• d-2 




1 

2 

1 3 
2 4 

2 ••• d-4 





(Ad) -.(IS) / (d) x(li3) 

(Ad) AlE) , (d) x(liS') 
y^DBInivm y^DB/nim 


(Ad) A2B) , (d) x(2B) 
(Ad) -,(2B) 


(Ad) -,(0S') 
y^DBJnjm 


(Ad) AlB') 
("'DBJnjm. 


n 


3 


3 


3 


3 



1 
2 
3 
4 

d-4 


1 

2 

1 3 

2 4 
1 3 5 

1 3 ••• d-3 


2 

3 

2 4 
3 5 

3 ••• d-3 


1 

2 

1 3 
2 4 

1 3 5 

1 3 ••• d-3 


1 

2 

1 3 
2 4 

2 ••• d-4 



TABLE 11: Summary of the allowed ranges of indices n and j for the 20 sets of codependent spherical coefficients associated 
with general CPT-even operators. The dimension d is even with d > 4, while n < d — 4. The index m satisfies the usual 
restrictions —j<m<j, so there are 2j + 1 coefficients for each j. 



be expressed as 13 equations involving sums over twelve 
sets of E'-type and eight sets of B-type coefficients, for 
a total of 20 sets of codependent coefficients. The chal- 
lenge is to use these equations and the symmetries of 
(kpY^^" to identify the independent sets of coefficients 
among the 20 codependent ones. The 13 equations are 
given explicitly below. In them, all sums are restricted 
to even dimensions d > 4. The maximum range of the 
indices spans 0<n<d — 4forn and < j < n + 2, 
—j < m, < j for the eigenvalue indices of the angular mo- 
mentum. Details of the allowed index ranges for each of 
the 20 sets of codependent coefficients are given in Table 
II. 

For the matrix kde, six sets of coefficients are needed, 
and the result is 



{iiDE)±r= a;'-^->"±ilS-m(p) 

dnjm 

^ 2j(j+l) ^ y'^DEJnjm ^ ''K'^DEJnjm ) ^ 

(52c) 

dnjm 



Y (j+2)! yy^DEJnjm ^ ''X'^DEJnjni ) ■ 

(52d) 



In this decomposition, the coefficients {K^ijE)njyn emerge 
from the trace component of Kde- 

Another six sets of coefficients are needed for the ma- 
trix khb- The corresponding equations are 



dnjm dnjm 

^ {i''^DE)nj2 + i'^DE)n]m^) ' (52a) ^ (('^ifs)^ J + ('^ffB)ijm'') ' (53a) 

dnjm dnjm 

^ {~ + (^DE)ijm'')' (52b) ^ K^L'b)!^^ + i'^HB)njm^)^ (53b) 
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dnjm 

^ 2j(j+l) ^ y'^HBlnjm + ^K'^HB )njm I ^ 

(53c) 

(/«^,b)±± = ^ c^'^-^-VisF.-^Cp) 

>^ Y (j+2)! (>l«'iJsJnjm ^ ^K'^HBJnjm ) ■ 

(53d) 

The coefficients (k^^)^*^^^ correspond to the trace com- 
ponent of khb- 

Finally, eight sets of coefficients are required in the 
expansion of kdb- The results are 

{iiDB)rr= ^ Oj'^-'-^" oYjm{p) (4'].)iS' (54a) 
dnjm 







( 




(2B) 
"' 




n 




3 



1 
2 
3 
4 


2 

1 

2 

1 

2 


3 

4 

3 5 
4 6 


2 

3 

2 4 
3 


5 






d-4 


2 


4 ••• d-2 


3 




d- 


-3 


total 


i(d^ -4d- 18) 


i(d + 3)(d- 


-2)(d- 


4) 



TABLE III: Summary of the allowed ranges of indices n and 
j for the independent S-type spherical coefficients associated 
with CPT-even birefringent operators. The dimension d is 
even with d > 4, while n < d — 4. The index m satisfies the 
usual restrictions —j<m<j, so there are 2j + 1 coefficients 
for each j. For a given dimension d, the number of coefficients 
of each type is given in the last row. 



dnjm 



dnjm 

X 



1 / + f^id) , .{IE) 

j-jj I ^ y^DBInjm ^ ''V'^DB/njm 



,(Jd) ^i^B') (d) .(lE')N 

^X'^DBJnjm ''\'^DB)njm ji 



(54c) 



{kDB)±r= 1^ a;'^-"->"±lljm(p) 



dnjm 



X 1 ( -L ,-((d) .{IE) 

\/ 2j(j+l) ( + HI^DB)nj 

^(Jd) ^ilB') , (d) 

+ K'^DBInjm ^ ''K'^DBInjm ) 



(54d) 



(«bb)±± = ^ a;''-^->"±2F,-„(p) 



dnjm 



Y JJ+^.\\'^DB>njm ''K'^DB) 



{d) ^,(2i^)^ 



(54e) 



The sets of coefficients (K^i)i^j^\ ('«D]3)i\m^ and 



^DB/njm 



('*DB)i">n' ^'''^ associated with the antisymmetric part 
oi kde, which corresponds to ko+ in Eq. (32). 

The next step is to use the symmetries of the ten- 
sors (fc^))'«^A*i'ai ■■a(d-4) determine the independent 
sets of coefficients. These symmetries are described in 
Sec. II, following Eq. (6). As in the CPT-odd case, the 
independent coefficients can be identified by expressing 
the symmetries as differential equations in momentum 



space. It turns out that all the symmetries of fc^^ are im- 
plicit in the above 13 equations, except for the condition 
of vanishing antisymmetrization on any three indices of 
(j^^d)-^K\nvai...oni_i,) _ In momentum space, this symmetry 
can be imposed by requiring 



= d''P{kFT''^''^ 



(55) 



where d'^ = d/dpn- In terms of k matrices, this require- 
ment becomes the four conditions 

VbikHBT'' = 0, 

^bikoBT^ = 0, 

|^(KDB)"' = e'^''Ve(ftDs)/, 

^(^HbT" = -£''^''Ve(KZ5B)/ (56) 

in the notation of the Appendix. Note the similarity to 
the Maxwell equations. 

The symmetry constraints (56) lead to relations among 
the 20 codependent coefficients in the expansions (52), 
(53), and (54). The procedure for determining these re- 
lations involves expressing the differential equations (56) 
in terms of the operators d/div, 3/ dp, and J±, and then 
inserting the expansions. Following this procedure, we 
find that Eq. (56) provides nine constraints on the eight 
-B-type coefficients and eleven constraints on the twelve 
ii'-type coefficients. However, these constraints are not 
all independent. 

Some calculation with the nine constraints on the B- 
type coefficients shows that there exist only two linearly 
independent sets of i?-type coefficients for Lorcntz viola- 
tion, which we denote by {kp^)nf2. and (^p )nf2.- The 
index ranges for these two sets are summarized in Table 
III. All the corresponding Lorentz-violating operators 
produce birefringent effects. In terms of these two sets. 
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the eight codcpcndcnt B-type coefficients appearing in 
the expansions (52), (53), and (54) are given by 



(-^^l)S = ('^-3-n)n(fc(^));r_l),„, 



(57a) 



= -(^ - 3 - n)(, + 2)ij ^ 
-id-3-n){d-2- n){lf)l%^.^, (57b) 



/ (d) .(IB) _ I ^ o^^'^^'■''^(2-B) 



(4'l3)S = + 2)(n + i){li^Ffnf2 



(57c) 
(57d) 

{>^'Xri = i^+m'^'f:rri, (57e) 



njm ' 



^HBlnjm 

M) \(0B) 



= -^(«(« + 3) + j(j + 



F /njm 



-W-i-n){^f)l%^, (57f) 

= + 2)(n + 3) - j(i + 1))(4'')S™ 

+ i(d-3-n)(^?^)g%^.^, 

= W + 3)(j + 2)(i - l)(4'^)i^^) 



njm 

(57g) 



+ (rf - 3 - n)(n + 2)(:^?^)g!\) (57h) 



The same procedure appUed to the eleven constraints 
on the -E-type coefiicients reveals that there are four lin- 
early independent sets of E'-type coefficients for Lorentz 
violation. In this case, however, some components of the 
corresponding operators are birefringent while others are 
nonbirefringent. For case of application in searches for 
Lorentz violation, it would be useful to have a decom- 
position that separates birefringent and nonbirefringent 
components. At leading order, this separation is achieved 
by the Weyl decomposition (33). Unfortunately, the rep- 
resentations from the Weyl decomposition are incompati- 
ble with the i?-type coefficient representations contained 
in the symmetry constraints (55). Although the latter 
indeed include some operators that are purely birefrin- 
gent and others that are purely nonbirefringent, there 
are typically also birefringent operators with coefficients 
that contribute to the nonbirefringent matrices Ke-, ko+, 
ktr+, and /ttr-- Among the four independent sets of E- 
type coefficients for Lorentz violation, there are thus ones 
that contribute to (cp)'^'^, ones that contribute to (Ji^^t^^ ^ 
and ones that contribute to both. However, we do find 
that only a small fraction of the coefficients contribute 
solely to C'^^i^'^ , Consequently, it is reasonable to seek a 
decomposition separating the nonbirefringent terms ap- 
pearing only in {cp)'^'^ from the birefringent terms con- 
tributing to and possibly also to {cpY^ . 

We first consider the nonbirefringent case. With zero 





, (d).(OE) 
V-F Jnjm 


n 


3 












1 




1 




2 





2 




3 




1 3 




4 





2 4 




d-2 





2 4 ••• 


d-2 


total 




i(d-h l)rf(d- 


1) 



TABLE IV: Summary of the allowed ranges of indices n and 
j for the independent spherical coefficients associated with 
CPT-even nonbirefringent operators. The dimension d is even 
with d > 4, while n < d — 4. The index m satisfies tlic usual 
restrictions —j < m < j, so there are 2j + 1 coefficients for 
each j. For a given dimension d, the number of coefiicients is 
given in the last row. 



leading-order birefringence, Lorentz violation associated 
with the quadratic action for electrodynamics is com- 
pletely characterized by {cp)^" via 

+ r,''^{cprP-v''''{cpD. (58) 

The symmetry constraints (55) applied to this expression 
yield the equation = d^''{cp)'^^'^ in the nonbirefringent 
case. The form of this constraint suggests solving via a 
scalar potential ^p, where 



{cp)"'' d'^d'^^p, 



(59) 



where the derivatives act in momentum space. This re- 
sult holds in the nonbirefringent case but is false in gen- 
eral. Expanding the potential ^p in spherical harmonics, 



dnjm 



we arrive at the minimal set (c^'*)i°^^ of nonbirefringent 
spherical coefficients for Lorentz violation. This expan- 
sion reveals that no i?-type spherical coefficients appear 
in the absence of leading-order birefringence. The in- 
dex ranges and component counting for the coefficients 
(^F^)ijm ^'"^ summarized in Table IV. 

Using the expansion (60), the coefficients (cp)^'^ can be 

expressed in terms of the spherical coefficients {cp^)n^- 
We can also find the contributions to the matrix compo- 
nents in the expansions (52), (53), and (54) in the absence 
of birefringence. Some of these components vanish and 



(60) 
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n 


3 







1 
2 
3 
4 

d-4 


2 

1 3 
2 4 

13 5 
2 4 6 

2 ■■■ d-2 


2 

1 3 

2 4 
13 5 

1 3 ••• d-3 


2 

3 

2 4 

2 4 ••• d-4 


total 


i(d=^-d-30) 


i(d-4)(d2+d + 3) 


i(d-4)(d2 -2d-9) 



TABLE V: Summary of the allowed ranges of indices n and j for the independent -B-type spherical coefficients associated with 
CPT-even birefringent operators. The dimension d is even with d > 4, while n < d — 4. The index m satisfies the usual 
restrictions —j<m<j, so there are 2j + 1 coefficients for each j. For a given dimension d, the number of coefficients of each 
type is given in the last row. 



others become equal, yielding 



id) .(IE) 



0, 



('^DB)njm 
, (d) .(2E) „ 



(Ad) n(OB) 
y^DEJnjm ^ 



(Jd) ^{0E) 
(Ad) 

(d) .(2E) (d) A2E) 

DEJnjm K'^HBJnjm- 



(Jd) 

{i^Sde) 



(61) 



These conditions lead to vanishing ke+ and Kq-, as re- 
quired by nonbirefringence. The nonzero contributions 
are given by 

(Jd) s(QE') -(«+2)(n+3)+3(3 + l) f (d)s(aE) 
y^DEJnjm - 3 y'^F ){n+2)jm 

+ (d_3_„)(d_2-„)(c^^)^°:^\ (62a) 

(Jd.) AOE) -2n{n+2)-j{j+l) ( (d)AOE) ,„r>,. 

. (d) x(2£;) ^ -0-+2)0-+l)j0-l) / (d)^(Oi;) ,^ 

/^(d) \(.0E') -2(»+2)(n+3)+2j(j+l) / (d)x(Og) x 
K'^HBJnjm ^ 3 ){n+2)jm' 

(Jd) AOE) ^ 2n(«+2)+j(j-+l) / (rf)x(Og) 

K'^HBlnjm ^ 3 ^"^F ^(n+2)jm' v"^"^^ 



The above analysis completely characterizes the nonbire- 
fringent Lorentz-violating operators associated with the 
quadratic action for electrodynamics. 

In the presence of birefringence, three additional in- 
dependent sets of E-type coefficients appear. We can 
therefore seek four mutually independent sets of spherical 
coefficients, with one given by (c^'')^*^^ and the remain- 
ing three covering the birefringent portion of coefficient 
space. It turns out that the three new sets of coeffi- 
cients for birefringence have spin weights zero, one, and 

two, and we denote them by (fc^'')i°^\ (^p )nj^^ ^^'^ 

{kp The index ranges and component counting 

for these three sets arc siimmarizcd in Table V. 

Imposing the symmetry conditions (56), a lengthy cal- 
culation yields explicit expressions in terms of (/c^^)^°^\ 

(^i^^)nj-m; ^^'^ (^F^)rij>Ti ^'^^ the twelve sets of codepen- 
dent i?-type coefficients appearing in the expansions (52), 
(53), and (54). We find contributions to the codependent 
coefficients arising from the trace components oi Kde and 
khb given by 

(Jd) JOE') (n+2){n+3)-jij+l) (Jd)JOE) 
\'^DE)njm — 3 ^^F Jnjm 

+ (d-3-n)id-2-n){kf)l%^ 

_ 2{d-3-n)((n+l)(n+2) + (j+2)U-l)) (-J'^) .(IE) 

3 ^'^F ){n-l)]m 

2{d-i-n){d^2-~n) (-J.d.).{2E) x 
3 y'^F )(n-2)jm^ \^-^^) 



{l^HBlnjm ^ -("+l)jO + l)(,Cp. j(„+2)im' ("^S) 

(Jd) A2E) (j+2){j+l)jU-l) ((d). jOE) X 

K'^HBJnjm =i 2 v''F )(n+2)jm' \0^n) 

(Jd) AOE') ^ , . „ , , ^,,(Jd).(OE) .„„.x 

{I^DBlnjm ^ (d-3-n)(n+l)[Cp )(„+i)j„j, (b^l) 

[l^DBJnjm ^ -id-3-n)j(j + l)[Cp )(n+l)jm- (^^Ji 



/ (d) x(OiJ') ^ -2((n+2)(n+3)-j(j+l)) (Jd).(aE) x 



F Jnjm ' 



The contributions to the remaining codependent coeffi- 
cients in Kde are 



(Jd) JOE) 2n(n+2)+j(j+l) ,j (d).(OE) 
y^DE/njm — 3 V'^F Jnjm 



+ 



2(d-3-n)(4(n+l)(n.+2)+(j+2)(j 



-i))(f'^\(m 

K'^F )(n-l 



+ 



2(d-3-n)(d-2-n 



■[f^F )(n-2)jr. 



(n-l)jm 

(64a) 
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- 2(d-3-n)((n+l)(n+2) + (j+2)(j-l))(n+l)(fcp. )^„_i)^„. 



— ^((i) (2E) 

- 2id-3-n){d-2-n)in+l){kp )(„J2bm' 



(64b) 



(IS) 



y^DE/njm - 2 ^""F )njm 

+ (d_3_„)Q-+2)(j-l)(2(„+2) =+0+2)0-l))(A;^'')|^t^)^.^ 

— ^(c?) (2,E) 

+ {d-3-n)id-2-n)(2in+2f + U+2)U-l)){kp )\n-2)3m^ 

(64c) 

while the contributions to the remaining codependent co- 
efficients in k,HB are 



(Jd) x(OE) 2n{n+2)+j{j+l) (,{d).(OE) 
K'^HBJnjm — 3 '-'''-F )nim 



F )njm 
+ 2{n+3){n+i){kp y^j^, 



K'^HBJnjm ^ -("+1)JU+1)1«F )njm 

-2(n+3f(n+A){kp )^^.J, 



(65a) 



(65b) 



, (d) .(2E) (j+2)(j + l)j(j-l) (,{d)\(01S) 
y'^HBJnjm — 2 V^F 

+ (2(n+3)^-jO+l))(n+3)(n+4)(A;;."')^^^^ (65c) 



fnjni 

■^^d) (2E) 



Finally, the contributions to Kdb are 

, (d) ^(OiS') ^ , . „ _^nCi.('')^(0^) 

{l^DBJnjm ^ (d-3-n)(n+l)(fcp )(„_i)j„ 



- (n+2)(n+3)(n+4)(fcp, );^^.^^ (66a) 



i'^DBJnjm ^ -{d-3-n)j{j + l)[lip )(^n_iy 



{d).(OE) 



+ {n+2f(n+3){n+i)(]^pf^''2^ (66b) 



2 -n(„+2)(„+3)(n+4)(:^?^(l^) 



' njva 



,-p(d) r2E) 
■ 2(d-3-n)(n+2)(n+3){kp )\n-i)jr 



(66c) 



(«Dl3)iS ^ (i+2)a-l)("+2)(n+3)(n+4)(A;i.^^)i'^:^) 

+ 2(d-3-n)(«+2)(n+3)=^(fc^ )(„-l)jm- (66d) 

The above twelve equations completely characterize the 
birefringent CPT-even Lorentz-violating operators asso- 
ciated with the quadratic action for electrodynamics. 

To summarize this subsection, we find that the coef- 
ficients controlling the CPT-even Lorentz-violating op- 
erators separate into one set of nonbirefringent E^-typc 
spherical coefficients, three sets of birefringent E'-type 
spherical coefficients, and two sets of birefringent B-type 
spherical coefficients. The component counting for each 
of these sets can be found in Tables III, IV, and V. 



Adding the totals from each table, we find a net total 
of {d + l)d{d — 3) independent spherical coefficients at 
each even dimension d. As expected, this matches the 
group-theoretic result obtained in Sec. II. Some proper- 
ties of all these coefficients are summarized in Table XVII 
of Sec. VIII. 



IV. SPECIAL MODELS 

For a given observational or experimental procedure, 
particular sensitivity may be achieved to specific combi- 
nations of spherical coefficients. Also, the nine basic sets 
of spherical coefficients contain a large number of inde- 
pendent components, which suggests that some types of 
analyses may be challenging to perform in the general 
context. It is therefore useful to identify special models 
that arc limiting cases and are relevant to specific mea- 
surements. 

In this section, we first present several convenient lim- 
iting cases. Five basic types of model arc considered. 
One is the minimal SME, which restricts attention to 
renormalizable terms. Another is the isotropic limit, in 
which rotational invariancc is preserved in a preferred 
frame. A third is the vacuum model, which encompasses 
the subset of effects detectable in searches using disper- 
sion or birefringence from astrophysical sources. The 
fourth is the nonbirefringent nondispersive model, for 
which the Lorentz-violating photon propagation exhibits 
several parallels to the conventional case. The fifth is the 
vacuum-orthogonal model, which is the complement of 
the vacuum model and contains coefficients that can be 
detected only in laboratory experiments. In a final sub- 
section, we identify particular subsets of the spherical 
SME coefficients that correspond to some special cases 
appearing in the literature. 



A. Minimal SME 

The first limiting case we consider is the pure-photon 
sector of the minimal SME, obtained by restricting at- 
tention to Lorentz-violating operators of renormalizable 
dimension d < A. Numerous measurements of the cor- 
responding coefficients have been performed to date [6], 
with most analyses using cartesian components rather 
than spherical ones. Discussions of the minimal SME can 
be found in the literature [7, 11], so we restrict the fo- 
cus of this subsection to establishing the relationship be- 
tween cartesian and spherical coefficients. Working in a 
fixed inertial frame with cartesian coordinates {t, x, y, z), 
we determine the linear combinations of spherical coeffi- 
cients that correspond to the cartesian coefficients. 

For the Lorentz-violating operators with d = 3, which 
are CPT odd, there are four cartesian components of 
i^Ap)'^ ■ They are linear combinations of the four spher- 



ical coefficients (fc^p')ooo' and (^^f)oito' where m 



„(3) ^{1B) 
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(ui3) \(0S) 
y'^AFJOOO 


(1,(3) \iiB) 


y'^AFJOW 


(1,(3) \(1B) 

v'^AF^l(-l) 




, (•■!) 
(V")IJIJ 










1 



































(kApY 








-fi 






TABLE VI: Matrix elements relating cartesian to spherical 
coefficients in the CPT-odd part of the photon sector in the 
minimal SME. Note that an overall factor of •y/l/47r is omit- 
ted. 



The elements of the matrix S connecting the parity-even 
coefficients kde and khb and the spherical coefficients 
are given in Table VII. The elements of 5* for the parity- 
odd coefficients kdb are given in Table VIII. These ta- 
bles also relate the 20 spherical coefficients with d = 4 
to the relevant 20 vacuum coefficients cfijj^j^, ^{E)jm' 
^{B)jm^ which are defined in Eqs. (84a), (84b), and (84c) 
in the context of the vacuum model presented in Sec. 
IV C. 

In the photon sector of the minimal SME, the non- 
birefringent operators arc controlled by the ten spherical 
coefficients (c^^)|j"^\ as shown in Sec. IIIB. Nine combi- 
nations of these have been measured in laboratory-based 
experiments. The remaining coefficient combination, 



0,±1. In terms of a 4 x 4 matrix S, we can write the 
relationship generically as 

: S • /Cgpherical) 



(67) 



where /Ccartesian and /Cspherical are four-dimensional col- 
umn matrices containing the cartesian and spherical com- 
ponents, respectively. Table VI gives the elements of the 
matrix S. This table also displays the correspondence be- 
tween the four spherical coefficients (fc^;-)oof \ (^!iF)oim 



and the four vacuum coefficients k 



(3) 



which arc de- 



fined in Eq. (84d) below in the context of the vacuum 
model discussed in Sec. IV C. 

As an example from Table VI, consider the cartesian 
coefficient (kAp)'^- In terms of spherical coefficients, the 
table gives 



(kAFf 



1 



1 



3 /.,(3) 
8 \'^AF)o11 



3 /iL(3) 

K'^AFJoii-l) 



(3) 

{V)ll 



1.(3) 
2 '^{V)l{-1) 



). (68) 



Note that the spherical coefficients are complex in gen- 
eral, but the property (fc^^].)[,if = ~[{k'f^ 



MF^l(-l)J 



con- 



tained in the phase condition (45) ensures a total of four 
real degrees of freedom. Indeed, the components (kAF)'^ 
and (^af)^ correspond, respectively, to the real and 

imaginary parts of the spherical coefficient {k'"AF)ari'^ ■ 

For the renormalizable CPT-even Lorentz-violating 
operators, all of which have (1 = 4, the operator kp con- 
tains 20 real independent constants. From Tables III, 
IV, and V, we see that these include two scalar singlets 
(cf')ooo^ and (c^^)^oo^ one triplet (c^^)fi^\ and three 

quintuplets (0^^)22™, (^F^)22m\ (4''^)22m- lu tcrms of 
the k matrices of Eq. (18), the two singlets can be viewed 
as corresponding to the two trace components, while the 
triplet corresponds to the antisymmetric part of kdb- 
The three quintuplets match to the traceless parts of kde 
and Khb and the symmetric part of Kdb- 

We can again write the relationship between cartesian 
and spherical coefficients using the generic form (67). 



k. ^ - ^^f^(4)^(0i^) 

*'■+ ~ \/4¥ F '000 



3(4^^)^)], (69) 



is Lorentz invariant and leads to a rescaling of the elec- 
tric and magnetic fields. This combination is therefore 
typically taken to vanish, which implies the condition 



/J4)^(0F) _ (4).(0E) 
y^F JOOO ~ '^y'^F )200 ■ 



(70) 



This property has been assumed in constructing the en- 
tries for Table VII. 

Using the matrix elements in Tables VI, VII, and VIII, 
it is comparatively straightforward to convert between 
cartesian and spherical representations of the photon- 
sector coefficients in the minimal SME. Note, how- 
ever, that the spherical coefficients represent angular- 
momentum eigcnstatcs and therefore have simpler rota- 
tional properties. The behavior of the spherical coeffi- 
cients under rotations is discussed in Sec. V. 



B. Isotropic models 

In any given inertial frame, a small subset of Lorentz- 
violating operators preserve rotational invariance. Re- 
stricting attention to these operators defines an interest- 
ing limiting case of the general theory. In these mod- 
els, sometimes called 'fried-chicken' models due to their 
popularity and simplicity, the isotropic inertial frame 
must be specified because observer boosts to other frames 
destroy the rotational invariance. One natural choice 
for the preferred frame is the frame of the cosmic mi- 
crowave background (CMB), but other choices are pos- 
sible. Note, however, that isotropy in the CMB frame 
implies anisotropy in a Sun-centered frame and in labo- 
ratory experiments. 

Within our analysis, imposing rotational invariance is 
immediate. The general isotropic model is obtained by 
imposing the condition that all spherical coefficients van- 
ish in the preferred frame, except those with j = 0. This 
condition drastically reduces the number of available co- 
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TABLE VII: Matrix elements relating cartesian to spherical coefficients in the CPT-cvcn and parity-even part of the photon 
sector in the minimal SME. We have assumed (c^'')ooo'' = 3(c^'')2oo'- Note that an overall factor of yT/lTr is omitted. 
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TABLE VIII: Matrix elements relating cartesian to spheric al coef ficients in the CPT-even and parity-odd part of the photon 
sector in the minimal SME. Note that an overall factor of -^/l/Air is omitted. 



efficients. The only nonzero coefficients are: 

\<^p )n — y^p )nOa ' 
y^F in — y^F InOQ ' 
{k^AF)n = {k'AF)nOO ■ 

Following standard convention [38], these isotropic co- 
efficients are identified by a ring diacritic. Note that 



the coefficient {c}p )n controls isotropic nonbircfringcnt 
Lorentz-violating operators in the preferred frame, while 
the others control leading-order birefringent effects. 

In the general isotropic model, the nonzero components 
of Uaf are 

, .d—Z—n„n 

ihF)o = E /r- ^^afU (72a) 
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dn 



477 n + 2 



(72b) 



For the k matrices, wo find tlic expressions 



{KDE)rr = ^ ' 



dn 



4tt 



[- (n+l)(n + 2) 



+ {d-3- n){d-2 - n) 

X ((4'^)n + (fcl.'^)(„-2))] , (73a) 



dn 

X ((4'^)(n+2) - (4'^0 

+ ((i-3-n)(rf-2-n) 

x{{^f)n + (kf\n-2))l (73b) 



{KHB)rr = ^ 



dn 



47r 



(-2)(n + 2) 



X ((4^0(n+2) + (fci'')n), (73c 



(AHB)+- = i:^Sr^(-i)(" + 2)^ 



X ((4))(„+2) + (^^'0n), (73d) 



{kdb)±t = {Ti){d - 3 - n)(n + 1) 



dn 



An 



x((4^0(n+i) + (fci.'^n-i)). (73e) 



The index ranges and component counting for the 
isotropic coefficients can be determined from the J = 
columns of Tables I, III, IV, and V. For a given even di- 
mension d, there are d/2 coefficients (4')" ('^~4)/2 
coefficients {kp^)n, for a total of (d — 2) independent 
components. For odd d, there are {d — l)/2 coefficients 

(^af)"- ^ summary of the properties of these coefficients 
is provided as part of Table XVIII in Sec. VIII. 



C. Vacuum models 

For observations of light from distant sources, far-field 
solutions apply. The electromagnetic fields are well ap- 
proximated by vacuum plane waves with lu = p. Al- 
though Lorentz-violating operators typically generate a 
nontrivial dispersion relation with u) ^ p, we can im- 
pose a; = p in the operators kp and kAp to obtain 
leading-order results. The resulting vacuum coefficients 



for Lorentz violation comprise the set of linear combina- 
tions of spherical coefficients relevant for photon prop- 
agation in the vacuum. They are of direct interest for 
studies of light from astrophysical sources, and they can 
also be important in Earth-based tests. 

Imposing u) = pin the expansions (9) and (10) for Jcaf 
and kp dramatically simplifies the spherical- harmonic 
structure. It turns out that the relevant combinations 
of coefficients are those associated with the Stokes pa- 
rameters of the eigenmodes for vacuum photon propaga- 
tion, because these eigenmodes characterize the effects 
of Lorentz violation on the overall photon polarization 
[11]. To find the eigenmodes, we first need plane- wave 
solutions of the modified Maxwell equations (11). 

At leading order, the covariant dispersion relation Eq. 
(30) becomes 

-A{p^{kAF)^.f^O. 



anjuPffPSP'^P'' 



(74) 



The two solutions at leadi 



nip; orrler can 



he written as 



pO ~ (1 - ^0 ± ^(^1)2 + (^2)2 + (^3)2 



(75) 



where the three rotationally invariant dimensionless com- 
binations 

= '^{cprp^pjuj' 



(^1)2 + (.,2)2 = lix^rP''\x^)^^^,p^psp^p-^/i0\ 



^ = -p>'{kAF)f./iO^ 



(76) 



contain coefficients for nonbirefringent CPT-even, bire- 
fringent CPT-even, and birefringent CPT-odd effects, re- 
spectively. The three combinations <;^, "j^, turn out to 
determine the Stokes parameters of the two eigenmodes. 
The momentum structure of Eq. (76) reveals that the 
vacuum coefficients for odd d arc simply obtained as the 
totally symmetric and traceless part of (A;^^)'*"i'''"(<^-3) . 
For even d, a similar result holds for the vacuum coeffi- 
cients {cp Y'^ controlling nonbirefringent effects. 

The procedure for finding the Stokes parameters 
(s^, s^, s^) = {Q, U, V) associated with the two solutions 
in Eq. (75) follows the same basic steps as in Rcf. [11]. 
Using the plane- wave equation (22), we find a solution 
with unique polarization for each of the signs in Eq. (75). 
It is advantageous to work in the temporal gauge, for 
which Aq = and E oc A. We can then construct the 
specific Stokes parameters for each eigenmode. Adopting 
the orthonormal spherical coordinate system described in 
Appendix A 2, we define a general Stokes vector s and 
its corresponding three Stokes parameters by 



s = 



I 



\ 



(77) 



This Stokes vector completely characterizes the polariza- 
tion of an electromagnetic wave. It is also useful to define 



22 



a fourth Stokes parameter: 

5° = + 1^1', 



(78) 



which corresponds to the intensity /. 

Due to their orthogonaUty, the specific Stokes vectors 
for the two plane- wave solutions differ by a sign. It there- 
fore suffices to construct a single specific Stokes vector <j, 
which can be associated with the faster of the two bire- 
fringent modes. In spherical coordinates, this vector is 
given by 



-i((«e+r-(/Ce+)'^^)-(«o-r 
-(Ke+r + KlKo-r-C/io-)'^*) 
-{{kAFf-{kAFr)/uJ 



(79) 

where we have adopted a convenient normalization. At 
leading order, the components of <; are consistent with the 
combinations (76). In spherical coordinates, the fourth 
combination is 



= hi^DE - kHB + ikoB + ikHE) + - 



(80) 



at leading order. This combination is analogous to the 
scalar Stokes parameter s°. 

The combinations and c,^ are invariant undcT rota- 
tions about the direction p of the photon 3-momentuni. 
like the general Stokes parameters. They therefore have 
zero spin weight. The remaining combinations and 

transform under rotations as a rank-2 traceless tensor 
in the tangent space of the sphere and so are combina- 
tions of components with spin weight ±2. As a result, 
the expansion in spin- weighted spherical harmonics of the 
Stokes parameters for Lorentz- violating photon propaga- 
tion takes the form [22, 23] 

djm 



djm 



{B)jm.l ' 

(81) 



djm 



This result is independent of the index n because at lead- 
ing order we can take w = y» in Eq. (76). We remark that 

the notation c|^|^.^ adopted here for one of the sets of 
vacinim coefficients differs from that of Ref. [23], where 
the quantity 



is used instead. The improved present notation c 



(82) 



re- 



fiects the absence of birefringence from the corresponding 

Lorentz-violating operators. Also, the factors of (—1)-' in 
Eq. (81) have been introduced for convenience and to 



match the definitions in Ref. [23]. The point is that vac- 
uum models are well suited for studies involving radiation 
from astrophysical sources, for which the source direction 
n is normally specified rather than the propagation di- 
rection p = —h. For these studies, it is therefore more 
natural to work with spin weight defined with respect to 
h instead of p. The correspondence 



(- 



s-^ jm 



(83) 



which follows from Eq. (A7), can be used to transform 
between the two pictures. 

The net result of the above discussion is that the vac- 
uum behavior is controlled by four sets of vacuum co- 
efiicients cj^j^.„, k['^J^.^, and k[^^^^. These are 

related to the general coefficients for Lorentz violation by 
taking the limit lu = p in the various expansions provided 
in Sec. III. Via this limit, we find the results 



il)jm ~ 2\ V 2\'^DE)njm I 2\'^HB)njm 
n 

, (d) ^(0i^') _ . (d) .{OE') . (d) .{OE')^ , . 

^y'^DE/nim K'^HBJnim ^ ^K'^DBJuim ) i 1,°^^; 



'^{E)jm 



id) 

(B)jm 



(('*DB)iiJ + i^''HB)n]m + ^(Kijb)!^^), (84b) 



((Ad:) x(2B) 
W^DEJnjm 



0-^2)! 

id) n( 
/r 



y^HBJnjm 



^('^db) 



DB Jnjm ) 



(84c) 



M -Sr-(_Uj+l( d f,{d).{OB) 1 .,(<^)^(1B)^ 



(84d) 



As expected, no frequency or wavelength dependence ap- 
pears in these expressions. These equations reveal that 
the vacuum coefficients are linear combinations of the 
general coefficients involving identical d, j,m indices but 
different n values. 



The index ranges for the vacuum coefficients c 



Sd) 



Ad) 



'-(E) jm' "'(B)jm^ ^nd k^^-^.^ Can be found using general 
symmetry arguments or the relations (84). For even d, 



(d) 

(/)jm' 



there are {d — 1)^ coefficients c|j|^.^, {d —1)^ — 4 coeffi- 



(ij)jm' (d-l)2-4 coefficients k^j^^-^, giving a 



'(/)jm 

cients k 

total of 3(<i— 1)^ — 8 independent components associated 
with CPT-even Lorentz violation in vacuum propagation. 

For odd d, the coefficients ^(y)^^ for CPT-odd Lorentz 

violation have [d — 1)^ independent components. Some 
properties of these coefficients are summarized as part of 
Table XVIII in Sec. VIII. 

The phcnomcnological effects controlled by each set of 



vacuum coefficients are different. The coefficients fc, 

dd) 



(d) 

{E)j7n 



and k(^B)jm associated with CPT-even operators that 
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lead to birefringence, with the propagating modes being 



linearly polarized. The coeiEcients 



id) 

(V)jm 



control CPT- 



odd birefringence, and the corresponding eigenmodes are 
circularly polarized. In contrast, the CPT-even opera- 
tors associated with the coefficients c^rs - are nonbire- 
fringent. These physical differences suggest it may be of 
value to introduce special vacuum models as particular 
limiting cases. For example, a nonbirefringent vacuum 

model involving only nonzero coefficients c|jj^-^ can be 
countenanced. This model has Lorentz-violating oper- 
ators only in even dimensions d = 4, 6, 8, . . ., with the 
number of independent coefficients being 9, 25, 49 
respectively. Another model of interest can be obtained 
as a hybrid of the general vacuum model and the isotropic 
limit, by introducing a preferred frame in which atten- 
tion is restricted only to j = coefficients. This vac- 



uum isotropic model involves only the coefficients c 



and k 



id) 



(d) 
(1)00 



''(y)oo- therefore contains exactly one Lorentz- 
violating operator for each value of d. 



D. Camouflage models 

The vacuum models considered in the previous sub- 
section are well suited for astrophysical observations. In 
this and the following subsections, we consider instead 
the subset of spherical coefficients that are challenging 
to detect via studies of astrophysical birefringence or dis- 
persion but that produce observable effects in laboratory 
experiments. We begin here by focusing on Lorentz- 
violating models without leading-order birefringent or 
vacuum-dispersive effects, which we call camouflage mod- 
els. The next subsection generalizes this treatment to 
identify explicitly all vacuum-orthogonal coefficients, in- 
cluding ones that control non-vacuum birefringence. The 
key idea motivating these constructions is that studies of 
astrophysical birefringence yield among the best sensitiv- 
ities to Lorentz violation in any sector, due to the accu- 
mulation of polarization changes that occur over the cos- 
mological propagation times. This also applies to studies 
of astrophysical dispersion, albeit typically at lesser sen- 
sitivity. As a result, searches for Lorentz violation can 
naturally be split into ones sensitive to operators gen- 
erating vacuum birefringence or vacuum dispersion and 
ones with sensitivity to other operators. 

We begin by restricting attention to operators without 
leading-order birefringence, which are associated with the 
spherical coefficients {cp^)njWn- Some key results for this 
case arc discussed in Sec. Ill B above. Following Eq. (58), 
the Lagrange density reduces to 



(85) 



Note that CPT invariance holds because kAF = 0, and 
that only even- dimensional operators with E-type parity 
contribute. 



Using Eqs. (52) and (62), we can write the k matrices 
in terms of the coefficients {cp^)nj^ by taking the equal- 
ity in (62). However, it is often more convenient to use 
the explicit expansion (59) and (60) via the generating 
function $. The results are summarized as 

dnjm 

x{d-2-n){d-3-n){cf)^°f2, (86a) 



,d—4—n^n 



dnjm 



X {{n + 2)(n + 3) - jU + I)) (0'^'%%^^, (86b) 

(Ke-)rr = -2(/«e+) + - = Yl <^'^~^~>'' oY^m{p) 

dnjm 

-{2n{n+2)+3(3 + l)) ,(d)UOE) ^ 
^ 3 ){n+2)jm^ 

{ke-)±r = {ke-)r± = ^ '^'^"^'"p" ±lYjmip) 
dnjm 



x{±l){n + l)^m^{cf)l%^,^, (86d) 



(«e-)±± = ^''~'~"P'' ±2Yjm{p) 

dnjm 



-1/ (3+2)! ( {d).(OE) . 

{ko+)+- = -{ko+)-+ = Y ^'-'->"oyim(p) 

dnjm 

X (-i)(d - 3 - n)(n + 1) {c^^X%^, (86f) 

iko+)±r = -iko+)r± = Y '^^~^~''P'' ii^jmip) 

dnjm 

Recall that the matrix ke- is traceless and symmetric, 
while ko+ is antisymmetric. For no leading-order bire- 
fringence, the relations 



koE = ke- + ktr- + ktr+, 
kuB = —ke- — ktr- + ktr+, 

koB = ko+ 



(87) 



also hold. 

The next step is to consider dispersive effects. For 
dimensions rf > 4, only a subset of independent combi- 
nations of the coefficients (c^'')|j"^^ are associated with 
leading-order vacuum dispersion. The nondispersive op- 
erators are precisely the ones of interest that define the 
camouflage models. At leading order, the condition that 
ensures no vacuum dispersion is 



0=Pi^P,.{cfT''\p^=o- 



(88) 
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V^F Jnjm 


n 


3 










1 


1 




2 
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3 


1 


3 


4 





2 4 


d-4 





2 4 ••• rf-4 


total 


w 


- 2)(d- 3) 



TABLE IX: Summary of the allowed ranges of indices n and j 
for the independent camouflage coefficients. The dimension d 
is oven with d> 4, while n < d — 4. The index m satisfies the 
usual restrictions —j<m<j, so there are 2j + 1 coefficients 
for each j. For a given dimension d, the total number of 
coefficients is given in the last row. 



This is satisfied if the generating function $f in Eq. (60) 
is of the form = P^^f- We can therefore define the 
camouflage coefficients (c'^^)^°^^ through the expansion 



E 

dnjm 



d-i- 



(v) 



This result leads to the comparatively simple relation 

( id).{OE) _ (-,{d)-.{OE) _ (-,(d).(OE) 

Jnjm — y'^F Jnim F )(n-2)im ) 



and k^v)jm- "^^^ number of independent vacuum coeffi- 
cients grows as for large d, whereas the total number of 
coefficients grows as d^. Consequently, the vacuum coefii- 
cionts comprise a comparatively small portion of the total 
coefficient space. In this subsection, we construct a com- 
plete set of independent coefficients spanning the comple- 
mentary part of the coefficient space. We refer to these as 
the vacuum-orthogonal coefficients. At leading order, the 
corresponding vacuum-orthogonal operators induce nei- 
ther dispersion nor birefringence for photon propagation 
in the vacuum. Instead, these Lorentz-violating effects 
can become relevant whenever the boundary conditions 
or the properties of macroscopic media differ from those 
for vacuum propagation. 

To extract the vacuum-orthogonal coefficients, we be- 



gin by considering the E-type vacuum coefficients 
and k, 



id) 



^(E)jm- Expanding Eqs. (84a) and (84b) in terms 
of the general spherical coefficients, we find that each 

is a combination of (c^^)i°^', {k'^p^)%^, (^p and 
{kp Careful consideration shows that the restric- 



Jnjni ■ 



tion k, 



\E)jm = can be achieved by writing {k'^p^)^^^2. 
as combinations of {k^p^)^^^ and {k^p^)^^^ and a new 

.(d). 



set of vacuum-orthogonal coefficients (k"p' )^°^' . The net 

Sd)^{aE) fT^KioE) 



result is the replacement of {k'^p')nj^ with (^p')^^ m 



the nine sets of general coefficients, implemented via the 
substitution 



in the limit of no leading-order birefringence or vacuum 
dispersion. 

The index ranges for the camoufiage coefficients 



{^p )njm are given in Table IX. For each even di- 
mension d, there are {d — l)((i — 2)((i — 3)/6 indepen- 
dent components. They represent combinations of the 

{d + l)d{d — l)/6 coefficients {c'p' fn^ thsii are comple- 



mentary to the (d — 1)^ vacuum coefficients c|jj^.^ intro- 
duced in Sec. IV C. Note that a subset of the camouflage 
operators are rotation invariant, so a hybrid camouflage 
isotropic model exists that has {d — 2)/2 independent co- 
efficients for each even d. 

The camouflage coefficients ((?^'')i°m are challenging 
to detect via astrophysical observations of birefringence 
or dispersion because their effects arise only at higher or- 
der. Along with the minimal-SME coefficients {c'p')^^^^, 
they are best sought via alternative methods such as lab- 
oratory experiments. A class of experiments sensitive to 
the effects of these coefficients is discussed in Sec. VII. 



E. Vacuum-orthogonal models 



,Ad).(OE) .-pid) (OE) Md).(OE) 
y'^F Jnjm ~^ y'^F Jnjm y'^F J (n+^^ 



(n+2)jm 
f-p'-d\ilE) ,-p(d) (^2E) 



where the numerical factors ai and 02 are 

(d-5-n)(2(n-|-4)^-Kj-|-2)0-l))+2(n-|-3)(n-|-4)(n-|-5) 



ai 



iU+i) 



a2 = (j+2)0 + l)j0-l) [(d-5-n)id-4-n){2{n+i)^+{j+2){j-l)) 
+4{d-4-n){n+3)(n+4)'+{2{n+3)'-j{j+l)){n+3){n+4)] . 

(92) 

The result (91) represents combinations of i?-type coeffi- 
cients that are complementary to the vacuum coefficients 
^(£)jm- "^^^ index ranges and counting for the vacuum- 
orthogonal coefficients (fc^^)^°^^ are given in Table X. 
The restriction c['t} . = 

{I)jm 

the coefficients (c^^)i°^\ 



yields a similar result for 



In Sec. IV C, we identified the subset of spherical co- 
efficients relevant for photon propagation in the vacuum. 



.id) 



Sd) 

'■{I)jm^ '^{E)jm^ '^{B)jm^ 



.id) 



y'^F Jnjm 



.-,(d)UOE) _ /-,(c^)^(OS) _ , /, (d)x(OE) 
y^F Jnjm y^F J(n-2)jm "1V'*'F J(n-2)]m 



+ b2{kp )^„_;)^.„ 



, , f-p^d.\i2E) 



(93) 
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- l)(d- 2)(d- 3) - 1 



TABLE X: Summary of the allowed ranges of indices n and j 

for the £-type vacuum-orthogonal coefficients {kp )„^^ as- 
sociated with CPT-even nonbirefringent operators. The di- 
mension d is even with d > 4, while n < d — 4. The index m 
satisfies the usual restrictions —j<m<j, so there are 2j + 1 
coefficients for each j. For a given dimension d, the number 
of coefficients is given in the last row. 




1 
2 
3 
4 

!-4 



total 



1 

2 

1 3 
2 4 



d-4 



id(d-2)(d-4) 



TABLE XL Summary of the allowed ranges of indices n and j 

for the S-type vacuum-orthogonal coefficients (kp )njrn 
sociated with CPT-even nonbirefringent operators. The di- 
mension d is oven with d > 4, while n < d — 4. The index m 
satisfies the usual restrictions — j < m < j, so there are 2j + 1 
coefficients for each j. For a given dimension d, the number 
of coefficients is given in the last row. 



where the numerical factors bi, 62, and 63 are 



7 _ {d-3-n){d-2-n)+2{d-2-n)n+n{n+l)-j{j+l) 

"1 - (d-3)(d-2) ' 

1, _ (ci-3-7t)(2(7t+l)(w+2) + (j+2)(j-l))+2(7t+l)(n+2)(7t+3) 

"2 — (<i-3)(d-2) ' 



^3 = 



_ {d-5-n){d-4-n)-{n+3){n+4) 



(d-3)(d-2) 



The coefBcients (c'i?^)i°m appearing in Eq. (93) are the 
camouflage coefficients introduced in the previous sub- 
section. Their index ranges and counting are given in 
Table IX. Note that appears in Eq. (93), so 

the expression (91) is needed to fully reduce icp^)nj^- 
Note also that Eq. (93) generalizes Eq. (90) to include co- 
efficients for birefringent operators that leave unaffected 
vacuum propagation at loading order. This shovirs explic- 
itly that these operators contribute both to birefringent 
and to nonbirefringent effects, an effect discussed in Sec. 
IIIB. 

Determining the combinations of spherical coefficients 
that are complementary to the vacuum coefficient k'"^B)jm 
is more straightforward. We can write 



k 



id) 

{B)jr. 



njm 



+ (i + 2)(i-l)(4^))J^)). (95) 

Taking into account the index ranges of {k^p^)^^^ and 

(^f')!^™' fi'^*^ t^^* ^ suitable set of coefficients com- 
patible with the restriction 

(d)s(lS) 



id) 

{B)jr. 



is obtained by 



of vacuum- 



replacing (fcV)rijm ^it^ ^ ^'^^ 

orthogonal coefficients. The combination leading to van- 



set (^F^)nfnl 



ishing k 



(d) 

(B)jr. 



(UdhilB) 
y^F /njm 



is given by 

)njm ){n+2)jr 



{d-2) (jS'^Kim 

0+2)0-1) I J(„+l)jm- 



(96) 



(94^ The index ranges and counting for the coefficients 
(^F given in Table XI. Note that the coeffi- 



cients {icp and (^p )^^^^ are absent from Eq. (93), 
so the corresponding operators have no nonbirefringent 
effects. 

Finally, we construct the combinations of spherical co- 
efflcients that are complementary to Although 



the vacuum coefficients k^(y-^j„^ have the comparatively 
simple form (84d), finding combinations that cover the 



coefficient space under the restriction k 



(d) 



involves 



some calculation. It turns out to involve two new sets 

(^AF)nj^, (^AF)nj2 of vacuum-orthogoual coefficients, 
which appear via the substitutions 

(M .(OB) (d-2-n)(»+3) //7:('^) x(OB) _ M^) (oB) n 
y'^AFJnjm ~^ d{d-2-n+j) {y^AFJnjm \'^AF)(n-2)jm) 



(M^i^B) 
K'^AFJnjm 



(d-) .(IB) 
(n-1] 

W .(OB) 



n+l\'^AF)(n-l)j 
-n+i \^ 



d—3—n+j 



"AF) 



+ i)jr. 



d (l.^'^UilB) 
AF 'njm ' 



(97) 

\'^AF)(n-l)jmJ 

(98) 



The index ranges and the numbers of indepen- 
dent components for the vacuum-orthogonal coefficients 



-.(OB) 



.(IS) 



ikAF)nj-ni ^nd {k^p)^^^ arc shown in Table XII. 

To summarize, in this subsection we have completed 
the decomposition of the original nine sets of general 
spherical coefficients into those that generate leading- 
order effects in the vacuum propagation of light and those 



26 













(11?) 
"' 


n 




j 



1 
2 
3 
4 




1 
1 



1 




2 

3 

2 4 




1 

o 
z 
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-l)(d- 


-2)(d-3) 


i(d + l)(d- 


-l)(d-3) 



TABLE XII: Summary of the allowed ranges of indices n 
and i for the vacuum-orthogonal coefHcients {kj^p)^^^ and 

{kAF)njm associatsd with CPT-odd nonbirefringent opera- 
tors. The dimension d is even with d > 4, while n < d — 4. 

The index m satisfies the usual restrictions —j < m < j, so 
there are 2j + 1 coefficients for each j. For a given dimension 
d, the number of coefficients is given in the last row. 



F. Connections to other formalisms 

Several specialized models involving particular 
Lorcntz-violating photon operators with d > 4 have 
been considered in the literature. The generality of the 
SME implies that any realistic model for the photon 
propagator compatible with standard field theory is 
encompassed via special values of the coefficients kp 
and kAF- In this subsection, we outline some of these 
connections. For a selection of specific models defined 
via field theory, we provide explicit limiting values of the 
SME coefficients for Lorentz violation that reproduce 
the models. We also offer here some remarks about 
the relationship between the SME and the kinematical 
approach to Lorentz violation, which is based on altering 
the transformation laws. Comments on the links between 
the photon sector of the SME and Lorentz-violating 
modifications of the photon dispersion relation outside 
the context of standard field theory can be found in Sec. 
VIA. 



1. Field-theoretic models 



that comprise the complementary subset. The four sets 



of vacuum coefficients are c 



k 



{d) 

{V)jm- 



(d) 



k 



They are discussed in Sec. IV C. The com- 



plement consists of the nine reduced sets of vacuum- 
orthogonal coefficients (si''^)i"f„', (A?^ )i°f„', {Hp 

(IB) 

'ijm ' 



M'i).(2E) f-pWjlB) M'^).(2B) ,T^'^K{OB) /7;('^)x(l 
y'^F )njm' y'^F Jnjm ' y'^F )njm' y'^AFJnjm' {'^AFlnj 
— iC*^) (IE) 

and {kAFj^jm- Note that these are nonzero only for 
d > 4. This decomposition is summarized as part of 
Table XVIII in Sec. VIII. 



j(<i)x(0i5) 



dis- 



Except for the camouflage coefficients (c?_p Onjm 
cussed in the previous subsection, all the vacuum- 
orthogonal coefficients control birefringence effects that 
cannot be detected at leading order via vacuum propaga- 
tion. To measure these coefficients, alternative methods 
such as laboratory experiments are therefore desirable. 
In contrast, all four sets of vacuum coefficients are de- 
tectable in astrophysical tests involving birefringence or 
dispersion, with the exception of the special d = 4 coef- 



ficients c|j|^^ . The latter have been extensively studied 



through a variety of different methods [12-17]. 

Note also that various hybrid models involving the 
vacimm-orthogonal coefficients can be countenanced. 
For example, a general vacuum-orthogonal isotropic 
model is obtained upon further restricting attention 

to the isotropic coefficients (c?^^) 



{"'F jriin 1 ana 



'F )nOO > y^F JnOO ' 

(^AF)nm- model has {d — 2)/2 nonbirefringent 

and (d — 4)/2 birefringent operators for even d, along 
with [d — 3)/2 birefringent operators for odd d. 



We begin by discussing models defined via a Lagrange 
density in field theory. Several spc;cialized models exist 
that involve photon fields with a small number of specific 
Lorentz-violating operators of mass dimensions d = 5 and 
in some cases also d = 6. We provide here brief comments 
identifying the match between these models and the SME 
coefficients for Lorentz violation. 

One such model is presented by Gambini and PuUin 
[60]. Lorentz violation in this model is controlled by the 
parameter x^p- To see the connection to the SME, we 
can make the field redefinition 



E + 2xlpV X E. 



(99) 



This model then is equivalent to taking the special limit 
of the SME with the nonzero coefficients being (fc^^)^-''^ 
and (fc^^V'"^^ given by 



"AF) 



2x1 



(100) 

This model involves only isotropic Lorentz violations, so 
it is a special limit of the isotropic models discussed in 
Sec. IV B. In terms of the coefficients for Lorentz viola- 
tion defined in Eq. (71), the match 



ikfF)2 



2V47r xlp, 



At: x'lp/5 



(101) 



provides a complete specification of this model within the 
SME. It involves three nonzero isotropic coefficients con- 
taining a single degree of freedom. A nonlinear general- 
ization of this model is obtained in Ref. [61], in which the 
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propagator component is a special limit of the isotropic 
models discussed in Sec. IV B with three degrees of free- 
dom. 

Another specialized model incorporating photon-sector 
Lorentz violation is introduced by Myers and Pospelov 
[62]. This model has a unit timelike background vector 
n'^ that defines a preferred frame and a parameter ^/Mp 
setting the scale of the Lorentz violation. Only one d = 5 
Lorentz-violating operator affects the photon propagator. 
In terms of SME coefficients, the model is obtained by 
taking the nonzero coefficients for Lorentz violation to be 

{kfp)^'"' = --^(n'^n^n" - ^n^n^'^r)"^"). (102) 

In the preferred frame, the vector takes the form n'' — 
(1,0,0,0). This model can therefore also be defined 
uniquely in terms of isotropic Lorentz violation and un- 
derstood as a special limit of the isotropic models dis- 
cussed in Sec. IV B. It is equivalent to taking as nonzero 
only the coefficient (fc^''^)o, with the specific choice 



3^V47r 
5Mp 



(103) 



A model focusing on a Lagrange density that is explic- 
itly gauge invariant and involves Lorentz-violating oper- 
ators with (i = 5 is considered by Bolokhov and Pospelov 
[63]. By construction, the component of the model rel- 
evant to the photon propagator is restricted to Lorentz 
violation affecting leading-order vacuum propagation. It 
consists of (i = 5 operators governed by a totally sym- 
metric and traceless parameter C^^''. This parameter 
has 16 independent components, corresponding to the 
16 vacuum coefficients ^(y)^^ among the 36 independent 

coefficients {k^^j^)^'^'' for Lorentz violation at d = 5 in 
the SME. In terms of the latter coefficients, the model is 
fixed by 



k 



(5) 

(V)3r. 



symmetric 
traceless 



(104) 



where the correspondence on the left-hand side is de- 
termined by matching Eqs. (76) and (81), and the re- 
striction to the 16 relevant coefficients in (fc^""^)'"^'' is 
obtained by imposing total symmetry and tracelessness. 
The 20 SME coefficients at (i = 5 that are absent 
from this model are the vacuum-orthogonal coefficients 

(^Af) nfrL O^'Af) nfrL and (^af)^' ^^^^^ ™- 

affected leading-order vacuum propagation but nonethe- 
less produce leading-order effects in suitable laboratory 
experiments. 

The series structure of the SME Lagrange density (8) 
implies that models defined via nonpolynomial but ana- 
lytic functionals of field operators can be matched to the 
SME by Taylor expansion. This includes models with 
apparent singularities, provided the expansion is taken 
about a nonsingular point. In the latter case, distinct an- 
alytic continuations may correspond to different values of 



SME coefficients, which may lead to unusual dynamical 

effects. 

As an illustration, consider the class of Lorentz- 
violating models called very special relativity (VSR) [64] , 
in which the Lorentz group is broken explicitly to the 
four-parameter subgroup SIM(2). It is uschd to intro- 
duce the operator A'^'* = /{n ■ d), where is a unit 
null vector that establishes a preferred lightlike frame. It 
follows that the combination N -T involving a tensor field 
operator T is SIM(2) covariant but violates Lorentz co- 
variance. We can construct a generic on-shell linear VSR 
electrodynamics preserving gauge invariance by specify- 
ing field equations of the form 



(105) 



where K^(N) is a nontrival 4- vector function of A'* that 
can also depend on the metric, the Levi-Civita tensor, 
and derivatives. Note that it is problematic to obtain a 
gauge-invariant action for these equations. 

To investigate the correspondence to the SME, it is 
convenient to work in momentum space, where the op- 
erator lA^ = / [n ■ p) is singular whenever n ■ p = 0. 
Matching the momentum-space version of the VSR elec- 
trodynamics (105) to the momentum-space version of the 
SME equations of motion (11) therefore involves expand- 
ing K^{iN)F^"{p) about a nonsingular point Q in mo- 
mentum space. Assigning coordinates Pq to Q, we find 
the expansion converges either fov < n ■ p < 2n ■ pg 
or for 2n ■ pq < n ■ p < 0, corresponding to two differ- 
ent analytic continuations of i^^. Each expansion can be 
matched to a set of SME coefficients, but the two sets of 
coefficients differ. This implies the VSR electrodynamics 
(105) is represented by two different limits of the SME, 
according to the value of the photon momentum relative 
to the preferred lightlike frame. 



2. Robertson-Mansouri-Sexl model 

In a different vein, some authors adopt a kinematical 
approach to Lorentz violation that is based on modifica- 
tions of the transformation laws. An older test model of 
this type that is encompassed by the SME approach is 
the kinematical formalism of Robertson, Mansouri, and 
Sexl (RMS) [65, 66]. This approach assumes that there 
is a preferred universal inertial frame U in which light 
propagates conventionally as measured using a definite 
set of rods and clocks. In other frames, which include 
any inertial frame E relevant for experiment, light can 
behave anisotropically with respect to the boosted rods 
and clocks. The RMS approach assumes that the lengths 
of rods and the ticking rates of clocks are invariant in 
inertial frames related to U by RMS coordinate trans- 
formations T^i,. These arc deformations of the Lorentz 
transformations involving three functions of the boost v, 
conventionally denoted as a{v), b{v), d{v). 

The RMS formalism can be translated into the SME 
framework. Consider first the preferred universal frame 
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U. Since light is conventional in U by definition, the 
Maxwell action must be valid in this frame, so that 

^tuTu = -h'^n'^F.^F,^. (106) 

This Lagrange density is Lorentz invariant, so all Lorentz 
violation in this frame resides in the physics describing 
the chosen rods and clocks. For the photons, we can 
match the RMS formalism to the SME by imposing the 
condition 

/^photon /^photon (^(\^\ 

^SME, U ~ •'-RMS, U' U^' J 

which eliminates all Lorentz- violating effects in the pho- 
ton sector of the SME as seen in U. 

For the RMS rods and clocks, any realistic set is made 
of constituent particles and fields. Since the SME de- 
scribes general Lorentz violation for all particles and 
fields, it follows that the properties of any definite set of 
rods and clocks can be derived from the full SME action, 
at least in principle. This action contains infinitely many 
coefficients for Lorentz violation outside the photon sec- 
tor, so there is plenty of room for anomalous behavior. 
However, to preserve the observed isotropy of light in U 
required by the RMS formalism for any possible choice of 
rods and clocks, only violations of Lorentz invariance in 
U that produce isotropic effects on the rods and clocks 
should be countenanced. We conclude that the match 
between the RMS formalism and the SME in the frame 
U requires restricting the SME to a subset of coefficients 
outside the photon sector. For definiteness in what fol- 
lows, we denote these SME coefficients collectively by 
{k}. Note that {k} can include anisotropic coefficients, 
provided they have no anisotropic effects on the chosen 
rods and clocks. 

Next, consider an experimentally relevant frame E that 
is moving with velocity v relative to U. In the RMS for- 
malism, the properties of light in E are obtained by per- 
forming an RMS coordinate transformation using T^n. 
This leads to a modified Maxwell action with 

(108) 

where 

(5rms)m. = yipa{T-')\(T-^)\ (109) 

is an effective metric that depends on the three functions 
a, b, d. Note that physically different choices for the rods 
and clocks in the frame U imply different invariant RMS 
transformations and hence different o, b, d. 

In contrast, the SME properties of light and the 
boosted rods and clocks in the frame E are obtained by 
performing a particle Lorentz transformation with the 
velocity v. Since 'Cg^g°"u is invariant under particle 
Lorentz transformations, light in the SME must obey 
conventional electrodynamics in the frame E too. How- 
ever, the rods and clocks involve Lorentz-violating oper- 
ators that change under the particle Lorentz transforma- 
tion. This produces a deformation of the standard rods 



and clocks in E that depends on the SME coefficients 
{k}. 

We thus find that the RMS formalism and the SME 
naturally generate two different coordinate systems for 
describing physics in the frame E. To relate the two, we 
can redefine the length and time intervals specified by the 
boosted SME rods and clocks to match numerically those 
of the boosted RMS rods and clocks, and we can choose 
the SME synchronization to match the RMS one. The 
redefinition can be implemented by scaling the spacetime 
coordinates. In the photon sector of the SME, the scaling 
produces an effective metric (.gsME. cs)^^l/ that depends 
on the coefficients {k}. Since measurements made with 
rods and clocks using either coordinate system now agree, 
this metric must match the RMS metric (109) in the 
frame E. This yields the result 

(ffSME, es)t,v{{k}) = {gmAs)iiu{a,b,d), (110) 

which provides a direct correspondence between the SME 
coefficients and the RMS functions. 

The above match shows that the RMS formalism can 
be understood as a special limit of the SME in which nor- 
mal light behavior together with Lorentz violation affect- 
ing rods and clocks isotropically are assumed in the frame 
U. This limit excludes infinitely many Lorentz-violating 
effects. Since the RMS formalism is a special-rclativistic 
test model, the gravitational sector of the SME must also 
be disregarded. Note that the three functions a, b, d can 
be expanded in powers of the velocity to yield a triple 
infinity of constant parameters, which can be absorbed 
into the multiple infinity of coefficients {k}. In the ideal- 
ized case of rods and clocks formed from a scalar particle 
with only one isotropic dimension-zero coefficient k for 
Lorentz violation in the frame U, explicit expressions for 
the three RMS functions in terms of k are given in Sec. 
HI C of Rcf. [11]. 

Another point of interest is that physically different 
rods and clocks are associated with different functions a, 
b, d, and so involve different combinations of the coeffi- 
cients {k}. Within the RMS formalism, measurements of 
a, 6, d in a given experiment cannot meaningfully be com- 
pared to those in another experiment unless physically 
identical rods and clocks are used in both. Note that the 
rods and clocks must also be in the same physical state, 
since state changes in the presence of Lorentz violation 
can deform physical properties. This major disadvantage 
of the RMS formalism is circumvented by the SME. The 
SME coefficients for Lorentz violation are specific to par- 
ticles and interactions and can therefore be reported in 
an experiment-independent way in a conveniently chosen 
frame, which conventionally is taken as the Sun-centered 
frame described in Sec. V. 

To illustrate the above reasoning with an explicit ex- 
ample, we can consider a scenario in which the photons, 
the chosen rods, and the chosen clocks each have proper- 
ties governed by different effective metrics. Following the 
RMS assumptions, the effective metric (^photon, \j)fii' for 
photons in the frame U must be taken as the Minkowski 
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metric, while the effective metric (grod. u)aii> for the rods 
and the effective metric (.9ciock, 13)^1^ for the clocks char- 
acterize the Lorcntz- violating physics in U. We examine 
here a simple model for which 

(^'photon, u)jut' ~ ^jut') 

(ffclock, u)^*!^ = '7^11^ ~ (Cclock, \j)tiu- (HI) 

The quantities (crod, u)/lii/ and (cciock, v)^!^ can be viewed 
as effective coefficients for Lorentz violation defined in the 
frame U. Depending on the nature of the chosen rods and 
clocks, these coefficients can be identified either directly 
with specific SME coefficients or indirectly as suitable 
combinations of SME coefficients. We also assume the 
conditions 

(Crod, u)jO = (Crod, u)ok = 0, 
(Crod, v)jk = ^(Crod, v)ooSjk, 



(Cclock, u)jO = (Cclock, u)ofe = 0, 
(Cclock, u)jfc = ^(Cclock, u)00^jfe, 



(112) 



which ensure isotropic properties of the rods and clocks 

in U, as required by the RMS formalism. For simplicity, 
the coefficients are taken to be independent of spacetime 
position or, equivalently, independent of frequency and 
momentum. This simple model therefore has only two 
degrees of freedom controlling Lorentz violation. 

Suppose for definitencss that the experimentally rele- 
vant frame E is moving with velocity v = {v, 0, 0) relative 
to U. Following the general reasoning above, we can ob- 
tain the SME description of the model in E by performing 
a standard particle Lorentz transformation A from U to 
E. The resulting effective metrics in E are given by ex- 
pressions of the form 



(113) 



In the frame E, light remains conventional but the rods 
and clocks are distorted. To match to the RMS formal- 
ism, we must therefore seek alternative coordinates in E 
in which the rods arc isotropic and of the same length as 
in U, and in which the clocks tick at the original rate in 
U. The appropriate coordinate transformation C leaves 
the origin of E in place but acts on the effective metrics 



5e 



1 — (j-lT ^-lT 



A-i^A-iC-i. (114) 



Since the components {giod)jk determine the rod length 
and the component (gciock)oo determines the clock ticking 
rate, and since both these quantities are assumed invari- 
ant in the RMS formalism, the required transformation 
C is fixed by demanding that 

[C"'^ '^A""^'^(gi.od,u)A ^C^^jjk = (5rod,u)ife) 

[C-i^A-i^(5eiock,u)A-'C-i]oo = (5clock,u)00. (115) 



With the usual definition 7 = l/\/l — v'^, we find that 
the nonzero elements of C take the form 



^0 _ / l - (cciock, u)oo7^(l + 1^^) 



1 — (Cclock, u) 



00 



^1^^ /l + (w u)oo7'(|+i'2) 



1 + |(Crod, u) 



00 



2 — (-^ 3 



1, 
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which amounts to performing different dilations of time 
and of space in the direction of the boost. 

With this information in hand, we can construct the 
RMS transformation T = CK from U to E. This provides 
the direct correspondence between the SME coefficients 
and the RMS functions a, 6, d for this simple model. We 
obtain 



1 / 1 - (cciock, u)oo7^(l + k^"^) 



1 — (Cclock, u)oO 



'1 + (Crod, u)007^(^ + 'V'^) 



1 + 5(Crod, u) 



00 
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where e is the RMS synchronization function in Einstein 
synchronization. Unlike the elementary single-coefficient 
example given in Ref. [11], the present model has a ^ 1/6. 
It also has li = 1, but allowing frequency dependence in 
(crod, u)oo can generate d ^ \ via the mixing of frequency 
with momentum resulting from the A boost and its con- 
sequent effects on C. A frequency or momentum depen- 
dence may arise directly from the inclusion of matter- 
sector operators of nonrenormalizable dimension [44], or 
indirectly from combinations of SME coefficients of renor- 
malizable dimension when the motions of the component 
particles in the rod are incorporated. 

Expanding the above results for a, 6, d to leading or- 
der in v'^ and to leading order in coefficients for Lorentz 
violation yields 



a w 1 -I- av^, a = -i - ^(cdock, u)oo, 
bPil + l3v^, /? = ^ + ^(crod, u)oo, 
d^l + Sv'^, 6 = Q. 



(118) 



The combination of the RMS parameters a. 13, 5 that can 
be tested in Michelson-Morley experiments is known to 
be [66] 



P + S- 



12 i^' 



(119) 



and the combination tested in Kennedy-Thorndike ex- 
periments is 



a 



(3 + I = — ^(Crod, u)00 — ^(Cclock, u)oO, (120) 
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while Ivcs-StilwcU experiments are sensitive to a. Evi- 
dently, two of the three types of experiments are required 
to disentangle all effects, even for this simple model. 
Various special cases can be considered. For example, 
if either {c^od, u)oo or (cdock, u)oo vanishes, then either 
Michelson-Morlcy or Ives-Stilwell experiments have no 
signal. If (crod, u)oo = (cciock, u)oo, all experiments have 
interdependent signals. 

The example verifies that the RMS formalism concerns 
Lorentz-violation residing in the matter sector rather 
than in the photon sector. Note that different choices 
of rods and clocks generically involve different values of 
(crod, u)oo and (cciock, u)oo and therefore are associated 
with distinct Lorentz- violating effects, confirming that 
results from RMS experiments performed using differ- 
ent rods and clocks cannot meaningfully be compared. 
Note also that (crod, u)oo and (cciock, u)oo are defined in 
the frame U, which typically differs from the canonical 
Sun-centered frame S in which SME coefficients are re- 
ported. A conventional particle Lorentz boost from U 
to S can be implemented to identify the relevant coeffi- 
cient combinations in S. Interpreting experiments in the 
RMS formalism requires a choice of U, which is often 
taken as the frame of the cosmic microwave background, 
in which case the relevant boost to S is of order f 0~^ and 
so to a good approximation (ciod. u)oo — (crod. s)tt and 
(cciock, u)oo = (cciock, s)tt- Alternatively, the frame U 
can simply be chosen to be the frame S, since to date no 
compelling evidence for anisotropic Lorentz violation in S 
has been identified. Other choices are also possible. The 
requirement that U be specified to fix the RMS formal- 
ism is a disadvantage that is avoided in the SME, where 
any universal frame U is acceptable and moreover the ex- 
istence of U is unnecessary for interpreting experimental 
data. 



3. Deformed special relativities 

Another kinematical approach involves requiring the 
invariance of all physics under some specified modifica- 
tion of the Lorentz transformations. Recent efforts along 
these lines are generically called deformed special rela- 
tivity (DSR), or in some cases doubly special relativ- 
ity or kappa-deformed relativity [67]. They posit that 
all physics is invariant under a set of deformed nonlin- 
ear Lorentz transformations, usually one that introduces 
a maximum energy scale. A generic DSR model is de- 
fined by replacing the 4-momentum with a modified 
4- momentum tt^, typically with a Planck-scale suppres- 
sion factor for the deformation [68]. The momentum- 
space Lorentz transformations act conventionally on tt^. 
This induces unconventional DSR transformations on p^, 
which are required to leave invariant the physics of the 
model. 

Since the SME contains arbitrary polynomial Lorentz- 

violating operators at all mass dimensions, it must be 
possible to express any nonsingular DSR model involving 



realistic fields in the SME framework. The deformations 
normally are assumed to preserve rotation invariance, in 
which case the match to the SME involves the isotropic 
coefficients discussed in Sec. IV B. Within the context of 
the present work, we can investigate this correspondence 
explicitly in the photon sector. 

Consider a generic DSR model in which the replace- 
ment of the 4-momentum is specified by the nonlinear 
transformation 

p^^7r^ = M^''{p)p, (121) 

acting in momentum space, where M^'^{p) is a nonsin- 
gular matrix that is a deformation of the identity. By 
definition, the Lorentz transformations A^'^ in momen- 
tum space act on tt^ as usual, t^'h = A^'^tTj^. It 
follows that the expressions 

^ = S^y., 5/ = (M-iAM)/ (122) 

specify the nonlinear DSR transformations S*^"^ of p^. 

In the photon sector, we can determine the DSR- 
covariant dispersion relation for the photon by applying 
the replacement (121) to the standard Maxwell disper- 
sion relation p^rj^^Pu = 0. This gives 

Pm(.95sr)'^V = 0, (123) 

where the effective metric 

(5DSR),.. = VpAM-^)f.'{M-%" (124) 

is defined in momentum space. By construction, the dis- 
persion relation (123) is invariant under the DSR trans- 
formations S^" in Eq. (122). 

In the SME context, the DSR dispersion relation (123) 
is recovered as a limiting case of the scalar covariant dis- 
persion relation (30). The match arises in the special 
limit with {kAF)ti = and with the momentum-space 
operator (fc^)^"'"^ given by 

We see that the nonlinear transformation M^''{p) gen- 
erates a subset of the Lorentz-violating operators in the 
SME, typically in the form of an infinite series. 

An alternative derivation yielding the same result can 
be performed at the level of field theory. We can con- 
struct the DSR-invariant modified action for this generic 
model by taking the Maxwell action in momentum space 
and implementing the replacement (121). The transfor- 
mation of the photon field under 5^"^ is taken as 

A^ip) ^ A'^ip') = 5/(p)A.(p), (126) 

corresponding to the replacement — ^ M^^A^. This 
gives the momentum-space Lagrange density 

(127) 
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where the effective metric {gusRjfiv in momentum space 
is given by Eq. (124). By constrnction, the Lagrange 
density (127) is invariant under tlic DSR transforma- 
tions Sf^"^ in Eq. (122). In the SME context, this same 
Lagrange density is obtained as the special limit of the 
SME photon sector of the form 



y^photon 
•'-SME 



^photon 



DSR 



DSR 



where the momentum-space operator (kFY"'"^ is given 
by Eq. (125). This SME construction appUes for a generic 
DSR model, and it confirms that the nonhnear transfor- 
mation Ad^ ^ (p) produces a series of SME coefficients for 
Lorentz violation. 

Although the model (128) transforms nontrivially un- 
der conventional particle Lorentz transformations, ob- 
servable physical effects cannot arise. This is because 
in any frame the inverse momentum replacement can be 
applied to all particles and fields to recover an action 
of the usual Lorentz-invariant form in terms of the mo- 
mentum tt'^, and any experiment naturally identifies this 
momentum as the physical one. The situation in this re- 
spect differs from that of the RMS formalism, where only 
the rods and clocks are assumed invariant under modi- 
fied Lorentz transformations while light behaves differ- 
ently, so physical effects can arise. It also differs from 
the SME, where distinct particles and fields can break 
Lorentz symmetry in different ways and only a subset of 
coefficients are unobservable. The absence of observable 
effects when the conventional momentum is adopted is 
a known characteristic of DSR models [69], which fol- 
lows from their definition as nonlinear momentum-space 
representations of the usual Lorentz transformations. 

More generally, the above considerations reveal that 
modifications of the dynamics are unphysical whenever 
they arise from a universal and reversible momentum 
substitution. Observable effects might in principle be 
possible in special models if singularities in the physics 
obstruct the recovery of the usual 4-momentum via the 
inverse momentum replacement, although the require- 
ment of physical singularities seems unappealing. How- 
ever, an alternative approach docs exist. Observable ef- 
fects in models with deformed Lorentz transformations 
can be obtained by imposing the deformed invariance 
only on a subset of particles or fields, while either conven- 
tional Lorentz invariance or a different deformed invari- 
ance holds for others. This idea has not been investigated 
in the literature, perhaps because the idea of two or more 
mutually incompatible invariances in nature runs counter 
to the DSR philosophy. Note that any such models are 
subsets of the SME, so constraints from SME coefficients 
apply. 

As one simple exotic example of a model with multiple 
invariances, one can consider 'spinning' special relativity 
in which the deformed symmetry differs for fermions and 
bosons. The quadratic action for a fermion would then 
be invariant under one transformation, while that for a 



boson would have a different invariance. The deformed 
transformations can also be chosen to depend explicitly 
on the representation, so that multiple invariances would 
be involved in the various pieces of the quadratic action. 
Different versions of this idea could be considered. For in- 
stance, 'flavorful' and 'colorful' special relativities could 
be constructed by choosing the deformations to vary with 
the particle species, or more specifically with the repre- 
sentation of the internal symmetry group. In any case, 
the full action generically breaks the individual invari- 
ances, thereby leading to detectable signals. As before, 
any subgroup of the deformed transformations that leaves 
invariant the full action is associated with unobservable 
effects. 



V. REFERENCE FRAMES AND ROTATIONS 

For comparative purposes, it is useful to adopt a stan- 
dard inertial frame in reporting measurements of coeffi- 
cients for Lorentz violation. The canonical frame used 
in the literature is a Sun-centered celestial equatorial 
frame [11, 34]. Cartesian coordinates in this frame are 
denoted (T, X, Y, Z) . The Z axis lies along the rotation 
axis of the Earth, while the X-Y plane coincides with the 
Earth's equatorial plane. The X axis is directed from the 
Earth to the Sun at the vernal equinox. One advantage 
of this conventional choice is that transforming between 
the Sun-centered frame and a laboratory frame is com- 
paratively simple. 

In a typical application, a measurement of coefficients 
for Lorentz violation is made in a laboratory frame of ref- 
erence. However, the rotation and revolution of the Earth 
imply this frame is noninertial, so the coefficients for 
Lorentz violation in the laboratory vary with time. These 
varying coefficients are related to the constant ones in 
the Sun-centered frame by an observer Lorentz transfor- 
mation, which predominantly involves rotations. In this 
section, we construct the relevant rotation transforma- 
tions to an arbitrary laboratory reference frame, includ- 
ing one based on a rotating turntable. We apply the re- 
sults to obtain the transformation of spherical coefficients 
for Lorentz violation between the laboratory and Sun- 
centered frames. Since the spherical coefficients repre- 
sent a decomposition based on angular momentum, their 
transformation is comparatively straightforward. More- 
over, the spin weight is unaffected by rotations because 
helicity commutes with the angular momentum J. The 
only index that changes under rotations is therefore the 
Jz eigenvalue m. 



A. Rotation matrices 

To construct the rotation transformations, we adopt 
Euler angles a, /3, and 7 that relate two arbitrary 
cartesian frames with coordinates {x, y, z) and {x' , y' , z') 
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through the rotation 



(129) 



This rotation can be visuahzed by starting with the two 
frames coinciding, rotating the second frame by 7 about 
the z axis, then by (5 about the y axis, and finally by 
a about the z axis again. Acting on coordinates, this 
combination of rotations can be shown to be implemented 
by the matrix equation 



( cos 7 sin 7 
-sin7 COS7 



V 







1 



cos /? — sin /? 
1 
y sin/3 cos/3 



V 



COS a 


sin a 






— sin a 


COS a 


•1 
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Note that this result implies that the net rotation can 
alternatively be viewed as a rotation by a about ^, fol- 
lowed by a rotation by /3 about the rotated y axis, and 
then by a rotation by 7 about the rotated z axis. 

To determine the rotation rules for spherical harmon- 
ics, consider a spin-weighted function g/ expanded in 
both coordinate systems: 

= E sYimm = E /i™ sY,^{^). (131) 

The spherical harmonics are related by 
which implies the transformation 



(132) 



(133) 



In these expressions, the quantities 



^^L'(",/3,7) 



(134) 

are the Wigner rotation matrices [70]. It can be shown 

that these matrices are independent of s, which is to be 
expected because J commutes with the hclicity operator. 
Consequently, all spherical coefficients arc rotated using 
the same set of matrices, regardless of their spin weight. 
The Wigner matrices are often written in the form 



(135) 



1)^^,(0, /3,0) are called the little 



where d^±,{P) 
Wigner matrices. The phases in this equation corre- 
spond to the two rotations about the z axis, while d^I^ , 



accounts for the rotation about y. Explicitly, the little- 
matrix elements are given by 



jO') ..fc+mW V(j+m)!0-m)!(j+m')!0-m')! 

"mm'V^'/ — y-^> (3-m-k)\{m-m'+k)\{3+m'-k)\k\ 
k 

X (cosf)'^(tan|)''+"', 

(136) 

where the sum is restricted to all k for which the argu- 
ments of the factorials are nonnegative. 



B. Laboratory frame 

We next apply the Wigner matrices to relate the coef- 
ficients for Lorentz violation in the laboratory frame to 
those in the Sun-centered frame. For this purpose, it is 
useful to work with a canonical laboratory frame [11]. 
Cartesian coordinates in this frame are denoted (x, y, z). 
The z direction is directed towards the zenith, and the 
X axis lies at an angle (j) measured east of south. The 
colatitude of the laboratory is denoted x- The orienta- 
tion of the laboratory with respect to the Sun-centered 
coordinates {X, Y, Z) is determined by the local sidereal 
time Tq. Since the Z axis points towards the celestial 
north pole while the X and Y axes lie in the equato- 
rial plane with right ascension 0° and 90°, respectively, 
it follows that the laboratory z axis points toward right 
ascension w^T®, where cj® ~ 27r/(23 hr 56 min) is the 
Earth's sidereal frequency. 

With these conventions, the angles (j), X) and ui^T^ 
represent three Euler angles giving the relevant net rota- 
tion between frames. From the perspective of the Sun- 
centered frame, the laboratory frame is obtained by ro- 
tating by (j) about Z, then by x about Y , and lastly 
by uj^T^ about Z. The Euler angles relating the Sun- 
centered frame to the laboratory frame are therefore 



a = uj®T^, P = X, 1 = ' 
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It follows that the explicit rotation relating the two sets 

of coordinates is 



cos X — sin X 

1 
sin X cos X 

(x\ 




cosa;0T0 smoji^T,^ 
— sina;0T0 coswqTq 
1 



Y 



(138) 



The spherical coefficients in the laboratory frame can now 
be expressed as Sun-frame coefficients through the rela- 
tion 



lab 



m' 

= E^^ 



mm 



Sun 



'''''d'^±'{-x)^^, (139) 
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where IC^m and /C|™ represent arbitrary spherical coef- 
ficients for Lorentz violation in the laboratory and Sun- 
centered frames, respectively. 

The result (139) is auspicious. A key signal in many 

experiments is the sidereal variation introduced by the 
rotation of the Earth. Here, this rotation is expressed in 
a simple form involving time-dependent phases e*'"'^®^® . 
Only the colatitude x appears in the little Wigner ma- 
trices c^TO'rTi(~x)) which arc time independent. Conse- 
quently, for typical applications these time-independent 
factors need be calculated only once for a given experi- 
ment at fixed X- 

Some experiments involve turntables rotating about 
the vertical axis. This situation can be incorporated into 
the above rotation by fixing the laboratory frame with re- 
spect to the turntable. This implies the azimuthal angle 
acquires a time dependence of the form cj) = WttTtt, where 
cott is the turntable rotation frequency and Ttt is mea- 
sured from a time when the x axis points south. Again, 
the time dependence enters through simple phases. 

As a simple example, consider the rotations of vector 
coefiicients in an experiment involving a turntable. Cal- 
culating the Wigner matrices for j = 1, we find that the 
rotation between spherical coefficients in the laboratory 
and Sun-centered frames is given by 








laboratory-frame coordinates by 





/ cos^ I 



sm 



2 X 



^smx cosx 
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sin 



2 I 



-^sinx 



cos^ ^ 
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in matrix form. 
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Sun 
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In some situations, it may be convenient to define a 
third frame that is fixed with respect to the appara- 
tus. The advantage of this third frame is that it may 
be chosen to simplify calculations. For example, a labo- 
ratory apparatus often has a symmetry axis, so adopting 
a third apparatus frame with one coordinate axis along 
the symmetry direction may be convenient. To incorpo- 
rate this in the above formalism, it suffices to identify 
a suitable apparatus frame and to determine the corre- 
sponding Euler angles relating it to the laboratory frame. 
The laboratory-frame and apparatus-frame spherical co- 
efficients are then related through Eq. (133). Thus, if 
the apparatus-frame coordinates {x' , y', z') are related to 



fx' 

y' 



I 



cos 7 sin 7 

— sin 7 cos 7 
1 



cos a, sin a 
— sin a cos a 
1 
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where a, /?, 7 are appropriate Euler angles, then the 
apparatus-frame spherical coefficients can be written 
in terms of the laboratory-frame spherical coefiicients 
through 



j^-^app 



lab 



(142) 



Assuming the orientation of the system is fixed in the 
laboratory frame, the Wigner matrices for this rotation 

are constant. 

As a simple illustration, consider a system with sym- 
metry axis oriented along the x axis of the laboratory 
frame. In this case, it may be beneficial to choose an ap- 
paratus frame having angular-momentum projection axis 
z' aligned with the symmetry axis. This can be achieved 
by taking x' = —z, y' = y, and z' = x. A suitable choice 
of Euler angles is then given by a = 0, /3 = 90°, and 
7 = 0. 



VI. ASTROPHYSICAL TESTS 

In this section, we discuss searches for Lorentz viola- 
tion involving observation of radiation from sources at 
cosmological distances. Due to the large baselines in- 
volved, searches for birefringent and dispersive effects can 
in principle achieve high sensitivities to almost all the 
vacuum coefiicients cj^]^.^, fcj^'^^.^, kfj^j^, and kl^^.^ 
introduced in Sec. IV C. 

Where relevant, vacuum birefringence provides consid- 
erably greater sensitivity than vacuum dispersion, as is 
shown below. It is therefore natural to separate the rel- 
evant Lorentz-violating operators into two classes. The 
first class is controlled by the vacuum coefficients k'f^^j^, 



.(d) 



"'(B) jm' '^iv)jm produces leading-order birefringence. 
The second consists of vacuimi operators causing disper- 
sion without leading-order birefringence, and is associ- 



ated with the coefficients c, 



for d> 4. Observations 



-{I)jm 

of dispersion are therefore well suited to measurements 
of c|^j^.^, while studies of birefringence are appropriate 
for measurements of the remaining coefficients. 

We begin in Sec. VIA with a discussion of dispersion 
tests. The basic theory and results are summarized, and 
new constraints are obtained using recent results. In Sec. 
VI B, we consider birefringence tests. A general treat- 
ment is first outlined, and then applications to point 
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sources and to the cosmic microwave background (CMB) 
are presented. We obtain new constraints from gamma- 
ray bursts on coefficients of mass dimensions five, seven, 
and nine, and we discuss some general features of the 
effects of Lorentz violation on the CMB. 



A. Dispersion tests 

Frequency-dependent photon velocities arise from 
Lorentz- violating operators with d 7^ 4, which cause wave 
dispersion. Astrophysical searches for vacuum dispersion 
seek these differences in the velocity of light at differ- 
ent wavelengths. Typical searches involve explosive or 
pulsed sources of radiation, such as gamma-ray bursts, 
pulsars, or blazars, that produce light over a wide range 
of wavelengths in a short period of time. Assuming the 
temporal structure of the emission is sufficiently well un- 
derstood, observed but unexpected arrival-time differ- 
ences can be interpreted as wavelength dependences in 
the velocity. Even if detailed frequency information is 
unavailable, limits can still be obtained from the pulse 
width because dispersion results in a spreading of wave 
packets. 

A number of astrophysical searches for a modified pho- 
ton dispersion relation have been performed. Most of 
these studies assume isotropic Lorentz violation. To 
make the connection between these approaches and the 
SME, we momentarily restrict our attention to the vac- 
uum isotropic model discussed in Sec. IV C, which is the 
relevant limit for astrophysical studies of isotropic vio- 
lations. Recall that this model has exactly one nonzero 
spherical coefficient for Lorentz violation at each d, con- 
sisting of c|jjgQ for even d and fc^y^Qg for odd d. In the 
isotropic limit, the velocity defect arising from Eq. (75) 
is given by 

-^^^ -^E^'"'(-4')oo±0' (143) 

d 

in terms of the photon energy E. The coefficients c'^^f^^ 
are associated with CPT-even operators producing dis- 
persion but no leading-order birefringence, while nonzero 
coefficients fc^yjoo i^iply birefringence and corresponding 
changes in polarization. Among the studies involving 
modified dispersion relations with isotropic Lorentz vio- 
lation, it follows that only those with odd d, and birefrin- 
gence or even d and no birefringence are consistent with 
linear effective field theory in flat spacetime [71]. 

Isotropic Lorentz-violating effects in modified disper- 
sion relations are sometimes described using a velocity 
deviation of the form 6v = ±^iE, where ^1 is a constant 
[62, 72]. This model is phenomenologically equivalent 

to the single SME coefficient k^^y^Q = ^i- How- 

ever, because this particular combination causes bire- 
fringence, its best constraints currently come from po- 
larimetry observations, which are discussed in the next 



subsection. An isotropic higher-order correction of the 
form Sv = has also been considered [73]. This case 
corresponds to the d = 6 coefficient c|^jQ(, = — \/47r ^2- 
Bounds on this term from the active galaxy Markarian 
501 are currently of order 10~^^ GeV~^ [23], though some 
evidence for nonzero dispersion from this source exists 
[74]. Other cases that arc sometimes considered involve 
an isotropic nonbirefringent linear defect 6v = £,iE [75] 
or an isotropic birefringent defect Sv = ±^2E^, where 
the sign indicates helicity [76]. Both these cases are in- 
consistent with the present general analysis and may be 
problematic. 

In principle, searches for Lorentz violation via disper- 
sion are sensitive to all coefficients with d ^ A. How- 
ever, considerably greater sensitivities to coefficients as- 
sociated with birefringence are typically accessible via 
polarimetry. This can be understood as follows. In a vac- 
uum dispersion study involving a source at baseline dis- 
tance L, the quantity of interest is the change 6t ~ 5vL 
in arrival time of the signal, which implies a sensitiv- 
ity to 5v given by St/ L. In contrast, for a polarimetric 
study of the same astrophysical source, the quantity of 
interest is the phase difference S(p — ESvL of the eigen- 
modes, which yields sensitivity S<j)/LE. Comparing these 
two results, we see that dispersive sensitivity depends on 
the difference in arrival time, while polarimetric sensi- 
tivity depends on the periodicity oc 1/E of the source 
radiation. Consequently, to achieve similar sensitivity, 
an astrophysical dispersion test requires a time resolu- 
tion comparable to the inverse frequency of the photons, 
which is infeasible. 

We can therefore conclude that astrophysical disper- 
sion studies are best suited to searches for the nonbire- 
fringent vacuum coefficients c^jj^-^. Note, however, that 
birefringence also causes a spreading of wave packets, due 
to the differences in velocity of the two birefringent eigen- 
modes. Moreover, all birefringent operators with d 7^ 4 
are also dispersive. A definitive interpretation of an ob- 
served arrival-time difference as a dispersive nonbirefrin- 
gent effect associated with the vacuum coefficients c| jj^-^ 
therefore requires a polarimetric study of the signal or 
elimination of possible contributions from the other vac- 
uum coefficients via independent studies. 

Setting all other coefficients to zero, the velocity defect 
including anisotropics is given by 

5^ ^ _,o ^ _ ^ ^Y.^^f^) c[%^. (144) 

djm 

Note that this involves only even-dimensional Lorentz- 
violating operators and that the d — 4 case involves no 
dispersion. For operators with d > 4, the sensitivities 
increase with energy and so high-frequency sources can 
be expected to yield the sharpest results. Note also that 
limiting attention to isotropic dispersion disregards a to- 
tal of (rf^ — 2d — 2) independent types of vacuum Lorentz 
violation at each d. 

The difference in velocity between photons of different 



35 



Model 



Coefficients 



Result 



System 



Vacuum 



E,„.oy,™(99.r,240°)c(f),„ 

E,>„oF.>.(147M20°)c[?;,^ 
E,^oF,>„(330°,-30°)f;,,^ 

E,,„oyi™(99.7°,240°)cW„ 
E,.n0y.>n(147M20°)f^^ 



< 1 X 10-^** GeV~^ 
3tl X 10-2^ GeV-2 

< 3.2 X 10"^° 

< 7.4 X 10-22 GeV-2 

< 3 X 10-1=^ GeV"-* 

< 2.6 X 10-23 GeV"* 



GRB 021206 
Markarian 501 
GRB 080916C 
PKS 2155-304 

GRB 021206 
GRB 080916C 



< 4 X 10"^® GeV-2 

lOl^ X 10-22 GeV-2 

< 1.1 X 10"^^ GeV"2 

< 2.6 X 10^21 GeV"2 



Vacuum isotropic 



(6) 

^(1)00 
I (6) I 



(8) I 
^(7)00 



< 9 X 10" 
< 9.2 X 10" 



GeV"' 
3 GeV" 



GRB 021206 

Markarian 501 
GRB 080916C 
PKS 2155-304 

GRB 021206 
GRB 080916C 



TABLE XIII: Constraints on spherical coefficients from astrophysical dispersion studies. The first five rows give constraints on 
the vacuum coefficients with d = 6, 8. The next five rows give the constraints on coefficients in the isotropic limit, for which 
there is exactly one nonzero coefficient for each d. Except for the limits from GRB 080916C and PKS 2155-304, which are 
obtained in the text, all results in the table are talcen from the analysis of Ref. [23], which used data from Refs. [74] and [77]. 
The bounds shown are at the 95% confidence level. 



energies leads to an arrival-time difference given by [23] 
i2 - ii « / — 7j — dz 

Jo ilz 

(145) 

where the source redshift is z and ti, t2 are the propa- 
gation times for photons with observed energies Ei, E2 
and velocities vi, V2- Also, 

= Holflril + zf + + zf + + zf + Oa]'/' 

(146) 

is the Hubble expansion rate at z, expressed in terms of 
the present-day Hubble constant i?o — 71 km s-^Mpc""'^, 
radiation density fi^ — 0.015, matter density fim ~ 0.27, 
vacuum density JIa ~ 0.73, and curvature density Vl^ = 
1 — — f7„,, — Qa- In Rc;f. [23] , these expressions are used 
to place direction-dependent bounds on combinations of 
the coefficients c|^j^.^ and c^^i-^j„i^ using observations of 
GRB 021206 [77] and of the blazar Markarian 501 [74]. 
These results are summarized in Table XHI. 

As an illustration, consider the recent measurements 
made by the Fermi Observatory on the source GRB 
0809 16C [78]. This exceptionally energetic source pro- 
duced a burst of photons with observed energies ranging 
to 13. 22^*;;™ GeV, all of which arrived within 16.54 s 
of the initial detection of low-energy photons. The high 
photon energies and the large redshift of 2; = 4.35 ± 0.15 
make this burst a sensitive probe of Lorentz violation. 
A conservative bound on the vacuum coefficients for 



Lorentz violation can be obtained using the la lower lim- 
its for the energy and redshift. Performing the integral 
(145) for d = 6 and assuming the lower energy is negligi- 
ble, we find the constraint 

^o>S™(147°,120°)cg.„ < 3.2 X IQ-^^ GeV'^ (147) 

jm 

on a direction-dependent combination of operators for 
Lorentz violation of mass dimension d = 6. Taking in- 
stead operators of mass dimension d = 8, we obtain the 
constraint 

^o:i"im(147°,120°)cg.„ < 2.6x10-23 GeV-^ (148) 

jm 

Note that these are one-sided bounds that suppose 
the higher-energy photons propagate more slowly and 
that disregard possible burst-timing structure from the 
source. Under the assumption that the observed time dif- 
ference is due to Lorentz- violating effects, a careful study 
of the leading edge of the high-energy photons might also 
permit the derivation of a lower positive bound for each 
of the above coefficient combinations. 

Another example is provided by the recent data ob- 
tained for the active galaxy PKS 2155-304 by the High 
Energy Stereoscopic System (HESS) [79]. This source 
has redshift 2: = 0.116 and a light curve spanning an en- 
ergy range of a few TeV, with time delays of a few tens of 
seconds. The reported analysis places a constraint of 41 
s TeV-2 at the 95% confidence level on dispersion effects 
quadratic in the energy. Performing the integral (145) 
as before and identifying the reported constraint with 
(^2 — ti)/{E2 — El) implies a conservative bound on a 
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direction-dependent combination of vacuum coefRcients 
for Lorentz violation with d = 6. We obtain 



E 



0^m(330° 



-30°) c 



(6) 



< 7.4 X 10-22 GeV" 



(149) 

at the 95% confidence level. A more complete study of 
the existing data from this source could yield additional 
constraints involving operators of mass dimension d > 8. 

The above analyses demonstrate that a single point 
source provides sensitivity to only a limited number 
of direction-dependent combinations of coefficients for 
Lorentz violation. Multiple sources are needed to ac- 
cess all coefficients for a given value of d. For example, 
there are 25 coefficients c^P} . and 49 coefficients c[ ^ ■ i 
and so even when birefringence is neglected we see that 
a corresponding number of sources at different locations 
on the sky is required to constrain fully these coefficients 
and the corresponding types of Lorentz violation. Data 
from gamma-ray bursts and other burst sources obtanied 
by existing telescopes, including Fermi, HESS, the Major 
Atmospheric Gamma-ray Imaging Cherenkov Telescope 
(MAGIC) [74], and the Very Energetic Radiation Imag- 
ing Telescope Array System (VERITAS) [80], or by fu- 
ture telescopes such as the Advanced Gamma-ray Imag- 
ing System (AGIS) [81], the Cherenkov Telescope Array 
(CTA) [82], and the High Altitude Water Cherenkov Ex- 
periment (HAWC) [83], could be combined to measure 
completely the coefficients c'^yrn ^'^^ various fixed values 
of d > 6. 

In contrast, in the limit of rotation invariance we re- 
cover the vacuum isotropic model, which reduces the 
number of coefficients to one for each d. This implies 
that a single source suffices to place constraints when 
only one value of d is considered at a time. For exam- 
ple, in the vacuum isotropic model a single constraint 
< 10-22 GeV-2 in any location on the sky suffices to 
exclude the suggestion of a signal for d = 6 Lorentz vi- 
olation from Markarian 501, whereas the general SME 
treatment requires at least 25 independent sources at 
this constraint level. For the limiting case of the vacuum 
isotropic model, the above bounds from Fermi reduce to 
the one-sided constraints 



and 



c[^joo < 1.1 X 10-19 GeV-2, 



cj^joo < 9.2 X 10-23 GeV" 



while the one from HESS reduces to 



|C(jjoo| < 2.6 X 10-21 GeV" 



(150) 



(151) 



(152) 



Both the isotropic and the anisotropic constraints ob- 
tained from GRB 080916C and PKS 2155-304 are also 
included in Table XIII. 



B. Birefringence tests 

In birefringent scenarios, the two eigenmodes propa- 
gate at slightly different velocities. This implies that the 
superposition of the modes is altered as light propagates 
in free space. Since the two modes differ in polarization, 
the change in superposition causes a change in the net 
polarization of the radiation. This provides a signature 
of Lorentz violation. In the present subsection, we out- 
line the theory of these polarization changes and discuss 
birefringence tests based on polarimetry using both point 
sources and the CMB. 



1. Theory 



The direction of a Stokes vector s — {s^ , , s'^)'^ m 
the abstract Stokes-parameter space uniquely character- 
izes the polarization of the associated plane wave [84]. 
A Stokes vector in the s^-s^ plane corresponds to linear 
polarization, while a Stokes vector along the axis rep- 
resents circular polarization. Other directions represent 
general elliptical polarizations. 

In this picture, birefringence can be understood as a 
rotation of the Stokes vector s = {s^ , s"^ , s^)-'' about the 
birefringence axis <; = (q^ , q'^ , (;^)'^ . The birefringence 
axis represents the polarization of the eigenmodes and is 
determined by the properties of the medium. The total 
angle of rotation is equivalent to the change in the relative 
phase between the two eigenmodes. With the normaliza- 
tion adopted in Eq. (81), we can write this rotation in 
differential form: 



ds/dt = 2E(; x s = -iT, ■ s, 



(153) 



where E is the photon energy and T,^'^ = —2iEe^*'\^ 
represents the rotation generators. 

Some basic features of birefringence in the context of 
Lorentz violation can be extracted from this picture [11]. 
If CPT is conserved, then ^ lies in the s^-s^ plane, and so 
the birefringent eigenmodes are linearly polarized. Lin- 
early polarized radiation therefore typically rotates out of 
the s^-s^ plane and becomes elliptically polarized. Sim- 
ilarly, circularly polarized radiation rotates away from 
the axis, becomes elliptical, and may eventually ro- 
tate through a linear polarization. In contrast, if CPT is 
violated, then the birefringent eigenmodes are circularly 
polarized with one being left-handed and the other right- 
handed. The rotation axis <; is therefore aligned with the 

axis in this case. As a result, linear polarizations re- 
main linear, but a change in the linear-polarization angle 
occurs. However, circular polarizations remain circular 
because they are eigenmodes of the propagation. 

In a typical application, one considers a distant source 
of polarized light and integrates the rotation (153) from 
emission to detection. This yields changes in polarization 
that can depend on both energy and direction of propa- 
gation. To search for birefringence, we can either model 
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the polarization at the source and seek discrepancies in 
the observed polarization, or we can test for unexpected 
energy dependence in the polarization parameters. 

Since we are characterizing Lorentz violation using 
spin-weighted spherical coefficients, it is convenient to 
reformulate the rotation of the Stokes vector in terms of 
spin-weighted Stokes parameters. We can decompose the 
components of the Stokes vector according to their spin 
weight with respect to the line of sight n. This yields 
S(o) = as a Stokes parameter of spin weight and 
S(±2) = T is^ as two Stokes parameters having spin 
weight ±2. 

In the Stokes basis with s = (s(+2)5 S(o)) S(-2))"^) the 
rotation generators are given in matrix form as 



^ = 2E\ -i 



^(0) -<?(+2) 



';(-2) -C(o) 



(154) 



Here, <; = (^(+2)) '?(o)) '^(-2))"^ is the birefringence axis in 
this basis, with components given by 



djm 



{V)jr. 



(155) 



djm 



Since the rotation of the Stokes vector is determined by 
the combination the effects of birefringence enter 
with an energy dependence of E'^~^. This implies that 
an increased sensitivity to coefficients with > 3 can be 
achieved by using higher-energy photons. We also see 
that unconventional energy dependence is a signal for 
Lorentz violation. Only the d = 3 case leads to energy- 
independent birefringence. 

For many astrophysical sources, cosmological expan- 
sion is significant during the time of flight and must be 
incorporated in the analysis. We implement this by ex- 
pressing the differential rotation in terms of redshift: 



ds 



iYiz ■ s 
(1 + z)H, 



(156) 



where T,z represents the rotation matrix at the 
blucshiftcd energy (1 + z)E and source direction n. To 
obtain the net polarization change, we then integrate this 
expression from source redshift z to 0. 

Some searches for Lorentz violation investigate 
Lorentz-violating operators of a specified mass dimen- 
sion d. With this assumption, the calculation of the net 
rotation is simplified because the energy integral is inde- 
pendent of the matrix multiplication. In the CPT-odd 
case with a single odd value of d, the rotations of the 
Stokes vector about the axis lead to a change in the 
linear-polarization angle ip without affecting the degree 
of linear or circular polarization. The rotation is diago- 
nal in the spin-weighted basis, and we obtain the simple 



result 



S(±2) 
S(0) = 5(0)2 



(157) 



relating the present-day polarization to the original po- 
larization at redshift z. Here, the change 6^z in polar- 
ization is given by the integral 



id) 

{V)jm- 



(158) 

The linear-polarization angle tjj at the present epoch is 
then related through 



(159) 



to the blueshifted angle tpz- 

In the CPT-even case with a single even value of d, 
the cigcnmodcs arc linearly polarized and the rotation is 
more complicated. It is convenient in this case to define 
the direction-dependent phase 



»4(n) _ ^(+2)(^) 



l^(+2)(' 



n 



(160) 



which controls the evolution of the polarization. The 

phase angle ^ is twice the polarization angle of the eigen- 
mode of propagation. Consequently, linear polarizations 
with angle tp = ^/^ ov ip = + 90° remain unaffected 
as the radiation propagates. Calculation shows that the 
redshift integral can be expressed using this phase and 
the angle 



^z=E' 



d-3 



f 

Jo 



(1 + ^) 



d-4 



-dz 



{B)jm) 

The net rotation is given by 

s = mz- Sz, 
where the Miiller matrix m, takes the form 



(161) 



(162) 



COS^ $2 

k sin 2$,e'« 



-i sin 2$2e~*^ sin^ ^^e"^*^ 
^sin2$^e-'« 



COS 2$2 
isin2$,e'^ 



COS^ $2 



(163) 



sin^ $^6^'^ 
in the spin- weighted Stokes basis. 

2. Point sources 



We next use the above theoretical results to obtain 
bounds on spherical coefficients from polarimetry of as- 
trophysical point sources. Some of the tightest existing 
constraints on Lorentz violation have been achieved in 
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FIG. 2: Spectropolaximetry bounds from GRB 930131 on 
spherical coefficients corresponding to Lorentz-violating op- 
erators of mass dimensions d = 5,7, 9. The source di- 
rection is given by the angles n = (98.2°, 182.1°). The 
solid curve represents the effective degree of polarization 
Ileff = \/ (s"*")^ + (s^)^. The effective polarization angle 
V'efF = tan~^(s^)/2(s^) in degrees is displayed with crosses. 
The shaded area is the disallowed region with Xlefi < 35%. 
All coefficients are in units of GeV*"''. 



FIG. 3: Spectropolarimetry bounds from GRB 960924 on 
spherical coefficients corresponding to Lorentz-violating op- 
erators of mass dimensions d = 5, 7, 9. The source di- 
rection is given by the angles n = (87.3°, 37.3°). The 
solid curve represents the effective degree of polarization 
Ilcff — \J (s^)^ + (s^}^. The effective polarization angle 
i/icff = tan^^(s^)/2(s"'^) in degrees is displayed with crosses. 
The shaded area is the disallowed region with Ileff < 50%. 
All coefficients are in units of GeV"''. 



birefringence searches involving high-frequency sources 
such as gamma-ray bursts [21, 72, 76, 85]. However, 
point sources have the disadvantage that a single line of 

sight n is involved, which provides sensitivity to only a 
restricted portion of the space of coefficients for Lorentz 
violation. As with astrophysical dispersion tests, multi- 
ple sources are therefore required to perform a compre- 
hensive search, even for a fixed value of d. For example, 
for the case of = 3 a multiple-source search involving a 
large number of radio galaxies [19] has placed a limit on a 

quantity Pq, = —2(k\p)a corresponding to the constraint 



E 



(3) 



< 6 X lO-''^ GeV 



(164) 



at the 95% confidence level in terms of spherical co- 
efficients. In the isotropic limit, this gives the limit 
< 2 X 10^^^ GeV, although sharper bounds have 



,(3) 
''(V)OO 



recently emerged from CMB polarimetry as discussed be- 
low. Multiple-source searches for the case of d = 4 have 
also been performed [11, 20, 21, 23]. For example, a 
search using 16 sources [11] places a limit that translates 
in the present context to the rotationally invariant con- 
straint 



' + |fc(sWP)<5xlO- 



32 



(165) 



at the 95% confidence level. For larger values of 
d, multiple-source birefringence analyses offer excellent 
prospects for systematic tests of Lorentz violation at ex- 
treme sensitivity. 



In this subsection, we illustrate the procedure and ob- 
tain first constraints on some spherical coefficients with 
larger d, by analyzing the recent evidence for significant 
polarized components in the radiation from the gamma- 
ray bursts GRB 930131 and GRB 960924 [86]. Obser- 
vations of gamma rays associated with these two sources 
suggest that they are polarized at levels of Hgaoisi > 35% 
and 11960924 > 50%, respectively. Ideally, if we knew 
the degree of polarization and polarization angles at the 
source, we could search directly for changes due to bire- 
fringence. Without this information, however, we can 
still place limits on decoherence effects caused by bire- 
fringence. The point is that significant birefringence 
would lead to large differences in observed polarizations 
at slightly different frequencies, effectively unpolarizing 
the radiation. Evidence for polarization can therefore be 
used to constrain the frequency-dependent birefringence 
caused by violations with dimension d> 3. 

Figures 2 and 3 show the calculated effective degrees 
of polarization and the changes in the polarization an- 
gles for the above two gamma-ray bursts in scenar- 
ios with nonzero spherical coefficients corresponding to 
CPT-odd Lorentz-violating operators of mass dimensions 
d = 5,7, 9. In constructing these plots, we assume the 
radiation is initially 100% linearly polarized, which is a 
maximally conservative assumption in the present con- 
text. The displayed results are then obtained by numer- 
ically calc ulating the ch ange in the effective polarization 
Ileff = i/ (s^)^ + (s^)^ smeared over observed frequen- 
cies, following the basic procedure outlined in Ref. [21]. 
The shaded regions in Figs. 2 and 3 show the ranges of 
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Model 



Coefficients 



GRB 930131 



GRB 960924 



< 7 X 10"^^ GeV"^ 

< 2 X 10"^^ GeV"^ 

< 6 X IQ-" GeV"^ 



< 4 X 10"^^ GeV"^ 

< 5 X lO"^'^ GeV"^ 

< 1 X IQ-" GeV"^ 



Vacuum 



1 Ej-rn ^(V)jv 
I J2jm oYimifi) k^y^.^ 



< 2 X 10-32 GeV-1 

< 7 X 10"^^ GoV" 

< 2 X 10"^" GoV" 



< 1 X 10-32 GoV^i 

< 2 X 10-2-* GeV-3 

< 4 X lO"^'^ GeV-^ 



Vacuum isotropic 



1.(7) I 

J9) j 
^(v)ool 



-37 



-37 



Vacuum 



(4) 

{E)jm 



y. 



(n) {k 

(8) 



(4) \ 
(B)jm) 



< 10 



< 10 



(n) (fcf ) 



< 10-29 GeV" 

< 10-20 GeV- 



< 10-29 GeV" 

< 10-20 GeV- 



TABLE XIV: Constraints on spherical coefficients from polarization observations of the gamma-ray bursts GRB 930131 and 
GRB 960924. The first three rows give constraints on the vacuum coefficients for the CPT-odd cases with d = 5, 7, 9. The 
arguments of the spherical harmonics are n = (98.2°, 182.1°) for GRB 930131 and n = (87.3°, 37.3°) for GRB 960924. The 
next three rows give the constraints on coefficients in the isotropic limit, for which there is exactly one nonzero coefficient for 
each d. The final three rows give the approximate sensitivities achieved for vacuum coefficients in the CPT-even cases with 
d = 4, 6, 8. Unlike the constraints in the first six rows, the results in the final three rows cannot be interpreted as definitive 
bounds because the amount of birefringence in the CPT-even case depends on details of the source polarization. All constraints 
are at the 95% confidence level. 




FIG. 4: Spcctropolarimctry bounds from GRB 930131 on 
spherical coefficients corresponding to Lorentz-violating op- 
erators of mass dimensions d = 4, 6, 8. The source direction 
is given by the angles n — (98.2°, 182.1°). The shaded area 
is the disallowed region with Ileff < 35%. All coefficients are 
in units of GeV*-'*. 



FIG. 5: Spcctropolarimctry bounds from GRB 960924 on 
spherical coefficients corresponding to Lorentz-violating op- 
erators of mass dimensions d = 4, 6, 8. The source direction 
is given by the angles n = (87.3°, 37.3°). The shaded area is 
the disallowed region with Ileff < 50%. All coefficients are in 
units of GeV*-''. 



coefficient space that are excluded by the observation of 
polarization in the radiation from these sources. Coeffi- 
cients lying in these regions would cause depolarization 
beyond what is observed. The resulting constraints for 
the vacuum coefficients and for the limiting case of the 
vacuum isotropic model are summarized in Table XIV. 

In the CPT-even case, the linear polarization of the 
source could in principle coincide with one of the eigen- 
modes of propagation. This situation occurs when the 
phase angle ^(n) is twice the initial polarization angle 
'ipo, as discussed above. Consequently, the results for any 



given point source contain unbounded regions of coeffi- 
cient space, and so definitive constraints on the spheri- 
cal coefficients cannot be obtained. However, the analy- 
sis does achieve high sensitivities to Lorentz violation in 
certain regions of coefficient space. Figures 4 and 5 show 
the portions of coefficient space excluded by the obser- 
vations of GRB 930131 and GRB 960924 for CPT-even 
operators of mass dimensions d = 4, 6, 8 [87]. The shaded 
areas in these figures represent disallowed regions. Their 
shape and extent demonstrates that part of the coefficient 
space is excluded at high sensitivity, while emphasizing 
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the need for simultaneous analysis of multiple sources to 
obtain definitive constraints. The approximate sensitiv- 
ities to the vacuum coefficients achieved from the two 
sources are summarized in Table XIV. 

The above examples show that high-frequency po- 
larimetry of gamma-ray bursts has the ability to probe 
Lorentz violation at extreme sensitivity. The use of mul- 
tiple sources offers the potential to constrain all of the 
vacuum coefficients . , fc^ll . , and fclfl . associ- 
ated with leading-order birefringent Lorentz violation at 
large d. In the CPT-odd case, the limits can completely 
constrain the single vacuum isotropic coefficient at each 
odd d because the direction of propagation is irrelevant. 

The power of this type of analysis is apparent in com- 
paring the sensitivities in Table XIV to those achieved 
via the astrophysical dispersion tests discussed in Sec. 
VIA. For d = 6, the polarimetry of gamma-ray bursts 
attains a sensitivity roughly a million times beyond that 
of the dispersion tests involving Markarian 501, despite 
the million-fold difference in energy. This confirms the 
suitability of dispersion tests to constrain instead only 
the vacuum coefficients c|^j^.^, which are associated with 
Lorentz-violating operators that have no leading-order 
birefringence. 

3. CMB tests 

In principle, an extended source can provide access 
to all birefringent Lorentz-violating operators and hence 
bypass the major limitation of point sources described 
above. A prime example of an extended source is the 
CMB. Since the CMB is the oldest observable untainted 
radiation and hence represents the longest available base- 
line, an analysis of CMB polarization might be expected 
to yield high sensitivities to Lorentz violation. However, 
this expectation may fail for operators at larger d due 
to the comparatively low CMB frequencies. Nonethe- 
less, the CMB docs provide interesting opportunities for 
lower-dimensional violations, and in particular it is the 
best available source for studies of coefficients with d = 3. 
In this subsection, we discuss and illustrate some of the 
unusual features that can arise in the CMB in the event 
of significant birefringence [22, 23, 88]. 

The accepted description of the CMB uses a decompo- 
sition of the temperature T and the Stokes parameters 
into spin- weighted spherical harmonics [59] , analogous to 
the discussion in Sec. IV C. Typically, power spectra are 
introduced to characterize the strength of each mode and 
the correlations between them, according to 

Cf'""' = ^ E«4™ (166) 

where Xi, X2 range over T, B, V and where are 
the coefficients in the spherical-harmonic expansion of 

X. While temperature anisotropics have been firmly es- 
tablished, the detection of polarization in the CMB by 



several experiments [89-92] is of more significance in the 
present context. 

In the conventional Lorentz-invariant picture, temper- 
ature and density fluctuations at recombination provide 
the necessary anisotropies to produce a net polarization 
[93]. In addition, these processes lead to a correlation 
between temperature and the iS-parity component of the 
CMB. The E component makes up only a tiny fraction 
of approximately 10~^ of the total radiation. Several 
measurements of this small degree of polarization have 
been made. The B-type polarization is expected to be 
even smaller and uncorrclatcd with temperature. The 
standard picture predicts no significant V polarization 
because Thomson scattering produces only linear polar- 
ization. 

The presence of Lorentz violation may alter many of 
the above properties. It can introduce unexpected types 
of polarization, and it can induce mixing between initially 
uncorrclatcd modes during the nearly 14 billion years of 
propagation. The associated violations of rotational sym- 
metry can also cause mixing across multipoles in j and 
TO, including for modes of the same polarization type. 
Here, wc provide a discussion of general features based 
on a numerical survey of some Lorentz-violating models 
exhibiting these effects. The calculations parallel those 
presented in Refs. [22, 23]. 

Some qualitative features can be determined directly 
using the intuition provided by the Stokes rotations de- 
scribed in Sec. VI Bl. For example, the vacuum coef- 
ficients k'fy^j^ are associated with CPT-violating oper- 
ators and produce local rotations of the Stokes vector 
about the axis. This causes a global mixing of the 
linearly polarized E and B modes. The result can be un- 
conventionally large B polarization, although no circular 
V modes can appear. In contrast, the vacuum coefficients 
^{E)jm ^^"^ ^{B)jm associated with CPT-even operators 
cause both mixing between E and B modes and also 
the emergence of V modes, since in this case the Stokes 
vector rotates out of the s^-s^ plane. Consequently, mix- 
ing between the three types of polarization are possible 
in this scenario, although details of the mixing depend 
strongly on the specifics of the Lorentz-violating opera- 
tors involved. 

Other key features that may be present for some 
types of violations include dependences on the photon 
frequency and birefringence varying with the direction 
of propagation. Only the d = 3 vacuum coeflBcients 
lead to frequency-independent rotations. Also, only the 
j = coefficients generate direction-independent rota- 
tions. Consequently, the isotropic vacuum coefficient 
(3) 

^(V)oo Pi'ovides a particularly simple special case. It 
causes a frequency-independent mixing that is uniform 
across the sky and that leads to a simple rotation be- 
tween E and B modes. Using Eq. (158), we estimate 
this rotation for CMB radiation to be 

(5^- ~ lO''^ degree/GeV. (167) 
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Coefficients 


Result 


Reference 




(15 ± 6) X 10"*^ GeV 


[22] 




(lot*) X 10-*=* GeV 


[23] 


1.(3) 


(6.0 ± 4.0) X 10-*=* GeV 


[94] 




(2.5 ± 3.0) X 10-*=* GeV 


[95] 




(12 ± 7) X 10"*=* GeV 


[22] 




(1.2 ±2.2) X 10-*=* GeV 


[96] 




(2.6 ± 1.9) X 10-*=* GeV 


[97] 




< 2.5 X 10-*=* GeV 


[98] 




(2.3 ± 5.4) X 10-*=* GeV 


[23] 




(-1.4 ±0.9 ±0.5) X 10-*=* GeV 


[99] 



TABLE XV: Constraints on spherical coefficients with d = 3 
from CMB studies. The table lists some existing l-tj results 
for the scalar magnitude jfc^^j defined in Eq. (168) and for 
the isotropic component k^y)oo^ ^ obtained via CMB anal- 
yses. Constraints on these quantities from other sources are 
compiled in the data tables of Ref. [6]. 



Simple frequency-independent rotations of this type have 
been considered by several groups [22, 23, 94-99], and the 
existing measurements are listed as part of Table XV. 
Figure 6 illustrates the type of mixing that results. We 
see that initial E power partially rotates into B power. 
Also, the initial TE correlation induces a TB correlation. 
Furthermore, since the B polarization is generated from 
the original E polarization, these two modes become cor- 
related and so a significant EB component emerges. 

Other isotropic rotations from CPT-odd operators of 
larger d, such as the one associated with the isotropic 
vacuum coefficient ^[yjQQ, lead to similar effects. How- 
ever, the frequency dependence introduced by coefficients 
at larger d implies that the amount of rotation depends 
on the photon frequency. Spectral signatures of this type 
should be accessible to observations having sensitivity to 
a wide range of frequencies. 

Any CPT-odd Lorentz- violating operator with nonzero 
j produces direction-dependent rotations. An example is 
shown in Fig. 7, where a comparatively large value of the 
anisotropic vacuum coefficient ^[y^j^g been chosen to 
illustrate the effects. While the local polarization rota- 
tions arc similar to those from the presence of a nonzero 
^(V)oo' anisotropics in this case cause the correlations 
to disappear globally. The effects saturate for large rota- 
tions, and the net result is roughly equal amounts of E 
and B. The overall degree of polarization is unaltered, 
and there is little correlation between any two modes. 

Birefringent operators that are CPT-even are necessar- 
ily both frequency dependent and anisotropic, resulting 
in similar behavior to the previous case. However, a dis- 
tinctive feature of CPT-even violations is the mixing of 
linear and circular polarizations. Figure 8 illustrates the 
generation of V polarization in the presence of a nonzero 

vacuum coefficient ^[tjloo- ^^^^ case, we find that the 
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FIG. 6: Correlation spectra for isotropic CPT-odd Lorentz 
violation. Circles show the effects of Lorentz violation due to 
the spherical coefficient ^(v^joo = 12 x 10-*^ GeV, while lines 
show the Lorentz-invariant case. Power is transferred from E 
to B, and an EB correlation is generated. Also, the initial 
TE correlation induces a TB correlation. As expected for 
CPT-odd Lorentz violation, no significant V polarization or 
correlation appears. The coefficients Cj have units of ^V? . 



anisotropic mixing causes a depletion of E polarization, 
which is rotated into roughly comparable amounts of B 
and V polarization. As in the case of nonzero fc^y^j^Q and 
indeed in all cases we have studied involving nonzero coef- 
ficients with j ^ 0, the anisotropic effects tend to deplete 
correlations when Lorentz violations are large. 

Our survey reveals that similar features as those illus- 
trated above also recur for other vacuum coefficients for 
Lorentz violation that control birefringent operators. We 
thereby find that generic signals of Lorentz violation in 
the CMB can incorporate one or more of the following 
basic features: (a) a depletion in the EE and TE spectra; 
(b) the introduction of unconventionally large B polar- 
ization; (c) the appearance of TB or EB correlations; 

(d) the development of significant V polarization; and 

(e) frequency dependences of the power spectra. 

None of these basic features is readily apparent in the 
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observed data, which suggests that any Lorentz viola- 
tion in the CMB must be small. A study of the high- 
frequency BOOMERANG data [91] places constraints 
on several of the vacuum coefficients with d < 6, find- 
ing some evidence for nonzero Lorentz violation [22]. A 
more recent analysis of the d = 3 coefficients using the 
five-year data from the Wilkinson Microwave Anisotropy 
Probe [90] finds no evidence for isotropic violations in- 
volving the coefficient fcfy\gn, but uncovers some support 



''(y)oo 

for nonzero anisotropic vacuum coefficients k 

J3) 



(3) 

(V)lm 



[23]. 



This study constrains the coefficients k^y-^j^ at the level 
of 10^*'^ GeV. For purposes of reporting constraints, the 
spherical coefficients with d — i can conveniently be sep- 
arated into the isotropic component fc^y-jQo ^md the scalar 
magnitude 



FIG. 7: Correlation spectra for anisotropic CPT-odd Lorentz 
violation. Circles show the effects of Lorentz violation due to 
a large value of the spherical coefficient fc^^'j-^g = 64 x 10"'*^ 
GeV, while lines show the Lorentz-invariant case. Power is 
transferred from E to B, and a loss of TE correlation is evi- 
dent. As expected for CPT-odd Lorentz violation, no signif- 
icant V polarization or correlation appears. The coefficients 
Cj have units of /iK^. 
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FIG. 8; Correlation spectra for anisotropic CPT-even Lorentz 
violation at photon frequency uj = 380 GHz. Circles show 
the effects of Lorentz violation due to a large value of the 



spherical coefficient fcj^''j2p = 16 x 10 



while lines show the 
Lorentz-invariant case. Power is transferred from E to B and 
V, while no significant correlations appear. The coefficients 
Cj have units of /iK^. 



Table XV lists existing constraints on these quantities 
obtained from studies of the CMB. Other limits are given 
in the data tables of Ref. [6] . Improved sensitivities can 
be expected from future CMB data to be obtained by 
various experiments and missions including, for example, 
the Planck satellite [100], the Q/U Imaging Experiment 
(QUIET) [101], the CMBPol mission [102], the E and 
B Experiment (EBEX) [103], the Experimental Probe 
of Inflationary Cosmology (EPIC) [104], and the Spider 
balloon observatory [105]. 



VII. CAVITY EXPERIMENTS 

Laboratory experiments provide alternative methods 
to search for Lorentz violation in electrodynamics and 
have the ability to probe many coefficients inaccessible in 
astrophysical searches. The most common Earth-based 
tests are contemporary versions of the classic Michelson- 
Morley [1] and Kennedy-Thorndike [2] experiments and 
are based on electromagnetic resonant cavities [12-17] . In 
this section, we discuss possible laboratory studies of the 
effects of Lorentz-violating operators of arbitrary mass 
dimension d. We outline a theoretical approach for de- 
termining the resonance frequency of a cavity in the pres- 
ence of Lorentz violation. The approach is illustrated by 
applying it in the context of the camouflage model, for 
which the Lorentz violation has no leading-order birefrin- 
gent or dispersive effects in astrophysical photon propa- 
gation. As a specific example of the techniques, we de- 
rive an explicit result for the fractional frequency shift 
of the TMqio mode in a cylindrical cavity with circular 
cross section, including the time dependence induced in 
the signal by the rotation of the Earth. While we focus 
specifically on resonant cavities, many of the ideas dis- 
cussed here are generic and can be applied in the context 
of other laboratory-based experiments, including space- 
based missions such as the Atomic Clock Ensemble in 
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Space (ACES) [106] or the Space and Time Anisotropy 
Tests (STAT) combining the former Space Time Asym- 
metry Research (STAR) and OPTIS missions [107]. 

A. General theory 

The strategy behind many laboratory tests of Lorentz 
invariance is a search for minute variations in some ob- 
servable with changes in the orientation or velocity of the 
apparatus. For cavity experiments, a suitable observable 
is typically the fractional frequency shift 8v jv induced by 
the Lorentz violation [11]. At leading order in coefficients 
for Lorentz violation, the fractional frequency shift takes 
the generic form bvjv = '^j„^M^jm^jrm where M^jm is 
an experiment-dependent matrix factor and ICjm repre- 
sents the relevant spherical coefficients for Lorentz viola- 
tion discussed in Sec. III. 

Normally the experiment-dependent matrix factors are 
determined in a cavity frame that is fixed with respect 
to the system in question. In this frame, the matrices 
are denoted M^'^ and are constant. However, since this 
frame is noninertial, the cavity-frame coefhcicnts JCJ^ 
for Lorentz violation vary with changes in the orientation 
and velocity of the apparatus with respect to the stan- 
dard Sun-centered frame. In practice, boost effects are 
suppressed by one or more powers of the typically small 
velocities ~ 10^^ involved, so for simplicity in what fol- 
lows we neglect boost effects and focus on violations of 
rotation invariance. 

The rotations relating spherical coefficients in the cav- 
ity, laboratory, and Sun-centered frames are given in Eqs. 
(139) and (142). Since the rotation between the cavity 
frame and the laboratory frame is often constant, it is 
convenient to define laboratory-frame matrices M}-"^ via 

= E^.w^^L(-7,-/3,-a)- (169) 

m' 

In terms of these and the spherical coefficients in the Sun- 
centered frame, the fractional frequency shift takes the 
form 

jram' 

which explicitly reveals the sidereal dependence. For ex- 
periments involving turntables, where the cavity and lab- 
oratory frames rotate relative to each other, the varia- 
tions resulting from the turntable rotation are incorpo- 
rated through the phase (j) = ^ttTtt- 

Manipulation of the modified Maxwell equations (16) 
leads to a perturbative estimate for the fractional fre- 
quency shift given by [11] 

V " -4M / "^'^^^ SD-B*. 6H) (171) 



where 

{U) = \j d^x [E* D + B* H) (172) 

is the unperturbed energy inside the resonator. This for- 
mulation allows for a general linear and lossless medium 
inside the cavity in addition to modifications due to 
Lorentz violation. The fields E, B, D, and H are un- 
derstood to be solutions of conventional electrodynamics 
in the absence of Lorentz violation, while 

6D = KDE ■ E + KDB ■ B + 2kAF X A, 

SH = khe ■ E + khb ■ B - 2{kAF)oA + 2kAFAo 

(173) 

represent leading-order perturbations to the D and H 
fields due to Lorentz violation. The result (171) assumes 
that the fields vanish outside the cavity volume V. For 
simplicity in what follows, we suppose that the resonant 
modes under consideration are nondegenerate. In the 
case of degenerate resonances, Eq. (171) yields a weighted 
average value for the fractional frequency shift. 

To express the fractional frequency shift in terms of 
the spherical coefficients introduced in Sec. Ill, we next 
convert to momentum space. However, some care is re- 
quired to avoid divergences arising from discontinuities 
at the boundary of the cavity. This technical issue stems 
from the differential nature of the operators kp and kAP- 
Although the fields are taken to be continuous through- 
out the interior volume V of the cavity, they may be dis- 
continuous across the surrounding surface S. As a result, 
the fields may not be strictly differentiable over all space, 
and the derivatives in Eq. (173) may be associated with 
^-function behavior on the surface S. The usual recipro- 
cal nature of position and momentum spaces then leads 
to momentum-space representations that fail to vanish 
sufficiently rapidly at infinite momentum, which intro- 
duces divergences in the integral. This issue is absent 
in the minimal SME with d = 3 and d = 4 because no 
derivatives appear in the constitutive relations for that 
case, but it is endemic for Lorentz- violating operators 
with d > 5 and requires a procedure to eliminate the 
divergences. 

One way to address this technical issue is to define new 
everywhere-differentiable fields E and B_ that are equal 
to E and B inside V but that may be nonzero in the 
region outside V where the original fields E and B van- 
ish. We then have two sets of fields, both satisfying the 
Maxwell equations inside V. Derivatives of the extended 
fields E_ and B remain finite everywhere, including on the 
surface 5* of the cavity, but E_ and B need not satisfy the 
Maxwell equations outside the cavity. Using both sets of 
fields, we can construct a finite version of the fractional 
frequency shift (171) by replacing the fields in Eq. (173) 
with their extended versions. This procedure removes 
the divergences at the cavity boundaries. 

Performing a Fourier transform, we obtain a conver- 
gent momentum-space expression for the fractional fre- 
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quency shift (171) given by 

— = -7777V / d?p{E* ■ KDE ■ E- B* ■ KHB ■ B 

V 4{U) J 

+ E* ■ kdb B-B* khe ■ E 
+ ^{ujkAF-p{kAF)o)-{E* xE)). (174) 

This rcsuh is independent of gauge choice, as expected. 
Since the fields and their extensions agree inside V, the 
cavity energy (U) can be calculated using either of the 
forms 

{U} = \j d^P {E* -D + B*- H) 

= \j d^PiE* -D + B* -H). (175) 

Note that both sets of fields are needed in the integral 
(174) for this procedure to work. The extended fields 
control divergences in momentum space, while the unex- 
tended fields restrict the integration to the volume V of 
the cavity in position space. 

The general procedure for determining the effects of 
Lorentz violation on the resonant frequency of a given 
cavity then involves the following steps. First, obtain 
the fields E, B, D, H in the cavity in the context of 
conventional electrodynamics, incorporating in the usual 
way any permittivity and permeability of the media in- 
volved. Next, construct extensions of the fields that 
are everywhere smooth beyond the cavity volume. Per- 
form the Fourier transform to derive the correspond- 
ing momentum-space fields. Then, rising the spherical- 
harmonic expansions of kAF and the k matrices described 
in Sec. Ill, calculate the cavity-frame matrices M^^J^ from 
Eq. (174). Finally, combine the results using Eq. (170) 
to extract the desired fractional frequency shift 5v/v dis- 
playing the explicit sidereal time dependences. 

B. Example: camouflage model 

A general analysis of cavity experiments incorporating 
all spherical coefficients for Lorentz violation is challeng- 
ing due to the large variety of effects. Moreover, although 
terrestrial tests can provide valuable independent and 
fully controlled checks on vacuum birefringence and dis- 
persion, the high sensitivities attainable in astrophysical 
searches for birefringence and dispersion make it reason- 
able to neglect effects from vacuum coefficients in the con- 
text of resonator experiments. Laboratory searches using 
cavities are therefore well suited to study the vacuum- 
orthogonal coefficients for Lorentz violation, which play 
no role in the vacuum propagation of light. 

In this subsection, we illustrate the treatment of the 
spherical coefficients in the derivation of the cavity- 
frame matrices Mj^. For simplicity, we focus primarily 
on Lorentz-violating operators that produce no leading- 
order birefringence or vacuum dispersion. For d = 4, this 



class of operators includes ones corresponding to the co- 
efficients c^(j)jm t^^t have already been widely studied in 
cavity experiments [12-17]. For general d, it corresponds 
to the special subset of the coefficients (c^')^°^^ formed 

by the camouflage coefficients (S^p^)njWn introduced in 
Sec. IV D. 

We begin by considering a scenario involving the coef- 
ficients (cp'')^*^^) controlling nonbirefringent effects. We 
seek an expression for the fractional frequency shift of 
the form 

^ _ V TU^'') (JdhiOE) , . 

y ~ 2^ ■^^'■(cF)njmy'^F )njm- \'-"^) 
dnjm 

To find the matrices . , note that the fractional 

frequency shift in the nonbirefringent case may be written 
as 

where $f is given by Eq. (60) and the derivatives act in 

momentum space, = d/dpi,. Separating explicitly the 
temporal and spatial components gives 

+ {B* xE-E* xBfVa^ 

-25«^B*-B)V„V6]$F. (178) 

Inserting the spherical-harmonic expansion (60) of $f 
and using the identities in Appendix A 2, we can deter- 
mine the contribution to the fractional frequency shift 
from each coefficient {c'p')^n^ ^^d hence identify the ma- 
trices m'^^K . . 

{CFjnjm 

To obtain an explicit integral expression for the ma- 
trices M.^'^'K ■ , we can eliminate the magnetic fields 
using the momentum-space Faraday law uiB = p x E. 
This determines the matrices Ad'}'^^ ^ ■ as integrals over 

\cp jnjm ^ 

the fields E and E_. In terms of the coordinate basis de- 
scribed in Appendix A 2, there are six field combinations 
contributing to the integrals. They can be specified as 

s° = {E+)*E^ + {E-yE_, 

= e:k^ + {E±rE,, 

1(0) = {ErTE,, (179) 

where 

E± = ^{Ee±iE^). (180) 



45 



The combinations s'^ and S(±2) represent smoothed ver- 
sions of the usual Stokes parameters s° and S(±2)- 
The combinations i(±i) and S(o) helicity-(±l) and 
hehcity-0 combinations, which provide two new trans- 



verse Stokes parameters vanishing when Er = = 0. 
In terms of these quantities, the expUcit integral expres- 
sion for the matrices A^l? is 



M 



id) 



d^pp 



n-2 



|(a;^ -/)-^|4±||(+2lSm(p)S(_,_2)(p) + -2ijm(p)S(-2) 



- ((n - 1)^2 + (rf _ 2 - n)p')./^^{+,Yj^{p)l^^,^{p) - -iY,^{p) l^^,){p)) 
+ {{n - - p2) _ (d _ 2 - n)(d - 3 - 3n)/ _ n{n - 1^) oYjm{p) s%p) 

+ {n{n-l)iJ^-{d-2-n){d-3- n)p'') qF,™ (p) |(o) (p) 



(181) 



where uj is the usual resonant angular frequency. 

The result (181) applies both to empty cavities and to 
ones containing a material medium, provided the medium 
is lossless. The form of the integral reveals that the sensi- 
tivity to Lorentz violation depends both on the shape of 
the cavity and on the properties of any material medium 
it contains. This means that different geometries or me- 
dia can be adopted to access different combinations of 
coefficients for Lorentz violation. For example, parity- 
breaking resonators may be used to access parity-odd 
Lorentz violations that cannot be detected at unsup- 
pressed levels using parity-symmetric systems [108]. The 
properties of the resonator enter implicitly through the 
conventional fields E, which arc determined by solving 
the conventional Maxwell equations inside the cavity in 
the presence of the medium, if any. We emphasize that 
the frequency lj is fixed for a given mode, but any mo- 
mentum p can contribute to the integral and hence to 
the matrices M 



id) 



(cp)njm' 

Next, we focus attention specifically on the combina- 
tions of the coefficients (c^'')|y^^ that govern nondisper- 

'^^^^nj^rl coefficients corre- 



sive effects. For d = A, the [c-p 

spond directly to the vacuum coefficients c|^j^.^, as de- 
can therefore 



scribed in Sec. IV A. The matrices M 



(4) 

(I)jm 
<(4) 



be expressed in terms of the matrices M ) ■ . We find 



M 



(4) _ 



= -M 



(4) 



+ \m 



(4) 



(1)00 ~ 4"' '(cf)OOO ^ 4"' '(cf)OOO' 



-2M 



(4) 

(cF)llm' 



■^(ct)22m- 



(182) 



With these results and Eq. (181), the sensitivity of any 
given cavity to camouflage operators for Lorentz violation 
can be determined. 

Typically, the above analysis is performed in the cav- 
ity frame, although many of the equations hold for an 
arbitrary inertial frame. Given the results for the cavity- 
frame matrices Ai , the corresponding expressions for the 
laboratory-frame matrices follow from Eq. (169) and from 
the orientation of the cavity in the laboratory. Applying 
Eq. (170), we finally obtain the fractional frequency shift 



„(4) 



jmm' 



(184) 



which exhibits the variation with sidereal time. 



C. Example: circular-cylindrical cavity 

As an explicit illustration, we calculate in this sub- 
section a number of elements of the matrix A4 for the 
TMqio mode of a cylindrical cavity with circular cross 
section. The cavity is centered at the origin of a cavity 
frame, with the z' direction aligned with the symmetry 
axis. For definiteness, let R be the cavity radius and 2R 
be its length, so that its ends lie at z' — ±R. Here, we 
consider the case of a vacuum cavity for simplicity. The 
solutions for the TMqio mode in the absence of Lorentz 
violation are given by 



For d > 4, the relevant matrices are associated with the 
camouflage coefficients (c?f?')i">2 instead. Using the re- 
lation (90), we find 



M 



(dF)njm 



M 



id) 

(cF)njm 



M 



id) 

(cj?)(ra-|-2)jm' 



(183) 



Jq{p'xqi/ R)p' inside V , 







outside y, 



(185) 



where xqi is the first zero of the Bessel function Jq, and 
where p' represents the radial unit vector. 
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The electric field E is discontinuous at the ends of the 
cavity, and its derivatives are discontinuous at the sides 
p' — R. We therefore must seek a differentiable extension 
E of the electric field of the type described in Sec. VII A. 
This is relatively straightforward to find in the present 
case because we have an analytic solution that could be 
extended to infinity. However, for the purposes of nu- 
merical calculation it is beneficial to construct instead a 
field E_ that vanishes outside a larger volume V' contain- 
ing the cavity volume V. With a field of this type, the 
Fourier transforms can be accurately determined by nu- 
merical integration over only a finite region V of space. 

An extension suitable for numerical work can be con- 
structed with the aid of the C smoothing function 
gq{x;a,b), defined by gg{x;a,b) = for a; < a < 6 or 
6 < a < a;, gq{x; a,b) = 1 for a < b < x or x < b < a, and 

^9(^5 a, b) - (^_^)29+i (g + l + n)!(9 + n)! 

(186) 

otherwise. This function continuously interpolates from 
at a; = a to 1 at X = 6, is constant outside that inter- 
val, and is g-times differentiable everywhere. Note that 
C°° functions of this type also exist, but for numerical 
calculations the above polynomial is easier to handle and 
suffices for our purposes. 

Using this smoothing function, we can define an ex- 
tended electric field by 

E = gg{p'; 2R, R)gq{z'; 2R, R)gg{z'; -2R, -R) 

X Jo{p'xoi/R)p'. (187) 

Taking the integration volume V' as a cube of side length 
4i?, we see that the extended field E_ matches E inside 
cavity volume V, vanishes outside V, and is g-times dif- 
ferentiable. The value of q that is required to ensure the 
finiteness of the integrals depends on the mass dimen- 
sion d of the Lorentz-violating operators being consid- 
ered. For a given dimension d, the matrices kde, i^HB, 
and koB involve d — 4 derivatives. Use of the Faraday 
law to eliminate the magnetic field B introduces another 
derivative. Consequently, the extended field E_ must be 
at least C"*""^, so choosing g > d — 4 should be sufficient. 

The next step involves obtaining the Fourier trans- 
forms of the field components E^^i, Eyi, E^', E_^i, E^,, 
E_^, . We implement these via fast Fourier transform over 
an X X grid in the volume V' . We use the trans- 
formed fields to determine the momentum-space Stokes 
parameters from Eq. (179), and we then perform numer- 
ical integration to evaluate the integral (181) and obtain 
values for M'^^J^-^^.^. Lastly, we use Eqs. (182) and (183) 
to determine the components of interest for the matrices 
A^^rv^^ and M^'^'^^^ . As expected, the results con- 

verge to stable g-independent values for large A, pro- 
vided g is sufficiently large. Table XVI summarizes the 
results obtained through this procedure for d = 4, 6, 8. 
Note that the symmetry of the mode implies no contri- 
bution from parity-breaking coefficients with odd j or 
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TABLE XVI: Nonzero matrix elements Ml}'^^ and 
^(d)cav ^j^^ fundamental mode of a circular cylindrical 

( C F)njm 

cavity. Values for the cases d = 4, 6, 8 are displayed. The 
units are R*~'^, where R is the radius and half the length of 
the cavity. 

from coefficients with m ^ 0, a result that is recovered 
numerically. 

Cavity experiments searching for Lorentz violation 
typically compare two identical cavities with different ori- 
entations. In the present example, the cavity-frame ma- 
trices A4 are then identical for the two cavities. However, 
their differing orientation implies their laboratory-frame 
values would differ. This leads to the slight frequency 
difference that constitutes the signal for Lorentz viola- 
tion. 



VIII. SUMMARY AND DISCUSSION 

In this paper, we derive and study gauge-invariant 
Lorentz- and CPT-violating terms associated with the 
effective photon propagator in the Lagrange density of 
the SME, allowing for operators of arbitrary mass di- 
mension d. We begin by showing that Lorentz viola- 
tion at mass dimension d is characterized by a set of 
{d + l){d — l){d — 2)/2 independent coefficients for the 
CPT-odd case and another set of {d + l)d{d — 3) inde- 
pendent coefficients for the CPT-even case. The compact 
Lagrange density (8) incorporates these effects for all d. 
It includes and extends the pure-photon sector of the 
minimal SME [7, 11]. 

The Lagrange density (8) implies the equations of mo- 
tion (11) for the photon field. An interpretation of these 
equations is elaborated in terms of electrodynamics in 
macroscopic media via the introduction of the operator 
constitutive tensors (14). We derive the covariant scalar 
dispersion relation (30), which must be satisfied by non- 
trivial plane-wave solutions. For the CPT-even viola- 
tions, we extend the widely used k matrices [11] of the 
minimal SME to general k operators (32) that are rele- 
vant for studies at arbitrary d. 

The unconventional properties of the Lorentz-violating 
eigenmodes can be characterized in terms of birefrin- 
gence, dispersion, and anisotropy. All CPT-odd oper- 
ators lead to birefringence. We study the conditions for 
leading-order birefringence from CPT-even operators via 
a Weyl decomposition of the constitutive tensor. For 
this case, we conjecture that nonbirefringent terms are 
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uniquely associated with the non-Weyl part of the con- 
stitutive tensor when the decomposition is expressed in 
terms of a suitable effective metric. For arbitrary Lorentz 
violation, we introduce a duality symmetry determined 
by the effective metric and argue that birefringence is a 
consequence of the breaking of this symmetry. 

In Sec. Ill, we obtain a complete classification of the 
coefficients for Lorentz violation at arbitrary d, using an 
SO (3) decomposition in terms of spin- weighted spheri- 
cal harmonics. A review of the spin-weighted spherical 
harmonics and derivations of some useful mathematical 
results arc provided in Appendix A, along with a dis- 
cussion of the relationships to angular momentum, he- 
licity, and parity. The SO (3) decomposition reveals a 
total of nine independent sets of spherical coefficients for 
Lorentz violation that control birefringence, dispersion, 
and anisotropy in the photon propagator. Table XVII 
lists these spherical coefficients and displays the ranges 
of their indices and their counting. The labels and indices 
on a given spherical coefficient identify the key proper- 
ties of the corresponding Lorentz- violating operator, and 
their interpretation is summarized in the paragraph con- 
taining Eq. (45). More detailed properties of these nine 
sets are compiled in Tables I, III, IV, and V. 

Table XVII reveals the existence of several classes of 
operators that are absent in the minimal SME, which 
is restricted to d = 3 and d = A. The CPT-odd case 
at = 3 involves two sets containing four independent 

spherical coefficients, while the CPT-even case at d = 4 
involves three sets with 20 coefficients, one of which is 
an unobservable constant. However, the CPT-odd case 
at d = 5 has 35 coefhcients distributed among three sets 

rather than two. The additional set {k^py^jj^ controls 
iJ-parity CPT-odd effects, which are absent in the mini- 
mal SME. Similarly, the CPT-even case at d = 6 has 126 
coefficients distributed among six sets rather than three. 

The additional sets [I^f (^f )nf2^ (^f )S 

ern iJ-parity spin-one operators and also spin-two oper- 
ators with both E- and _B-type parities, all associated 
with qualitatively new effects. 

For applications to observation and experiment, it is 

valuable to define various limiting cases of the general 
theory. Sec. IV presents several of these limits, while 
Table XVIII summarizes their specific content in terms 
of spherical coefficients for Lorentz violation. The most 
widely studied limit to date is the minimal SME, which 
involves a total of 23 independent nontrivial spherical co- 
efficients. The explicit connections between the spherical 
coefficients and the usual cartesian coefficients for the 
minimal SME is provided via Tables VI, VII, and VIII. 

Another useful limit involves the specification of a pre- 
ferred frame in which all Lorentz violation is isotropic. 
The corresponding isotropic or 'fried-chicken' models arc 
discussed in Sec. IV B. The isotropic requirement elim- 
inates all but three sets of spherical coefficients and re- 
duces the growth of coefficient numbers to be linear in- 
stead of cubic at large d. For CPT-odd isotropic oper- 



ators with d, = 3,5,7... there are only 1, 2, 3 . . . types 
of Lorentz violation, all of which are bircfringcnt. Simi- 
larly, for CPT-even isotropic operators with (i = 4, 6, 8 . . . 
only 2,4,6... coefficients arise, of which 0,1,2... are 
associated with leading-order birefringence. Note that 
isotropic CPT-even birefringence is a physical feature 
only for d > 6. 

A third useful limiting subset of the general theory is 
obtained by restricting attention to Lorentz- violating op- 
erators that are nonbirefringent and also are nondisper- 
sive in the vacuum at leading order. Sec. IV D constructs 
the corresponding camouflage models. These models in- 
volve effects that are challenging to detect in astrophys- 
ical studies of birefringence and dispersion. They are 
described by the single subset of coefficients (<?^'')|i"^^ 
appearing at even d and hence associated with CPT-even 
Lorentz violation. Table IX summarizes some properties 
of these coefficients. For d = 4 only a Lorentz-invariant 
trace appears, so for most purposes it suffices to take 
d > 4. The camouflage coefficients then govern effects 
outside the minimal SME. For operators of larger mass 
dimension d = 6, 8, 10 . . . there arc 10, 35, 84 . . . indepen- 
dent effects, with the number of independent coefficients 
growing rapidly as the cube of d for large d. Also, for 
d = 4,6,8... there are 1,2,3... isotropic camouflage 
coefficients, which therefore also belong to the general 
isotropic model. 

The extreme sensitivities to Lorentz violation available 
via studies of birefringence and dispersion of astrophysi- 
cal sources provides motivation for a further refinement 
in the classification of spherical coefficients, based on sep- 
arating those coefficients that affect the propagation of 
light in the vacuum from all others. We refer to the 
former as vacuum coefficients and to the complement 
as vacuum-orthogonal coefficients, and we identify the 
latter by a negation diacritic ^. This classification is 
summarized in Table XVIII. The vacuum-orthogonal 
coefficients appear only for d > 5, so they represent 
qualitatively new effects that are absent in the minimal 
SME. Moreover, the numbers of vacuum and vacuum- 
orthogonal coefficients grow as and d"^ for large d, re- 
spectively, so the vacuum-orthogonal coefficients repre- 
sent most of the coefficient space for large d. For exam- 
ple, the vacuum coefficents govern all the physical effects 
in the minimal SME for which d = 3 and d = 4, but for 
d = 5 they span only 16 of the 36 possibilities and for 
d = 6 only 67 of 126. 

The vacuum coefficients are constructed in Sec. IV C. 
At leading order, they are identified by requiring that 
the radiation fields are plane waves. This restricts at- 



tention to four sets of coefficients c|jj^- 



(E)jm^ "'(B) jm^ 

are associated with CPT- 
even Lorentz violation that is nonbirefringent at leading 



^(Vbm- coefficients cj^}^.^ 



order but dispersive for d > 4. The coefficients k 



and k 



id) 

(B)jrn 
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control CPT-even bircfringcnt effects that are 



also dispersive for d > 4. Only k 
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TABLE XVII: Summary of spherical coefficients for Lorentz-violating operators of arbitrary mass dimension. The first column 
specifies the CPT property and the total number of independent operators at each d. The corresponding spherical coefficient 
sets are listed in the second column. The remainder of the table provides the allowed ranges for the indices d, n, j and the 
number of independent components for each coefficient set. 
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TABLE XVIII: Summary of limiting cases. The first column specifies the limit, while the corresponding spherical coefficient 
sets arc listed in the second column. The remainder of the table provides the allowed ranges for the indices d, n, j and the 
number of independent components for each coefficient set. 



49 



CPT-odd effects, which arc also bircfringcnt and disper- 
sive. For any given d, the numbers of each type of vacuum 
coeflicients are roughly the same. For example, at d = 6 
there are 25 coefficients in c^^j^.^ and 21 each in k'"^E)jm 
and k'g)^^. 

The vacuum-orthogonal models are introduced in Sec. 
IV E. There are nine sets of vacuum-orthogonal coeffi- 
cients, of which four arc identical to and five arc reduced 
subsets of the nine sets in the general analysis. One of 
the latter consists of the camouflage coefficients. Some 
properties of the remainder are provided in Tables X, 
XI, and XII. Except for a single constant scale factor 
in (1 = 4, the vacuum-orthogonal coefficients appear only 
for d > 4. They govern Lorentz-violating operators that 
produce no leading-order birefringence or dispersion in 
vacuum propagation, although bircfringcnt or dispersive 
effects can appear in other physical contexts. For exam- 
ple, all 20 vacuum-orthogonal coefficients for d = 5 and 
49 of the 59 for d = 6 are associated with bircfringcnt 
and dispersive effects in suitable circumstances. The re- 
maining 10 for d = 6 are camouflage coefficients, which 
govern nonbirefringent effects. 

The remainder of the paper is concerned with applica- 
tions of the theory to observation and experiment. Re- 
sults of measurements are conventionally reported in the 
canonical Sun-centered frame, so for some applications it 
is necessary to perform transformations relating spheri- 
cal coefficients in the Sun-centered frame to a laboratory 
or other frame. These transformations are provided ex- 
plicitly in Sec. V for the case where the boost component 
is negligible. 

Section VI describes some applications in the astro- 
physical context. We focus on astrophysical studies 
of vacuum birefringence and dispersion. Birefringence 
studies offer extreme sensitivity to many vacuum coef- 
ficients, while dispersion tests access the remainder at 
lesser but nonetheless impressive sensitivities. Disper- 
sion constraints are discussed in Sec. VIA. We obtain 
expressions applicable to both isotropic and anisotropic 
dispersion, and we use the recent measurements of GRB 
080916C made by the Fermi Observatory to obtain new 
constraints on Lorentz violation involving operators of 
mass dimension six and eight. A summary of existing 
constraints from astrophysical dispersion tests is pro- 
vided in Table XIII. 

Birefringence constraints arc considered in Sec. VI B. 
The general theory is outlined, and the net rotation 
induced by arbitrary spherical coefficients is described 
quantitatively in terms of Stokes parameters. We use po- 
larimetric data from the gamma-ray bursts GRB 930131 
and GRB 960924 to set tight constraints on spherical 
coefficients associated with Lorentz-violating operators 
of mass dimensions four through nine. A summary of 
current limits from GRB polarimetry is provided in Ta- 
ble XIV. We also consider constraints obtained from 
polarimetric studies of the CMB. The maximal photon- 
propagation distances and the photon frequencies make 



the CMB particularly well-suited for measurements of 
spherical coefficients with d = 3. A survey is per- 
formed to categorize the mixing of polarizations in the 
CMB induced by Lorentz violation. Table XV compiles 
some existing limits from the CMB on both isotropic and 
anisotropic Lorentz violation. 

Section VII discusses applications in the laboratory 
context, focusing on the use of resonant cavities to search 
for Lorentz violation. These systems offer sensitivities 
to spherical coefficients that are challenging to access 
in studies of vacuum birefringence and vacuum disper- 
sion. A general theoretical procedure for deriving the 
fractional frequency shift in a cavity is presented in Sec. 
VII A. In the following subsections, the results are ap- 
plied to nonbirefringent Lorentz violation and in par- 
ticular to the camouflage coefficients, for which some 
relevant experiment-dependent factors are explicitly ob- 
tained. For the specific case of a circular-cylindrical cav- 
ity, we provide in Eq. (184) the fractional frequency shift 
including the explicit time dependence. 

The analysis in this work demonstrates that a com- 
prehensive search for Lorentz violation is best performed 
with multiple types of measurements. The most sensi- 
tive tests use astrophysical birefringence, which provides 
access to some vacuum coefficients. Astrophysical disper- 
sion offers high sensitivity to the remaining vacuum co- 
efficients. However, the bulk of effects involves vacuum- 
orthogonal coefficients. Direct sensitivity to these re- 
quires non-vacuum boundary conditions and hence lab- 
oratory studies. Disentangling the various birefringence, 
dispersion, and anisotropy effects requires a variety of 
laboratory experiments involving different boundary con- 
ditions and different media, and also the exploitation of 
signals from different rotations and boosts. Even if at- 
tention is limited to coefficients with comparatively low 
values of d, considerable room remains for investigations 
via both astrophysical observations and laboratory ex- 
periments. 
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APPENDIX A: SPHERICAL HARMONICS 

Angular-momentum eigenstates are irreducible repre- 
sentations of rotations, so tensors in three dimensions 
can be decomposed into components with definite or- 
bital angular momentum and spin. For spin-zero scalars, 
the spherical harmonics Yjm provide a basis for this de- 
composition. For nontrivial tensors, more general tensor 
spherical harmonics are needed to incorporate the spin. 
Among the most widely used are the spin- weighted spher- 
ical harmonics [57, 58], denoted sYjm- In this Appendix, 
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we briefly review the definitions and basic features of 
these spin-weighted spherical harmonics. We also obtain 
some mathematical properties that are used in Sec. III. 



1. Spin- weighted spherical harmonics 

A key concept underlying the spin-weighted spherical 
harmonics is the notion of spin weight. To introduce 
this idea, consider the problem of characterizing radia- 
tion propagating inward toward the Earth from a distant 
point source in the sky. The electric-field vector E is 
oriented perpendicular to the line of sight, so it lies in 
the tangent space of a sphere surrounding the Earth. In 
spherical polar coordinates, the angular components of 
E are Eg and E^. However, the alternative components 
E± Qc Eg ^ iEcj) can be considered instead. These have 
the advantage of transforming elegantly as E± e^^^ E± 
under a rotation of the local coordinates by an angle 5 
about the line of sight. The irreducible combinations E± 
are said to be spin- weighted functions of spin weight ±1. 
More generally, a function «/ is said to have spin weight 
s if it transforms according to gf — > e"'*'^^/ under a local 
rotation about the line of sight. 

The generator of rotations is the the angular- 
momentum operator J . Denoting the radial unit vector 
as n, the generator of local rotations about the line of 
sight is the operator h ■ J. This is the helicity opera- 
tor with respect to the line of sight. We see that the spin 
weight can be understood as the eigenvalue of the helicity 
operator. For instance, the irreducible combinations E± 
in the above example are components of definite helicity 
with respect to the line of sight. In this particular exam- 
ple, the light propagates in the 3-momcntum direction 
p = —h, so the helicity operator n • J is equivalent up to 
a sign to the helicity operator with respect to the momen- 
tum, p ■ J. However, this eqiiivalence fails for non-radial 
propagation. So in general there are two basic options 
for defining spin weight, as the eigenvalue up to a sign 
of either n ■ J or p ■ J. The choice of definition can be 
made based on suitability for the problem at hand. 

The helicity operator p j or h - J commutes with both 
the squared total angular momentum and the com- 
ponent Jz, where by convention we choose the projection 
axis to be the z direction. It is therefore possible to intro- 
duce simultaneous eigenfunctions for all three operators. 
We show in the next subsection that these simultaneous 
eigenfunctions are the spin-weighted spherical harmon- 
ics sYjm- However, typical discussions [57-59] of sYjm 
are based on raising and lowering operators for the spin 
weight, which connect different harmonics. In this sub- 
section, we define these operators and provide explicit 
expressions for gYjm- 

Acting on a function of spin weight s, the raising 
and lowering operators are given by 



The raising operator S acts on a function of spin weight 
s to yield a function of spin weight ,s + 1 . Similarly, the 
lowering operator 9 decreases spin weight by one. The 
spin- weighted spherical harmonics sYjm of spin weight s 
can be generated from the usual scalar spherical harmon- 
ics Yjm = oYjm by raising or lowering the spin weight \s\ 
times: 



Y 



Y {j+s)\ O-'jrri) O ^ 6 J, 



(A2) 



As discussed in the following subsection, the index j cor- 
responds to the eigenvalue of the squared total angular 
momentum = j{j + 1) and m to the eigenvalue of J^. 

The functions sYjm are nonvanishing for index values 
I'm] < i) as usual. Since the helicity is limited by j, they 
are also nonvanishing for j > \s\. They satisfy a relation 
analogous to that obeyed by the usual scalar spherical 
harmonics, 



Y Y 



/(2ii + l)(2j2±l) 
V 47r(2j3-|-l) 

X 0li2(-Sl)(-S2)b3(-S3)) 

X 0li2mim2|i3TO3) ssYj^ms, (A3) 

where the symbols {jij2'niifn2\j3'mz) represent the 
Clebsch-Gordan coefficients. This implies orthogonality 
of harmonics of equal spin weight, 



^Y-rnin) sYj,^.{n) sin 



(A4) 



Note, however, that this orthogonality does not extend 
to sYjm of different spin weight. The harmonics sYjm 
also satisfy the completeness relations 



^sY*,n{n) sYjm{n ) = d{n - n ) . 



jm 



Here, we adopt phase conventions ensuring 



XT* 



(_l\s+m Y-, . 



under complex conjugation and 

sYjm{-n) = {-ly -sYjm{fi) 



(A5) 



(A6) 



(A7) 



under parity. 

Although Eq. (A2) can be used to generate the spin- 
weighted spherical harmonics, in practice it is often easier 
to use the explicit expression 



sYjm{e,ct>) 



{j+my.{j-my. 

i-K (3+sy.{j-s)\ 



gim,^gij^2j| 



/ = - sin*** e{d0 ± i esc 63^) sin^' 6 «/. (Al) 



where ( ™ ) denotes the binomial coefficients. However, 
this expression can be numerically troublesome for high 
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values of j. One strategy for numerical applictions is 
to use Eq. (A8) to generate the low-j harmonics but in- 
stead to use recursion relations derived from Eq. (A3) 
to extract the higher-j ones. We find that two re- 
cursion relations are useful for this purpose. Taking 
si — 0,ji = l,mi = 1, and m2 = ±j2 in Eq. (A3), 
we obtain the recursion 

«^i(±i) = TV'^^e±'^sin^,y(,._i)(±,.^i), (A9) 

which relates harmonics with j = \m\. Also, taking 
si = 0, ji = l,mi = in Eq. (A3) leads to the recursive 
formula 



/((2j)''-l) 



0-l)^(2j-l)(2j-3) sm-2)m 



which permits the calculation 
lower-j ones with the same m 
determining numerical values 
harmonics then involves first 
ues for j = \s\, followed by 
harmonics with j = |m|, and 
fill in all remaining values up 



(AlO) 

of higher-j harmonics from 
A practical procedure for 
of spin-weighted spherical 
using Eq. (A8) to find val- 
using Eq. (A9) to obtain 
finally using Eq. (AlO) to 
to the desired maximum j. 



2. Covariant angular momentum 

In this subsection, we further explore the relation- 
ships between angular momentum, helicity, and the spin- 
weighted spherical harmonics. For definiteness, we con- 
sider helicity with respect to the momentum direction p, 
but the following discussion remains valid if n is substi- 
tuted for p throughout. Note also that we adopt a metric 
with positive signature when working with 3-dimensional 
tensors. 

We introduce standard angles in spherical polar c;oor- 
dinates, so that the momentum direction can be written 
as 

p — sin 9 cos <j>ex + sin 6 sin <t>ey + cos 6 e^, (All) 

where e^, e-y, form the cartesian basis vectors. The 
orthonormal basis vectors in the spherical polar coordi- 
nates are taken as 



er = e=p, e0 = e , = 



(A12) 



where ee and e,j, are the coordinate unit vectors asso- 
ciated with the spherical coordinates 6 and (j). We also 

define complex helicity-basis vectors 

er = e^=p, e±= = ^{e0±ie^), (A13) 

where an explicit choice of phase has been made. Other 
phase choices would change some of the relations below. 



With this definition, a rotation of the local coordinates 
about p by an angle 5 generates a phase shift in the 
helicity-basis vectors, e± — > e^^^e±. 

Using the helicity basis, we can readily decompose 
any tensor into components of definite helicity and spin 
weight. For example, a vector V has 0-helicity compo- 
nent Vr = V-i-r and ±l-helicity components = V^: = 
V-e^, the latter corresponding to spin-weight = =f1. The 
spin weight and helicity of any tensor component is read- 
ily obtained by simple counting, since each =F contravari- 
ant or ± covariant index adds ±1 to the spin weight. 

To construct operators acting to raise or lower the 
spin weight, it is convenient first to introduce directional 
derivatives in momentum space with respect to arbitrary 
basis vectors Co, according to 



The metric is 



da=ea- d, daP = Ba- 



9ab — ' ^b-) 



(AM) 



(A15) 



and we define covariant directional derivatives Va with 
connection 



(A16) 



as usual. For example, acting with Va on the components 

of a vector V gives 

VaV' = daV'' + r\^V'\ VaVb = daVb-r^,,Vc. (A17) 

In the helicity basis {e+, e^, e_}, we find the metric takes 
the explicit form 



/ 1 ^ 

9ab = 9""= 10 , 

\l oj 

while the nonzero connection elements are 



(A18) 



±r 



±± 



= -F= 



■P 



= {V2p ta,ne)-\ 



(A19) 



where p = \p\. The point is that the derivative V+ op- 
erating on an arbitrary tensor acts to create tensor com- 
ponents whose spin weight is increased by one, while V_ 
lowers the spin weight by one. 

Our goal is a decomposition in angular momentum, so 
we seek a covariant description of the angular-momentum 
operators. In momentum space, the orbital angular mo- 
mentum L has components 



-le pbVc 



and obeys the algebra 

[La, Lb] = —iSabcL''. 



(A20) 



(A21) 



Here, e"^'^ is the totally antisymmetric tensor, which in 
the helicity basis satisfies 



-l-r— 



(A22) 
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In particular, we have 

= 0, L± = ±pV±, (A23) 

which implies that orbital angular momentum is perpen- 
dicular to p, as expected. 

To complete the description of angular momentum, it 
is useful to define a covariant spin operator S with com- 
ponents Sa- Acting on a vector Va, we define 

5„n = -S\^V,, SaV' = S^^V, (A24) 

where 

S',, = is\,. (A25) 

This provides a covariant formulation of the general spin 
relation 

(Vi . S)V2 = iVi X = {Vire\Sa{V2)b). (A26) 

The spin operator acting on more general tensors obeys 
rules like those of the connection. For example, we obtain 

SaT\ = S\,T'', - S\,T\. (A27) 

Since the Sa are covariant operators rather than matrix 
operators, they obey slightly modified commutation re- 
lations 

[Sa,Sb] = -ieabcS''. (A28) 

With the above definitions, we can introduce a covari- 
ant operator for the total angular momentum as 

J„ = L„ + 5„. (A29) 

Like the covariant derivative, this operator has the ad- 
vantage of maintaining the explicit tensor nature of ex- 
pressions when operating on tensor components. It obeys 
the modified commutation relations 

[Ja, Jfc] = -zeafccJ". (A30) 

An interpretation of these operators is as follows. The 
radial angular momentum = Sr is the helicity, which 
is the negative of the spin weight. The operator J+ raises 
the helicity and J+ raises the spin weight, while J~ low- 
ers the helicity and J_ lowers the spin weight. These 
operators are equivalent to the raising and lowering op- 
erators Eq. (Al) via the correspondence 

J+ = -g/V2, J_ = 8/72. (A31) 

The spin- weighted spherical harmonics gYjm are eigen- 
functions of the commuting operators 

-J^ = s, J2=j(j + 1), J,=m. (A32) 

The ladder operators J± can be shown to commute with 
and Jz- In fact, J± commutes with V • J for any 



covariantly constant vector V, so that VaH = 0. Also, 
the commutator 

[J±,-Jr\ = ie±r^J^ = tJ±, (A33) 

implies that J+ and J_ raise and lower the spin weight, 
respectively, while leaving j and m unaltered. The lad- 
der nature of the J± can thus be seen either from the 
differential geometry or from the operator algebra. 

Since the spin-weight ladder operators J± are analo- 
gous to the conventional ladder operators that raise and 
lower the eigenvalue of the z component of the angu- 
lar momentum, many of the usual techniques and results 
apply. For example, the standard normalization scheme 
gives 

J± sYjm = -^^{jij + l)-s{8±l)) s±lYjm, (A34) 

where the phase is chosen to match the conventional def- 
initions when expressed in terms of S and 8. Repeated 
application gives the useful expressions 

oY,^ = i-iy2l^\/^^f^^{J^p (A35) 

and 

±Hy,-„ = i-iy2\^\/^^^^{J^)\^\ ,Yjm, (A36) 

which are equivalent to Eq. (A2). 

In practical terms, the above formalism leads to the 
simple results 

V. = d/dp, 

V± = ±(J± - S±)/p, (A37) 

which will be used extensively in Section III. Note, how- 
ever, that some care is needed when rising the results 
(A37) since they are non-tensorial relations, valid only 
in the helicity basis. For example, the ± indices on the 
right-hand side are noncovariant once applied, so further 
manipulations with covariant operators such as Va, Sa, 
J a are inappropriate. 

As a simple illustration of the above formalism, con- 
sider the covariant momentum-space laplacian acting on 
a scalar iIj: 

Va W = (Vr-V + V+V+ + V-V")^ 

= {w + W+ - ^+)^^ - W- - 

(A38) 

Using the results 

{J+J+ + J-J-)ij = J^ij, (A39) 
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we can reexpress VaV"V ^ follows: 
v„w = + ^(^+v+ + V.) - i(J_ V- - V.)) V 

i. dp'^ ^ p dp p'^ (A40) 

We thereby recover the familiar expression for VaV^V 
involving the angular momentum. 

It has been noted in the literature that the operators 
3 and 6 are equivalent to covariant derivatives in the 

two-dimensional tangent space of the sphere [57]. In our 
present language, this can be made apparent by defining 
a new covariant derivative as 

Va = Va - iSarbS^/p, (A41) 

which in the helicity basis gives nonzero connection ele- 
ments 

f±±± = -r=^^± = (\/2 p tan6i)-^ (A42) 

The derivative Vq corresponds to the projection of the 
usual three-dimensional geometry onto the embedded 2- 
sphere. In terms of this derivative, we can express the 
ladder operators as 

J± = ±V±/p. (A43) 

This relates our three-dimensional picture to the two- 
dimensional one. 



3. Parity 

In a typical application of the above formalism, a ten- 
sor is decomposed by considering its components in the 
helicity basis, which are spin-weighted functions, and ex- 
panding them using spin- weighted spherical harmonics. 
In many cases, it is useful to decompose fmther the re- 
sults in terms of parity properties [57]. This subsection 
contains a brief summary of the latter procedure, along 
with some figures providing insight into the structure of 
the resulting modes. 

Let t represent an arbitrary tensor component with 
spin weight s. Its expansion then takes the form 

t{p) = J2t3m sYjmip)- (A44) 

jm 

By exchanging all the -|- and — indices on t, we obtain 
another tensor component of spin weight —s. Denoting 
this component as i, we can expand it as 

= X]*J>" -sYjmip)- (A45) 

jm 

The two components t and t are parity conjugates, inter- 
changing under parity according to 

t{p)^{-iyt{-p). (A46) 



In terms of spherical coefficients, the parity transforma- 
tion gives 

tjm ^ i-iy+' ijm. (A47) 
It is then useful to define 

tjm = i-^nEjm - iBjm). (A48) 

The result is a splitting of t and t into modes with so- 
called electric-type parity or i?-typc parity, 

Ejm ^ {-lyEjm, (A49) 

and ones with magnetic-type parity or S-type parity, 

Bjm ^ i-iy+^Bjm, (A50) 

where the nomenclature is borrowed from radiation the- 
ory. 

As a simple example, consider a scalar function 

S{p) = Y^Ejrr^oYjm{p). (A51) 

jm 

Under parity, S transforms according to 

S{p) - S{-p), (A52) 

which in conjunction with Eq. (A7) confirms that scalar 
functions contain only i?-type components. Figure 9 il- 
lustrates the angular distribution obtained by considering 
each mode Eqo, Eiq, Eh, E20, E21, E22 of a real scalar 
field in turn. For example, the real scalar field S{Eii) 
associated with real En is 

S{En) = EnioYn - on(-i))- (A53) 

For each of the six images in Fig. 9, the three arrows per- 
pendicular to the spherical surface represent rectangular 
right-handed coordinate axes with z vertical. The solid 
disks represent positive values of the scalar, while the 
rings represent negative ones. The radius of each disk or 
ring is proportional to the magnitude of the scalar at that 
point. The images reveal the angular and parity symme- 
tries of each scalar mode. Consider, for example, the Eio 
scalar mode. The figure shows that this mode has ex- 
tremal magnitudes at the poles and equator and is sym- 
metric under rotations about z, all of which are features 
of a distribution with j = 1 and m = 0. Also, the sign 
of the mode changes on interchanging any two antipodal 
points, matching its negative parity. Note that the distri- 
butions are plotted for real amplitudes Ejm- Inclusion of 
a phase in an amplitude rotates the corresponding figure 
about z. 

An example with both E- and B-type content is pro- 
vided by the components 

V± = ^i±Ejm + iBjm) ±lYjm (A54) 
jm 
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FIG. 9: Angular distributions for electric-type scalar compo- 
nents. 



of a vector field V in the tangent bundle of the sphere. 
Figure 10 displays the angular distributions for a real 

vector field for each separate mode Eqq, -Eio, En, E20, 
E21, E22- For example, the components of the real vector 
field V(£'ii) associated with real En are 

V+iEn) ^ EniiYn - iY^_,)), 

V-{En) = -Eni-iYn - -i^c-i)). (A55) 

In the figure, the vector tangent to a sphere at a given 
point represents the corresponding mode. The angular 
and parity symmetries of each vector mode can be seen 
by inspection. For instance, the Eio vector mode has ex- 
tremal magnitudes at the poles and equator and is sym- 
metric under rotations about z, as expected for a mode 
with j = 1 and m = 0. Also, inspection reveals that ap- 
plying the parity operation reverses the fiow of vectors on 
the sphere, as is appropriate for a mode with E-type par- 
ity and j = 1. Note that the E'-type modes are curl free 
while the -B-type modes are divergence free, matching the 
usual properties in electrodynamics. Note also that for 
given values of j and m the E- and i?-mode field lines are 
perpendicular everywhere, reflecting the orthogonality of 
the modes. 

For a symmetric traceless 2-tensor, we can consider the 
components 

T±± = J2(^jrn ± iBjm) ±2Yjm- (A56) 
jm 




FIG. 10: Angular distributions for electric- and magnetic-type 
vector components. 



Figure 11 provides representations of each the six modes 
E20, E21, E22, B20, B21, B22 for a real tensor. The 
figure takes advantage of the spectral theorem applied to 
the decomposition of a real symmetric traceless tensor T 
that lies in in the tangent bundle of the sphere and has 
nonzero components T±±. In terms of its two orthogonal 
eigenvectors v and w lying tangent to the sphere, we can 
write 

T = vi^v-wi^w, (A57) 

where the eigenvalues of T are \v\^ and — |wp = — I'up- 

The tensor T can therefore be represented at each point 
on the sphere by the two vectors v and w. We denote 
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FIG. 11: Angular distributions for electric- and magnetic-type 
tensor components. 



the vector v hy a thick hne and the vector by a thin 
one, noting that the vector orientations are unspecified 



by Eq. (A57) and so the vectors arc best represented 
by unoricntcd hne segments. The angular and parity 
symmetries of each tensor mode are visible in the figure. 
For example, the extremal magnitudes of the E20 tensor 
mode occur at the poles and equator, and the mode is 
symmetric under rotations about z, as expected for a 
mode with j = 2 and m = 0. The symmetry under 
interchange of two antipodal points is a consequence of 
the positive parity of this mode. The orthogonality of 
the E and B modes for given values of j and m implies 
that each E-B pair of modes is related by interchanging 
a plus with a cross at each point. As an aside, we remark 
that this visualization has a parallel in general relativity, 
where the orthogonality of the plus and cross modes for 
gravitational radiation also arises from a symmetric real 
2-tensor tangential to the direction of propagation. 

We can also consider various types of pseudotensors. 
All these acquire an additional sign under parity, which 
implies that the roles of the E- and B-type coefficients are 
interchanged in the parity decomposition. For example, 
pseudoscalars contain only B-type components. Also, 
the parity decomposition of a pseudovector V' leads to 
components with spin weight ±1 of the form 

Vi = Y.{±Bjm + iEjm) ±lYjm (A58) 

jm 

instead. Similarly for a pseudotensor T', we have 

T±± = Y^iBjm ± iEjm) ±2Yjm- (A59) 
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